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1. SUMM ARY 

This report describes an investigation that was conducted to select and develop an 
improved analytical method for predicting the stability and control characteristics of elastic 
airplanes. * 

The investigation included consideration of the free-flight conditions of large airplanes 
in their “clean” configuration. Landing, takeoff, ground effects, stability augmentation, 
and control surface movements were not considered, nor were prediction of air- 
plane performance, flutter, or structural loads. The study was confined to a fligL' 
envelope extending from low subsonic speeds to Mach 5 and from sea level to 30 000 meters 
(93 360 ft) altitude. Within this envelope, rectilinear and curvilinear reference flight paths 
were taken and analyses made of arbitrary, large, and small perturbations of airplane motion 
about the reference flight paths. Airplane structural motions of dynamic-elastic and 
quasi-static-elastic characters were included, and their effects on stability determined. 


The approach taken in the investigation was to develop in order: 

• equations of motion 

• stability criteria 

• stability derivative prediction methods 

• stability characteristics prediction methods 

A major result obtained in the investigation has been the unification and development of 
aerodynamic, structural, and dynamic technologies into an overall plan for calculating 
elastic airplane stability characteristics. The plan involves a computing program system 




There are two essential features to this plan: 

1 . The system elements and the total system are well suited to computer programming. 

2. Provision is made for introducing empirical data (experimental measurements, 
handbook results, etc.) into the system. 

Other results and conclusions concerning the analysis methods and their application to 
large, flexible airplanes of the SST and 707-320B type are listed below. 

1 . The lumped parameter concept is the most practical way to represent an airplane for 
analysis. This involves paneling the airplane so that the aerodynamic and structural 
relationships among panels can be expressed in terms of influence coefficients. A 
geometry definition program is desirable to mechanize and help guide the paneling. 

2. Lifting surface theory gives good results and is the most suitable method for 
determining the aerodynamic influence coefficients. 

3. For most airplane configurations, equivalent beam structural models are adequate 
for determining structural influence coefficients. For very low aspect ratio wings or 
for increased accuracy, however, influence coefficients from large, sophisticated, 
finite element structural programs should be used. 

4. The stability criteria that apply to rigid airplanes apply directly tc elastic airplanes. 

5 . Stability characteristics can be influenced strongly by inaccuracies in estimating che 

rigid stability derivatives. These effects on design decisions are of the same order of 
importance a* 1 quasi-static-elastic and dynamic-elastic effects. 

6. The effect of airplane loading conditions (mass distribution) on an elastic airplane’s 
stability characteristics can be relatively large, sometimes reversing the sign of the 
parameters that represent flexibility effects. 

7. The evaluations of control effectiveness and control surface angles required to trim 
are influenced strongly by flexibility; however, since viscous and nonlinear 
aerodynamic effects are also important in those evaluations, the applicability of 
lifting surface theory based on potential flow as a prediction technique is limited. 

8. Airplane stability characteristics are, in general, moderately affected by flexibility. 
The quasi-static-elastic formulation is usually adequate for most prediction tasks, 
dynamic effects being generally modest. However, some dynamic check casi , 
should be ran, since there is no assurance that dynamic effects will be small for any 
given configuration. 

It is recommended that further development of prediction methods follow the plan* 
mentioned above. Lifting surface theory and beam analysis should be used to determine 
aerodynamic and structural influence coefficients and to calculate stability derivatives and 
stability characteristics. For dynamic stability investigations, if the center-of-gravity 
disturbance is characterized by the small perturbation equations of motion, the character- 
istic equation rooting technique may be used. The stability of large perturbations is best 
evaluated by time history calculations. 

*Mechanization of this plan is currently under way under the direction of personnel of the 
Non-Steady Phenomena Branch, NASA-Ames. 



2. INTRODUCTION 


The stability of an airplane is its tendency to persist in a particular reference motion 
(for example, steady, level flight) when it has been disturbed from that motion. Primary 
factors affecting stability are the changes in the aerodynamic forces and moments acting on 
the airplane that occur with changes in the airplane’s motion and orientation. They are 
expressed as stability derivatives evaluated at the reference motion condition, such as the 
change in airplane lift coefficient with change in angle of attack, BCk/do-i^f = Cj^. 

In the past, the effects of structural flexi ity on airplane stability were accounted for 
by modifying the stability derivatives. For example, a change in airplane lift due to a change 
in angle of attack might produce twisting of the wing or bending of the fuselage, resulting in 
a different lift change for a particular angle-of-attack change than would occur had the 
airplane been rigid. The rigid airplane stability derivatives, were then corrected to include 
these quasi-static-elastic effects. 

The effects of structural dynamics have generally been a consideration only in flutter 
predictions. The center of gravity of the airplane was (and is) considered to be in a state of 
steady, level flight while the structure is disturbed. For example, neutral flutter stability is a 
constant-amplitude, oscillatory structural motion resulting from a disturbance; however, this 
motion is regarded as having no effect on the steady, level motion of the center of gravity of 
the airplane. This is a satisfactory representation of the motion provided the frequency of 
the structural motion is well separated from the natural frequencies of the overall motion of 
the airplane, e.g., its short period longitudinal natural frequency. In this case the two 
motions are not sufficiently coupled for an exchange of energy between them. 

Certain aerodynamic and structural approximations have also been used. Aerodynamic 
surfaces -wings and tails— may have sufficiently large aspect ratios that lifting line theory 
may be used to predict aerodynamic loads and the structure may be treated essentially as an 
assemblage of beams. These and other approximations were satisfactory and led to relatively 
simple methods for predicting the influence of structural flexibility on airplane stability. 

However, the advent of large airplanes operating in the transonic and supersonic flight 
regimes has led to configurations for which some of the acceptable approximations of the 
past are of questionable validity or are obviously invalid. The frequencies of the structural 
motion have been sufficiently reduced by the increase in both airplane flexibility and cruise 
dynamic pressure that coupling with overall motion of the airplane is attendant. The use of 
low aspect ratio aerodynamic surfaces invalidates the lifting line aerodynamic approxima- 
tion and reduces the applicability of equivalent beam structural models. Different methods 
of analysis must be introduced. 



i Bisplinghoff and Ashley (ref. 1, chapter 9) have laid the ground work for new, less 
restrictive methods by presenting equations of motion that have the required degree of 
generality. Milne (ref. 2) also presents the development of more general equations of motion 
that integrate conventional stability and aeroelastic methods. Milne’s: rather extensive work 
also includes the application of the equations to the problem of slender airplane trim state 
and longitudinal stability. Most major airplane companies are also involved in developing 
new, less restrictive methods. However, results of these studies are usually not available in 
the open literature. 

With this background in mind, an investigation was conducted for the purpose of 
developing an improved method for predicting the stability and control characteristics of an 
elastic airplane. Objectives of the study were to: 

1. Develop equations of motion that have sufficient generality to handle large, flexible 
airplane stability and control problems, including dynamic-elastic effects. 

2. List and evaluate assumptions and restrictions introduced and determine how they 
may influence the prediction of airplane stability and control characteristics. 

3. Develop stability critt" ; " applicable to flexible airplanes. 

4. Develop improved methods for predicting elastic airplane stability derivatives using 
current and/or improved aerodynamic and structural techniques. 

5. Develop an improved approach for predicting elastic airplane stability characteristics 
using the results of 1 through 4 above. 

6. Document the results of the study in a precise, understandable form. 

The work was accomplished under the technical direction of the Non-Steady 
Phenomena Branch, Ames Research Center, Moffett Field, California. Members of the 
Aerodynamics and Structures Staff of the Commercial Airplane Division of The Boeing 
Company at Renton, Washington conducted the investigation as a joint effort. Frequent 
coordination meetings and reviews with Ames representatives were held during the course of 
the contract. 

The scope of the investigation included consideration of large, flexible 
operating in the flight envelope of fig. 1. Only “clean” configurations were studied; landing, 
takeoff, ground effects, stability augmentation, and control surface movements were not 
considered. Rectilinear and curvilinear reference flight paths were assumed, and arbitrary, 
large, and small airplane motion perturbations about the reference flight paths investigated. 
Structural motions of dynamic-elastic, quasi-static-elastic, and rigid nature were included in 
the study. 

Only stability and control characteristics were included in the investigation; no 
consideration was given to the prediction of airplane performance, flutter, or structural 
loads. Both static and dynamic stability characteristics were studied. Stability criteria were 
also examined, but no handling-qualities work was included. 



h. Altitude 
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The investigation was fui er restricted to the consideration and development of 
analytical methods; experimental methods were of interest only to the extent that they 

b 

served as a standard for comparing and evaluating the analytical approaches. Although the 
number of theoretical methods considered was large for the initial part of the investigation, 
it was rapidly narrowed for the major task of comparing and evaluating the several most 
promising approaches. 

The approach taken in the investigation was to develop in order: 

• equations of motion 

• stability criteria 

• stability derivative prediction methods 

• stability characteristics prediction methods 

Three mathematical models of an elastic airplane were considered; rigid (generally 
frozen in the reference shape), equivalent elastic*, and dynamic elastic. 

The rigid model admits no structural deflections from the shape in the reference 
motion, it serves as a base and provides a means for evaluating theoretical methods by 
comparing predictions with data obtained from essentially rigid wind tunnel models. The 
rigid model was sufficiently accurate to describe most aircraft prior to the introduction of 
large (weights of over 100 000 lb (45 359 kg)), swept-wing jet aircraft flying at dynamic 
pressures above 400 psf ( 1 9 152 N/m^). 

The equivalent elastic model assumes that all structural deflections are of a 
quasi-static-elastic nature. Air and inertia loads are considered to be in phase with the 
deflections. No structural dynamic effects are included. This model has also been referred to 
as static-elastic or quasi-static-elastic in the literature. For this model, flight test and 
flexible-model wind tunnel data are used as the basis for validation. The equivalent elastic 
formulation is satisfactory wherever there is a reasonable frequency spread between 
structural and control modes. 

The dynamic-elastic model is the most complex of the three. It takes into account 
dynamic motions of the structure as well as in-phase deflections. Usually this case is handled 
by considering 10 to 80 structural vibration modes. A variation of this approach that lends 
itself to computer calculation is the residual-flexibility formulation, which considers the 
conect phasing of the lower frequency vibrational modes but assumes that higher frequency 
modes are in phase with the loads. This approximation gives to the residual-flexibility 
method advantages for compr^r mechanization that may result in more accurate answers 
than would be obtained v.nh the straightforward inclusion of modes with their 

*The terms quasi-static-elastic, ,/atr*. , v • equivalent elastic are used interchangeably 

in this report. 


6 



mathematically correct phasing. The dynamic-elastic formulation is usually required for 
cases where the structural frequencies are near maneuver and control frequencies 
(approximately within a 2: 1 ratio). 

The equations of motion were developed first for the rigid airplane, treating first the 
general equations of motion, then airplane motion perturbations about the reference flight 
path. Next, the flexible airplane equations of .motion were developed, using the lumped 
parameter concept in describing structural flexibility and structural equilibrium. Perturba- 
tion motions were united with the structural motions for both equivalent elastic and 
dynamic-elastic cases. Both residual-flexibility and completely elastic approaches were taken 
in developing the equations to describe the dynamic-elastic airplane case. Assumptions and 
restrictions introduced into the developments were carefully evaluated as to their effects on 
stability and control applications of the equations. 

The approach taken in establishing stability criteria was essentially one of examining 
rigid airplane criteria to determine whether they do or do not apply to flexible airplanes. 
This part of the study was divided into consideration of static and dynamic stability cases. 
The usual mathematical approaches of characteristic equation rooting and time history 
traces to assess dynamic stability were taken. Some consideration was given to energy decay 
methods for possible future application. 

Stability derivative prediction methods for rigid airplanes can involve the use of a 
variety of aerodynamic theories and methods. The approach taken here was to compare and 
evaluate against experimental data those methods based on lifting line theory; lifting surface 
theory; and handbook compilations of theoretical, empirical, and test data. For the elastic 
airplane, the structural methods considered were finite element theory and equivalent beam 
analysis. Structural methods were compared for the purpose of finding the one most 
suitable for preliminary design application. 

Methods for the more extensive task of predicting the stability characteristics require a 
unification of the previously mentioned areas. Static stability characteristics were 
determined by using the stability derivatives. The dynamic stability calculation methods 
investigated consisted of approximate empirical handbook formulas, roots of the small 
perturbation characteristic equations, and time histories from solution of the large 
perturbation equations of motion. 

The approach taken for this part of the study was to investigate the regions of 
applicability of the available methods and to assess the relative importance of static-elastic 
and dynamic-elastic effects on airplane stability. Areas in the flight envelope of fig. 1 where 
specific calculations and comparisons were made to validate the methods are shown by cross 
hatching. Application of the methods to other areas within the envelope was justified by an 
evaluation of the assumptions and restrictions incorporated into the governing equations. 



The report consists of four separately bound volumes: a summary report and three 
appendixes. These also serve as a handbook which describes and discusses the pertinent 
aeroelastic methods. 

The summary report presents the results and conclusions of the study with discussions 
as required for the reader to gain an understanding of the subject. It should be useful for 
managers, those new to the field or experienced only in related fields, and any others 
desiring an overview of the subject. 

The appendixes contain those results of interest primarily to the specialist, including 
the detailed steps and discussion of the various derivations and developments as well as the 
variations in approaches and applications of the methods. Appendix A treats the 
development of the equations of motion and stability criteria. Appendix B develops and 
evaluates methods for determining longitudinal and lateral-directional stability derivatives. 
Appendix C, using the results of the previous work, discusses and evaluates the pertinent 
factors important to the prediction of airplane stability and response characteristics. 

Matrix notation and methods have been used in developing and presenting many of the 
equations in the summary report and succeeding parts of the study. The reader not familiar 
with the notation is advised to spend the short time necessary to leam the basic ideas and 
symbols of matrix methods. For the specialist, a thorough understanding is essential. 
Reference 3, or one of the many other available books on matrix methods, is recommended. 

It is also assumed that the reader has some background in the field of stability and 
control. If this is not the case, or if review is desired, Etkin’s book (ref. 4) is suggested for 
study. 
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3. SYMBOLS 


This list includes the symbols found in the Summary and appendixes. In different 
technologies some of the symbols have different meanings. For example, e means downwash 
angle to an aerodynamicisi, but strain to a structural engineer. In these cases the several 
del .\itions have be m listed after the symbol. 


General 

41 Aspect ratio nondimensional 

[ A ) Steady aerodynamic influence coefficients matrix, meters^/radian 

[ 6 A ) Unsteady aerodynamic influence coefficients matrix, meter^-seconds/ 

radian 

[Ajl, IA 2 ), IA 3 ], Aerodynamic matrices, newtons, newton-meters 
fA 4 ],[A 5 ) 

a Root of characteristic equation, second'* ; lift curve .slope, radian'* 


b 

Ci / 2 

C 2 

C D 

C D, 

c L 


Speed of sound, meters/second 

Vertical tail elastic to rigid lift ratio, nondimensional 

Acceleration, meters/second^ 

Wingspan, meters 

Cycles to damp to half amplitude, nondimensional 

Cycles to double amplitude, nondimensional 

Drag coefficient, D/qS, nondimensional 

Induced drag coefficient, Dj / qS, nondimensional 

Lift coefficient, L / qS, nondimensional 

Rolling moment cpefficient, M x / qSb, nondimensional 
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C N 

*-n 

C P 

C T 

Cy.Cy 

(cl 

Ko) 

fcj 

lc R i 


C R 


c 


c ref 

D 

D; 

A 

[D] 

d 

{dj} 

E 


Pitching moment coefficient, My /qSc, nondimensional 

Normal pressure force coefficient, N /qS, nondimensional 

Yawing moment coefficient, M z / qSb, nondimensional 

Pressure coefficient, (P-I^Vq^, nondimensional 

Thrust coefficient, T/qS, nondimensional 

Side force coefficient, Fy/qS, nondimensional 

Flexibility matrix with reference point fixed, meters/newton 

Flexibility matrix with reference point fixed and with reference 
point rows and columns removed, meters/newton 

Flexibility matrix with reference point free, meters/newton 

Residual flexibility matrix, meters/newton 

Wing chord, meters 

Root chord, meters 

Mean aerodynamic chord, meters 

c for the 707 and c R for the SST, meters 

Drag, newtons 

Induced drag, newtons 

Transformation matrix from fluid to stability axis system, 
nondimensional 

Elastic displacement, meters 

Column matrix of elastic displacement components at the i^ 1 
element, meters 

Matrix of elastic displacement perturbation, meters 

Total airplane perturbation energy, newton-meters; Young’s modulus, 
newtons/meter^; induced drag efficiency factor, nondimensional; 
energy, newton-meters 
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F 

F 

{F} 

< F A> 

< F i> 

{ f t) 

lF e ) 

f U 

r 

M 

{ f A> 

{*t} 

G 

GW 

G 

8 ij 

t 

8 i 

h 


Internal energy density, newton-meters^/kilogram 

Energy decay parameter, nondimensional 

Force, newtons; surface stress vector, newtons/meter^ 

Total force matrix, newtons 
Aerodynamic force matrix, newtons 

Flexibility matrix relating changes in panel centroid deflections to 
unit loads, meters/newton 

Generalized forces at i“ l element, arbitrary dimensions 
Thrust force matrix, newtons 

Flexibility matrix relating panel slopes to unit loads, radians/newton 

Aerodynamic influence coefficients (subsonic), newtons/radian 

Perturbation force, newtons; perturbation surface stress vector, 
newtons/meter^ 

Perturbation force matrix, newtons 
Aerodynamic perturbation force matrix, newtons 
Thrust perturbation force matrix, newtons 
Shear modulus, newtons/meter^ 

Gross weight, newtons 

Structural influence functions in diadic form with reference point 
free, meters^/newton 

•Aerodynamic influence coefficients (supersonic), newtons/radian 
Acceleration due to gravity, meters/second^ 

Unit base vector, nondimensional 

Altitude, meters; specific enthalpy, newton-meters/kilogram; center- 
of-gravity position, nondimensional 
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"Ill 

h n 

(h n - h) 

lL 


^xx’ ^xy’ ^xz 

lyy. Jy z , 1 ZZ 

m, [i] 


>H 

AAA 

Ej,k 

U, k 

J 

K 


K u 


K 


N 


K; 


K; 


K; 


/ 

K B(W) 

K W(B) 
f K ) 


Maneuver point position, nondimensional 
Neutral point position, nondimensional 
Static margin, nondimensional 

Velocity of panel normal to the streamwise direction, meters/second 
Moments and products of inertia, kilogram-meters^ 


Identity matrix, nondimensional 
Horizontal tail deflection, degrees 
Unit base vectors, nondimensional 


Torsional constant, meters^/radian 

Angular deflection at the exposed horizontal tail due to a unit load 
at the tail, radians/newton 

Structural stiffness coefficient, newtons/meter 

Ratio of aircraft nose lift to aircraft wing lift, nondimensional 

Effective change in vertical tail angle of sideslip due to a unit change 
in rolling acceleration measured at the exposed vertical tail, degrees/ 
radian/second^ 

Effective change in vertical tail angle of sideslip due to a unit change 
in yawing acceleration measured at the exposed vertical tail, degrees/ 
radian/second^ 

Effective change in vertical tail angle of sideslip due to a unit change 
in side acceleration measured at the exposed vertical tail, degrees/ 

meter/second^ 

Effect of lift carryover on the body due to the wing, nondimensional 
Effect of lift carryover on the wing due to the body, nondimensional 
Stiffness matrix with respect to fixed reference point, newtons/meter 


1 1 


‘•ill l i I 
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{Klj 

(K) 

m 

k 

(k],Ck3 

L 

i 

«v 

M 

M 

(Mj 

Cm3 


nij, ir^, 


m 

(m] 

M 

N 

N 

n 


Element stiffness matrix, newtons/meter 

Stiffness matrix with respect to free reference point, newtons/meter 

Generalized stiffness matrix with free reference point, newtons/ 
meter 

Thermal conductivity, newton-meters/second-meter-degrees Celsius; 
elastic constant, newtons/meter 2 ; Strouhal number, nondimensionat 

Corrector matrix for influence coefficients, nondimensional 

Lift, newtons 

Moment arm, meters; characteristic length, meters; pressure difference 
across surface, newtons/meter 2 

Wing c re f/4 to horizontal tail c rc f/4, meters 

Wing c re j-/4 to vertical tail c rc f/4, meters 

Direction cosines, nondimensional 

Mach number, nondimensional; mass of the airplane, kilograms 

Moment, meter-newtons 

inertial matrix, kilograms, kilogram-meters 2 

Generalized mass matrix, kilograms 

Direction cosines, nondimensional 

Perturbation moment, meter-newtons 

Mass matrix, kilograms 

Diagonal mass matrix, kilograms 

Yawing moment, meter-newtons 

Normal force, newtons 

Load factor, nondimensional; number of elastically connected mass 
elements used to represent the airplane, nondimensional 



<er 


n l» n 2’ n 3 


Direction cosines of the normal surface, nondimensional 


M 


p 


P,Q,R 


IP} 

P 


P,q,r 


Qi 


{Q} 


q 


^i 


q 


if) 


Unit vector normal to the surface, nondimensional 
Diagonal matrix of panel unit normal vectors, nondimensional 
Period, seconds 

Components of the angular velocity u> in the body axis system, radians/ 
second 

Total pressure, newtons/meter^ 

Aerodynamic panel pressure forces, newtons 

Static pressure, newtons/me ter roll rate, radians/second 

Perturbation components of angular velocity u>p in the body axis 
system, radians/second 

Generalized force, arbitrary dimensions* 

Matrix of generalized aerodynamic and thrust forces, arbitrary 
dimensions* 

Pitch rate, radians/second; rate of internal heat energy addition, newton- 
meters/second 

Generalized coordinates, arbitrary dimensions* 

Dynamic pressure, newtons/meter^ 

Pitch rate, qc re j/2V C j, nondimensional 

Matrix of generalized coordinates, arbitrary dimensions* 

Matrix of generalized coordinates of elastic free vibration, arbitrary 
dimensions* 

Cantilever eignvectors, nondimensional 


*The units of a generalized fprce times the generalized coordinates must be newton-meters. 



Universal gas con&tant, newton-meters/kilogram-degrees Kelvin; 
magnitude of position vector, meters; region of XY plane not covered 
by the airplane or wake, nondimensional 

Reynolds number, nondimensional 

Position vector at an initial instant of time, meters; body force per 
unit volume, newtons/meter^ 

Reference distance, meters; magnitude of the position vector, meters 


Yaw rate component, rb/2V c ^, nondimensional 

Position vector relative to the body axis system, meters; position 
vector relative to the fluid axis system, meters 

Position vector of the center of gravity relative to the fluid axis 
system, meters 

Position vector relative to the stability axis system, meters 

Position vector relative to inertia! space, meters 

Position vector of the center of gravity relative to the inertial space, 
meters 

Position vector in the undeformed airplane relative to the body axis 
system, meters 

Matrix of airplane position and orientation perturbations, meters, 
radians 

Reference area, meters^; airplane’s projection on the XY plane, 
nondimensional 

Diagonal matrix of panel areas, meters^ 

Complex frequency function, 1 /seconds 

Kinetic energy, newton-meters; thrust, newtons; time, seconds 

Time to damp to Vi amplitude, seconds 


Time to double the amplitude, seconds 



-1/T r 

->/T s 

t 

t* 

U 

u, v,w 

u> v, w 
u i 

A 

U 

{ U M U p} 

V 



i v p} 

<V 


w 

{X} 

X, Y, Z 
x, y, z 


Rolling convergence mode root, 1 /seconds 
> 

Spiral mode root, 1 /seconds 
Time, seconds; airfoil thickness, meters 
Nondimensionalizing time factor, seconds 
Potential energy, newton-meters 

Components of velocity V c in the body axis system, meters/second 

t 

Perturbation components of the velocity in the body axis system, 
meters/second 

Generalized coordinates, nondir. ensional 

t • 

Forward velocity component, u/V c ^, nondimensional 
Generalized elastic displacements, meters 
Lyapunov function, nondimensional; volume, meters' 5 
Equivalent airspeed, meters/second 

Velocity vector of the airplane center of gravity, meters/second 
Velocity vector, meters/second 

Perturbation velocity vector of the airplane center of gravity 
meters/second 

Matrix of airplane linear and rotational rate perturbations, meters/ 
second, radians/second 

Matrix of airplane linear and rotational acceleration perturbations, 
meters/second radians/second^ 

Weight, newtons; airplane's wake projection on the XY plane, 
nondimensional 

Matrix of panel centroid distances to the reference point, meters 

l , l 

Body-fixed-axis system (app. A); fluid axis system (app. B) 

I * . t . , t • , i i . I i * l i I i > > . • 

I I s . i . < * » 
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X B’ Y B’ Z B ; 
X B> yB’ Z B 

X o< Y o> Z o 

x'.y'.z'; 

Y 

rAyJ 

Z R 

{Z} 

I 1 
{} 

1 j 

N 

f 1 T .{ } T 
M - 1 

IK 111 
(0) 

{i} 

(1 

Greek Symbols 
a 

a R 

a ref 

m 


Body-fixed-axis system 

Axis system fixed to a material point 
Earth-fixed-axis system 

Side force, newtons 

Matrix of spanwise panel widths, meters 

Vertical displacement of structural reference point, meters 

Matrix of vertical displacements of each panel from equilibrium, 
meters 

Square matrix 
Column matrix 
Row matrix 
Diagonal matrix 
Transposed matrix 
Matrix inverse 
Determinant of a matrix 
All zero elements 
Column matrix of ones 
“Jump” in enclosed quantity 

Angle of attack, radians 

Angular rotation of structural reference point, radians 

•*» 

Angle between X body axis and V c ^ , radians 
Matrix of panel slopes, radians 
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0 


Angle of sideslip, radians 


r 

r G 

y 

A 

6 


{«} 

{^s} 


e 


€ 


a 


e 

e 

e s 

e ix, G iy, e iz 

0 


(M^ - 1 ), nondiinensional 
Circulation, meters^/second 

Structural influence functions with reference point fixed in diadic 
form, meters^/newton 

Flight path angle, radians: ratio of specific heats for air, 
nondiinensional 

Finite change in some parameter, nondimensionai 

Control surface deflection, radians; arbitrarily small number, non- 
dimensional; Dirac’s function, nondimensionai; thickness ratio, 
nondimensionai 

Matrix of displacements relative to a space-fixed inertial system, 
meters 

Matrix of flexible displacements relative to the structural axis system, 
meters 

Downwash angle, radians; arbitrarily small number, nondimensionai; 
strain, meters/meter 

Change in downwash angle at the stabilizer per unit change in wing 
angle of attack, defeat, radians/radian 

Damping ratio, nondimensionai; nondimensio.ialized coordinate, 
nondimensionai; dummy variable, nondimensionai 

Efficiency factor, nondimensionai; coordinate, nondimensionai; 
dummy variable, nondimensionai 

Euler angle, radians 

Perturbed Euler angle, radians 

Streamwise rotation of panel, radians 

Node rotations, radians 
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Rate of, change of Euler angle, radians/second 



£ » 9 *£ 


Rotational rate of paneled airplane about axis of rotation, radians/ 
second 

Rigid-body rotation about center of gravity, radians 
Angle mode matrix, radians/meter 

Eigenvalue, nondimensionai; taper ratio, nondimensionai; bulk 
modulus, newtons/meter^; Lame’s constant, newtons/meter^; sweep 
angle, degrees 

Roots of characteristic equation, 1 /seconds 

Reduced mass parameter, nondimensionai; Lame’s constant, newtons/ 
mcter^; extent of influence region, nondimensionai 

Cantilever mode shape matrix, nondimensionai 

Matrix of all cantilever modes, nondimensionai 

Poisson’s ratio, nondimensionai 

Coordinates, nondimensionai; dummy variables, nondimensionai 
Constant, 3.14159. . ., nondimensionai 
Density, kilograms/meter^ 

Normal stress, newtons/meter^; density ratio, nondimensionai; real 
root of characteristic equation, I /seconds 

Rotation of structural reference axis system, radians 

Rectilinear translation of structural reference axis system, meters 

Coefficient of viscosity, kilograms/meter-second; shear stress, 
newtons/meter^; time, nondimensionai 

Total. velocity potential, meters^/second; Euler angle, radians 

Normalized natural free vibration modes of the airplane, nondimensionai 
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0 

0 

[ 0 ] 

[0] 

J 

%< 

<P 

0(t) 

* 

>P 

f 

S? 

w 

co n 

"p 

Subscripts 

A 

a 

ac 

b 


Perturbation velocity potential, meters; perturbed Euler angle radians 

i 

Rate of change of Euler angle, radians/second 
Free-vibration mode shape matrix, nondimensional 
Rigid-body mode shape matrix, nondimensional 
Stress diadic, newtons/meter^ 

Normal mode of generalized coordinate, nondimensional 

Velocity potential, nondimensional 

Arbitrary positive function of time, arbitrary dimension 

Euler angle, radians 

Perturbed fiuler angle, radians 

Rate of change of Euler angle, radians/second 

Inertia diadic 

Phase angle, radians 

Frequency, radians/second; imaginary part of a pair of complex roots, 
i /seconds 

Undamped natural frequency, radians/second 
Perturbed angular velocity, radians/second 

Aerodynamic; airplane; aileron 
Aerodynamic 
Aerodynamic center 
Body reference axis system 
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eg 

CP 

D 

E 

E 

Eff 

EqEl 

exp 

F 

HB 
h, ht 
I 
i 

L.E. , LE 

Is 

P 

R 

r 

S 

sp 


Center of gravity 
Center of pressure 
Dutch roll mode 

Equivalent elastic (Formulation II); elevator 
Equivalent elastic (Formulation I) 

Effective 
Equivalent elastic 
Experimental 
Flutter 

Handbook methods 
Horizontal tail 
Inertia relief 
Lower surface 
Leading edge 

Lifting surface theory method 
Phugoid mode 
Rigid; rudder 

Rolling convergence root mode 
Spiral root 
Short period 


s 


Stability axis system; spiral mode 



si 

t 


u 

v, vert, V.T. 

W 

VVB 

WBT 

WT 

0 


oo 


Sea level 
Tip; total 
Upper surface 
Vertical tail 
Wing 

Wing-body 
Wing-body-tail 
Wind tunnel 

At « = 6 jg = ijj = 0°; initial state 

Steady state motion variables; trimmed condition 

Undisturbed condition 
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4. ASSUMPTIONS 


Assumptions used in developing the equations and methods are listed here for 
reference. Where appropriate in the summary report, pertinent assumptions used in 
obtaining a result or equation are given. However, discussions of the assumptions as they 
come into the developments are given in the appendixes. Further descriptions and 
justifications are included in those discussions. 


23 



General Assumptions 



Airplane mass and mass distribution are constant with time 
No thermoelastic effects considered 
No electromagnetic effects considered 
Symmetric airplane 

Variation of air density with altitude is negligible 
No gust effects considered 
Gravitational forces on the field are negligible 
Small perturbation theory 
Large perturbation theory 


Origin of coordinate system is at the center of mass 
Arbitrary perturbations 
Aerodynamic Assumptions 


© 

Potential flow theory 

© 

Thin body 

© 

Slender body 




High aspect ratio 

Prandtl boundary layer approximation 

Perfect gas, thermally nonconducting and chemically nonreacting 

Isentropic # flow 

Steady flow 
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Unsteady flow 
lnviscid flow 
Quasi-steady flow 

Aerodynamic influence coefficients for nonzero sideslip 
Continuum flow 
No finite shock waves 



Velocity field is irrotational 


Structural Assumptions 



Hooke’s law applies 

Only small strain and displacement gradients are considered 



Structural damping is negligible 



Structural perturbations can be represented by normal modes 
Completely elastic math model of elastic airplane 
Residual elastic math model of elastic airplane 
Equivalent elastic math model of elastic airplane 
Rigid math model of elastic airplane 

Airplane displacement vector field is such that the center of gravity 
does not displace or rotate 



X component of elastic deflection is negligible 
Y component of elastic deflection L negligible 

i i I . ! ii i , 

The structure can be adequately represented with beams 

* lit.,. 1 1 : 

Inertia of each finite mass element about its center of gravity is 
negligible i . 
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Dynamic Assumptions 



Free flight only 

© 

No spinning rotors 

© 

Steady-state curvilinear flight 

© 

Steady-state rotation is small 

© 

Zero-lag thrust derivatives 

© 

Ct - is negligible 

L 9 

© 

Cy • ,C Y - ,C| ;• , and C " are negligible 
* Pj * r| Yj Y| 

© x 

Cr, is negligible 

■ q 

© 

Steady-state rectilinear motion 

© 

Stick-fixed-and-unaugmented airplane 

© 

Thrust perturbation forces are negligible 

© 

Steady state, wings level, and zero sideslip 

© 

Level flight' (steady state) 
» 

© 

Linear aerodynamic stability derivatives 

© 

' Two-degree-of-freedom longitudinal motion 
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5. EQUATIONS OF MOTION 


The equations of motion for rigid airplanes and flexible airplanes are presented in this 
section in the forms deemed most suitable for assessing airplane stability characteristics. The 
stability of an airplane is its tendency to persist in a steady reference motion when it has 
been disturbed from that reference motion. The steady reference motion must be defined 
and the perturbation motion analyzed. Thus, four sets of equations of motion are presented. 
Two govern the reference and perturbation motions for rigid airplanes and two govern the 
reference and perturbation motions for flexible airplanes. 

The perturbation equations of motion are presented for three different order-of- 
magnitude approximations of the size of the perturbation motion variables. These are 
termed arbitrary, large, and small perturbation equations of motion. The smail perturbation 
equations are those which have had greatest application in the analysis of the stability of a 
steady reference motion. Arbitrary and large perturbation equations may be used in 
investigations involving large disturbances that upset the airplane and in the study of 
maneuvering flight. However, since neither upset nor maneuvering flight are within the 
scope of this study, the small perturbation equations of motion are given primary emphasis. 

The equations of motion presented for rigid airplanes are those which are familiar to 
engineers working airplane stability and control problems. They are essentially those 
developed by Etkin (ref. 4); however, here they are expressed in a body-fixed-axis system 
that is convenient for the application of aerodynamic influence coefficient methods. This is 
in contrast with the usual formulation in the stability axis system (ref. 4). In addition, the 
formulation provides for very general steady-reference motions including curvilinear flight. 

There is considerable diversity in the manner in which the equations of motion for 
flexible airplanes have been formulated. Modifications have been made in the past to 
accommodate new methods for predicting the aerodynamic forces on the airplane. This is 
apparent from a review of works on this subject such as those by Bisplinghoff and Ashley 
(ref. 1 ), and Milne (ref. 2). 

The formulation used here is called a lumped parameter formulation. It facilitates the 
use of aerodynamic influence coefficients that relate a change in aerodynamic force on a 
small region of the airplane’s surface, a panel, to an average change in flow incidence at 
another panel. Th ; s aerodynamic representation is particularly well suited for making 
empirical corrections to account for separated flow, viscous phenomena, and other 
aerodynamic phenomena that are not readily predicted theoretically but are of considerable 
concern to the stability and control engineer. 



The basis for the lumped parameter formulation is described briefly in this section and 
in more detail in app. A. A complete description may be found in Bisplinghoff and Ashley 
(ref. 1 ). It usually is used in flutter analysis, although it hag also been used to predict the 
longitudinal stability of flexible airplanes. It is the basis for the prediction of the 
longitudinal stability characteristics of the SST configuration appearing in this volume and 
in app. C. The formulation appearing in this section and in app. A is more general. This 
extended generality is included in the expectation tnat an aerodynamic influence coefficient 
method for wing-body combinations in nonsymmetrical motion can be developed. It also 
incorporates the useful approximation called residual flexibility, which is introduced by 
Schwendler and MacNeal (ref. 5). 

The presentation appears in the form of a derivation. This form of presentation is used 
merely to introduce the principles on which the equations are based and to delineate the 
most important approximations which are included in their formulation. The derivation is 
not complete, although it may appear extensive in the flexible airplane case because of its 
inherent complexity. A detailed derivation appears in app. A; related derivations concerning 
aerodynamic and structural theories are provided in app. B. 

Figure 2 summarizes the various forms for the equations of motion appearing in this 
section and gives a brief description of the approximations that characterize them. 


5.1 Equations of Motion For a Rigid Airplane 

5.1.1 General equations of motion.— The equations that govern the motion for a rigid 
airplane follow directly from fundamental principles of mechanics. These are the laws of 
conservation of momentum, which state that the rate of change of linear momentum is 
equal to the total force applied to the airplane and that the rate of change of angular 
momentum about the center of gravity is equal to the total applied force couple (or 
moment) about the center of gravity. They are stated analytically as 


0 

dv c - r- 

M — j~=M g + JF as 

(i) 

(D 

s 

j ^ r ^ 


@ 

(</> # 0>) =y rxFdS 
S 

(2) 
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Formulation 


Character 



Nonlinear; dynamically couple^ 
longitudinal and lateral-directional 
motions; no simplifying approxima- 
tions 


Nonlinear; dynamically coupled 
longitudinal and lateral-directional 
motions; no simplifying approxi- 
mations; no restrictions on mag- 
nitudes of perturbation variables; 
stability analysis by time histories 


Nonlinear; dynamically coupled 
longitudinal and lateral-directional 
motions; no restrictions on mag- 
nitude of perturbation variables, 
except that orientation angle per- 
turbations admit sm a as a, 
cos ass 1, Sin 0^0, and cos0=s 
1, etc; equations of motion for 
reference motion are separable; 
stability analysis by time histories 


Linear; dynamically uncoupled; 
equations of motion for reference 
motion are separable; stability 
analysis by characteristic roots 


FIGURE 2. - MODIFIED FORMS OF EQUATIONS OF MOTION 





where 

M = total airplane mass 

V c = velocity of the airplane c.g. relative to earth-fixed-axis system (x ,y' ,z ), fig. 3 

- » 

g = gravity force per unit mass * 

F = surface force per unit area 

S = airplane’s total surface 

ip - inertia tensor expressed as a diadic 

= i I x i + j Iyj + k I z k - i I xz k - k I zx i 

w = rate of rotation of airplane relative to e^rth-fixed-axis system 
T = position of an arbitrary Q point relative to airplane c.g., fig. 3 


w 



FIGURE 3. - AXIS SYSTEM FOR THE UNDEFORMED SHAPE OF AN ELASTIC AIRPLANE 


The time rate of change d/dt is that apparent to an observer in the earth-fixed-axis system 
x,y',z\ which is assumed to be an inertial reference frame. Thus, letting 3/3t represent the 
time rate of change apparent to an observer in a body-fixed-axis system x,y,z and letting a 
dot (•) represent the partial derivative 3/3t, the equation of motion may be expanded as 
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(3) 


M (V c + u x V c ) = Mg + F a + F t 


ip • w + uix (~ -w) = M a + M t 


(4) 


where 

f A + F T V s fdS 
and 

m a + m t = 

and where and are the total force and moment due to aerodynamics and Fy and Mi- 
are the total force and moment due to the thrust of the propulsion units. 


♦'D 


r x F dS 


Equations (3) and (4) represent the general equations of motion for a rigid airplane. 
The x,z plane of the body axis system has been made to coincide with the plane of 
symmetry of the airplane so that the inertia diadic contains only the I xz product of inertia. 
Certain assumptions have been used in writing these equations. They are: 

Free flight only 

Origin of coordinate system at the center of gravity 
Airplane mass and mass distribution constant with time 
No consideration of electromagnetic effects 
No spinning rotors 
Symmetric airplane 



Equations (3) and (4) may be expanded as six scalar equations, but before they may be (g 4 
so expanded the gravitation force vector must be written in terms of the body axis system. 

A convenient form is obtained in terms of Euler orientation angles. These angles orient the 
airplane in space relative to the earth-fixed-axis system, as shown by fig. 4. Let z'in the 
earth-fixed system be directed positively down toward the earth’s center, and let x' have 
some specified direction. Let axis system xi,yj,zj be initially codirectional with the x'.y'.z' 
system but with origin at the airplane’s c.g. Introduce the rotation $ about the z\ axis, 
which rotates a system X2,y2> z 2 suc h that X 2 is in the plane of airplane symmetry. Then 
rotate a X 3 ,y 3 ,Z 3 system about yj through the angle 6 so that X 3 coincides with x. Finally, 
rotate about X 3 through the angle © so that y 3 coincides with y. This completely orients 
the airplane relative to the earth-fixed-axis system. One may choose alternate definitions of 
the Euler angles; the definition presented here is the one used by Etkin (ref. 4, p. 100). 
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FIGURE 4.- EULER ANGLES 


The gravity force vector is in the direction of z\ i.e./® = gk\ The transformation of gk* 
to body-fixed axis is given by 

gk* = - g sin 9 i + g cos 9 sin# j + g cos 9 cos <p k 

Now, the general equations of motion may be written as six scalar equations in terms of the 
body axis system as 


I 

xx 


P - 


I 


xz 


MU + M (QW - RV) 
MV + M (RU + PW) 
MW+M (PV + QU) 
(R + QP) + (I zz - I yy ) QR 


= - Mg sin 9 + F a + F_, 

\ x 

= Mg cos 9 sin # + + F T 

y y 

= Mg cos 9 cos # + F^ + F t 


= m a +m t 

x X x 


(5a) 

(5b) 

(5c) 

(5d) 
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(5e) 


*yy « + ‘xz < p2 * e2 > + «zz ' V QR ' % + M T y 
■zz » - ! xz < P - «“> + <'xx * : zz» PR = M A + “T 


(5f) 


These are six equations in the eight unknown quantities, U, V, W, P, Q, R, 0 , and 9 . The 
aerodynamic and thrust forces and moments in the right-hand members are functions of 
these eight quantities or their derivatives with respect to time and, possibly, control 
variables. Th<= control variables are regarded as known functions of time. 

A complete c et of equations is obtained by introducing the following two kinematic 
relations obtained from the Euler angle definitions: 


9 = Q cos 0 - R sin 0 (6a) 

0 = P + Q sin 0 tan 6 + R cos 0 tan 9 (6b) 


The heading angle for the airplane may be obtained from the additional kinematic relation 

0 = (Q sin 0 + R cos 0) sec 9 (7) 


Equations (5) and (6) are eight equations in eight unknowns. They may be integrated if the 
aerodynamic and thrust terms are specified functions of the lotion variables, but they are a 
nonlinear system of ordinary differential equations. If initial data, consisting of U, V, W, P, 
Q, R, 9 , and 0 specified at an initial instant of time t 0 are given, then in general it is 
possible to integrate equations (5) and (6). This results in a determination of the variables at 
times later than t Q . Because of tne nonlinearity of these equations, however, integration is 
possible only by mechanical quadratures, except for exceptional cases consisting of steady 
motion of the airplane. 

Fortunately, the cases of steady motion that satisfy equations (5) and (6) represent 
solutions that are of prime interest. They are the subject of the following paragraph. The 
stability of an airplane in these steady motions, i.e., the tendency for the airplane to persist 
in the steady motion when disturbed from it, is evaluated using perturbation forms of the 
general equations of motion, equations (5) and (6). These perturbation forms are discussed 
in pars. 5.1.3 and 5.1.4. Paragraph 5.1.4 also introduces expressions that describe the 
dependence of the aerodynamics on the motion variables. 



5.1.2 Steady reference motion.- In a steady reference motion all of the derivatives 

with respect to time in equations (5) vanish. The equations of motion must be satisfied for 
the airplane to be in equilibrium. Therefore, the steady velpcities (denoted with a subscript 
one as Uj, Vj, Wj, Pt, Qj, Ri) must satisfy 


I 

xz 


M(Q 1 W 1 - RjVj) = - Mg sin 0 + F A + F T 

x i x i 

M(R 1 U 1 - PjWj) = Mg cos d l sin 0.^ + F A + F^, 

y x y 

M ^ p i v i " ^i u i^ = Mgcos 5 i cos ^i + f a + f t 

Z 1 Z 1 

- I xz < ¥V< I zz- I yy ,< ¥ ( l =M A x + M T 

l x l 

(P 1 2 -K 1 V(I xx -I 22 )P 1 S 1 =M a + M t 


I Q.R, + (I -I )P n Q =M a + M-, 
xz H l 1 ' yy xx' 1^1 A T 


(8a) 

(8b) 

1 

(8c) 


(8d) 

(8e) 

(8f) 


and the kinematic relations 

= ^1 C0S ^1 " R l sin 0 l 

+ (Qj^ sin + Rj^ cos 0^) tan 9 ^ 


(9a) 

(9b) 


Rates of change of the Euler angles 9 1 and 0 i have been admitted, but it must be 
recognized that this leads to time rates of change in the components of gravity force in the 
right-hand members of equations (8a), (8b), and (8c). The left-hand members of these 
equations are invariant in time. Thus, in the case of motion in which 9 j and O j are 
nonzero, to maintain steady motion the aerodynamic forces and moments, or those due to 
thrust, must be controlled in such a manner as to balance the gravity force changes. In a!! 
cases for which 9 j and 0 ( are zero the aerodynamic and thrust terms are constant during 
the steady motion. 
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Four of the most important steady reference flight conditions for assessing the stability 
and control characteristics of large airplanes are level flight, climbing flight, turning flight, 
and pullup. The steady velocity components and Euler angles and the equations of motion 
for those four cases are as follows: * 

a. Steady, level, rectilinear flight: 

Vj = Pj = = 0^ = 0 , 0^ = constant 


Mg sin 

6 1 + F A 

+ 

T 

= 0 


X 1 

x i 



f a 

+ f t 

= 0 


y l 

y i 


Mg cos 

V f a 

+ f t 

= o 


\ 

Z 1 



m a 

+ m t 

= 0 


x i 

X 1 



s 

> 

+ M 

= 0 


Z 1 

Z 1 



t Steady, climbing flight: ... 

The equations are identical to equations (10), although it must be noted that 
atmospheric density variation will lead to unsteady motion for constant control 
settings. However, this unsteady motion may be regarded as negligible for shallow 
climbs over moderate time periods. 

c. Steady, turning flight: 

= 0 , 0^ = constant and 0^ = constant 

MQjWj = - Mg sin 0 1 + F A + 

X 1 X 1 

M - PjWj) = Mg cos (h sin <t>, + F A + Frp 

y l y l 

- MQ^U^ = Mg cos 0^ cos ^ + F a + F,p (11) 


- ! xz Ql P l + ^zz - V <Vl - M A + M T 

X 1 X 1 

I X Z < P 1 2 - R 1 2 ) + ' I XX - 1 Z Z > P 1 R 1= M A +M T 


‘xZ + Oyy - l XX> P !<h - M A + M T 


where the kinematic relations, equations (9), may be used to eliminate Qj and Pj by 
introducing 
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Qj =Rj tan0 x 

tan 0, 


p = _ R i 

1 1 COS 0 . 


d. Steady pull-up: 


= Pj = = 0 and 0^ = 0 

= - Mg sin ^ + F A + F T 

X 1 X 1 

f a + f = 0 

y l y l 

- MQjUj^ = Mg cos 0 X + F a + F t =0 

z i h 


a 


m a 

+ m t 

= 0 

X 1 

*x 


m a 

y l 

+ M T 
y l 

= 0 

m a 

+ M T 

= 0 

z i 

Z 1 



The kinematic relation, equation (9a), gives 


VS 

To maintain a steady pull-up it is necessary that 


Mg cos 0 X Q x = F a + F t 
X 1 X 1 


and 


Mg sin 0 X Q x = F a + F T 
Z 1 Z 1 
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5.1.3 Perturbation equations of motion. - Many of the major objectives of a stability 
and control evaluation of an airplane are achieved by considering motion deviating from a 
specified steady reference motion of the airplane. Th& deviation is often termed a 
perturbation, and the variables that describe the motion are taken to have values consisting 
of the sum of the reference value and a perturbation value. The notation of par. 5.1.'.' is 
used to identify the reference values of the variables, while a lowercase letter or a subscript 
p is used to identify the perturbation values, e.g., 


U = U x + u 


The perturbation equations of motion are obtained from the general equations of 
motion, equations (5) and (6), by replacing the velocities, accelerations, and Euler angles 
appearing in those equations by the sum of a perturbation and a reference value. The 
reference values correspond to steady motion so that their derivatives with respect to time 
vanish. 


The substitution of perturbation and reference values for the values of the variables 
appearing explicitly in equations (5) and (6) does not constitute an approximation. 
However, an assumption is made regarding the aerodynamic and thrust terms. These terms 
are functions and are, in general, functions of the velocity of the airplane and its 
acceleration. It is assumed that the functionality is such that it is possible to write 



(13) 


Clearly, this eliminates functionality that is transcendental as in the case of the gravity 
forces that have transcendental (harmonic) dependence on the Euler angles. 

The perturbation form for the equations of motion may now be written as 


Mu + - RjVj + QjW + W^q + qw - R^v - V^r - rvj 


= - Mg sin (0 X + 0 ) + F a + F t 


+ f* +fn 


(14a) 


x. 
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(14b) 


Mv + M(R^Uj - P^W^ + R^u + U^r + ur - P^w - W^p - wp) 

h + V sin (0 1 + 0 p 


= Mg cos{Q, + 0J sin (0, + 0„) + F a + F t *+ f A + f T 


y x y x y y 

Mw + M(P 1 V 1 - + PjV + V^p + pv - Q^u - U^q - qu) 

= Mg cos (0 + e ) cos(0 1 + 0 ) + F a + F t + f A + f™ 
y y Z 1 z z 

‘xx P - ha (f + QjP + + OT) + <1„ - V (Q l r + R l q + rq) 


= m A +m 

X X 


I q + I (2P,p + p - 2R,r - r ) + (I -I ) (P-r + R..p + rp) 
yy xz ' r 1 ' ' xx zz' ' 1 1^ w 


m A +m T 

y * y 


l zz * " X xz < p ~ Q l r " R l q - qr) + (I yy " W (P l q + Q 1 P + pq) 


m A +m T 
z z 


(14c) 


(14d) 


(14e) 


(14f) 


Equations (14) are called the arbitrary perturbation equations of motion. The tenn 
“arbitrary” is used because no restrictions have been placed on the admissible magnitudes of 
the peiturbations. They govern the motion of a rigid airplane to tl same order of 
approximation as the general equations of motion, equations (5), except that the 
aerodynamic and thrust terms have been separated into reference and perturbation terms. 
That separation led to a simplification of equations (14d), (14e), and (14f). Since the values 
of the reference motion variables satisfy the equations of motion themselves, those terms 
involving only reference motion vanish. This simplification is not possible in equations 
(14a), (14b), and (14c) because of the sine and cosine functional dependence of the gravity 
force components. 

Large perturbation equations of motion are obtained by restricting perturbations of 
the Euler angles 0 and 0 to a magnitude less than about 7.5 deg. Then the trigonometric 
functions of these angles are approximated by 
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cos 0 « 1 
p 


sin 0 « 0_ 
P P 


(15) (G9 


COS 0 *s 1 

P 


sto V»p 


Introducing the approximations of equations (15) into equations (14a), (14b), and 
(14c) leads to 


Mu + M(Q 1 w + W x q + qw - R^ - - rv) 

= Mg 0 cos 0 + f A + ip 
p L x x 


(16a) 


Mv + M(RjU + U^r + ur - P^w - W^p - wp) 

= Mg (-0 sin 0. sin 0 + 0 cos 0 cos 0 ) + f. + f 
P 1 1 P 1 1 L v, 


(16b) 


Mw + M(P^v + VjP + pv - Qu - U^q - qu) 

= - Mg (0 p sin cos 0j^ + 0 p sin 0 X cos 0 ^ + f A + f T 


(16c) 


Equations (16) along with equations (14d), (14e), and (140 represent large 
perturbation equations of motion. It should be noted that terms involving products of 0p 
and 0p have been deleted. This is consistent with the approximations represented by 
equations (15), which are truncated power series expansions of sine and cosine functions. 
The system of equations is linear in terms of the perturbation Euler angles but nonlinear in 
the perturbation velocities. 

Small perturbation equations of motion are obtained by deleting products of the 
perturbation motion variables. Equations (14) reduce to 

Mu + M(Q.w + W q - R n v - V.,r) = Mg 0 cos 0, + f * + f™, . (17a) 

a 1 1 1 P 1 A v 
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Mv + M jr - P x w - Wjp) = Mg (- 0 p sin ^ sin 0 X 

L • 7, cos 0.) + f A + tr- 


im 



M\v + M(PjV + VjP - QjU - Uj q = - Mg ( sin x cos 

+ 0 sin 0. cos 0j) + f A + f T 
P z z 


( 17 c) 


I^P ' l xz + Q 1 P + P l q) + (I zz " I yy ) ( Q i r + R i q) = m A + m T (17d) 

X X 

>yy 4 + 2 'xz (P 1 P - B l r) + (1 xx - ‘zz> < P l r + R l« * m A + m T < 17e > 

y y 

] zz ■ ‘xz (p ■ Q l r ■ E q) + (I yy ' W ' p i* + QyP) = m A + m T < 17f) 

z z 


and the kinematic relations, equations (9), become 



0 p = -Qj 0 p sin 0j + q cos 0 X - R x 0 p cos 0 X - r sir 0 X 
0 p = P + (Q 2 0 p cos 0 X + r sin 0 X - Rj 0 p sin 0 X 


+ q cos 0 X ) tan 6 ^ + (Q^ sin 0^ + Rj cos 0^) 0 


P 


( 18 a) 


( 18 b) 


The terms in the small perturbation equations of motion appear to be of consistent 
order of magnitude. They are a set of linear ordinary differential equations provided the 
aerodynamic and thrust terms are written as linear functions of the perturbation motion 
variables. The form of the aerodynamic and thrust terms is the subject of par. 5.1.4 and app. 
B. 
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;An additional simplification is possible if consideration is given to the oroers of 
magnitude of the reference motion variables. In the analysis of airplane response, the 
perturbation motion variables have small but finite values, the limits on the validity of the 
application of the small perturbation equations are set primarily by the limit of validity of 
linear aerodynamic theory. In many cases, this limitation may "pply equally to the 
prediction of the aerodynamic forces associated with the reference motion; hence, in these 
cases the reference motion rotation rates and sideslip velocity must be limited to the sairr 
orders of magnitude as those of the perturbation variables 

Equations (17^ apply to the case of stability and control analysis when the 
aerodynamic terms are obtained from a linearization of the aerodynamics about some 
condition where the aerodynamic phenomena are basically nonlinear. The perturbation 
motion variables in this case are then limited in magnitude by the range of validity of the 
linearization. Aerodynamic terms in this case are obtained empirically from test data. 

Many important stability and control evaluations are performed on cuses where the 
magnitudes of the reference motion variables are sufficiently small to permit their products 
with the perturbation variables to be neglected. This is the case where stability and control 
characteristics of transport aircraft in their cruise condition are concerned. For these casts 
the small perturbation equations may be reduced to 

Mu + MW 1 q=Mg 0 cos 0 x + f A + f T (19a) 


Mv - MW.p--Mg (0 sin 0., sin 0 , - 0 cos 0 , cos 0.) + f. + f rp 

1 P * i I P 1 i A 1, 


(19b) 


y y 


M\v - MtLq = -Mg (0 sin 0 , cos 0 , +0, sin 0 1 cos 6 .)+ f A + fip (19c) (D4) 

l pi 1 ip i l a L 


z z 


I p - I r - m . 
xx xz A 


+ m n 


x 


(19d) 




+ m r 


(I9e) 
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and the kinematic relations reduce to 


0 = q cos 0, - r sin 0. 

p 1 1 

0 = p + (r sin 0 j + q cos 0 j) tan 0 ^^ 


( 20 a) 

(20b) 


5.1 .4 Aerodynamic derivatives. — The concept of an aerodynamic derivative or, more 
commonly, a stability derivative is developed in app. B. It relates a component of 
aerodynamic force or moment acting on an airplane to a motion variable or a parameter that 
describes a change in airplane shape such as a control surface deflection angle. The change in 
-airplane lift with change in airplane angle of attack is expressed as 

r = J- d -k 

L a ~qS 8a 

The notation of a partial derivative is appropriate, since all other motion variables and 
shape parameters are held constant while the angle of attack is allowed ; uy. Besides, the 
lift is not, in general, a linear function of angle of attack or any of the other variables that 
influence .he airplane lift. The stability derivative Cl , therefore, varies with the airplane’s 
flight condition. For perturbation motion of the airplane about a reference motion, the 
stability derivative is evaluated for the reference motion and may be considered to be 
constant in the perturbation motion. This local linearization, discussed more fully in Sec. 7 
of this volume and in detail in app. B, is a widely accepted practice (refs. 4 and 6). The basic 
mathematical theory is given by reference 7. 

Another widely accepted practice is to introduce a stability axis system (ref. 4, p. 103). 
This is a body-fixed-axis system having X s forward and in the direction of the freestream 
velocity in the reference motion of the airplane. The body axis system X,Y,Z is rotated 
about the Y axif from the stability axis system by an angle a re f, as shown in the following 
sketch: 




where Vj is the freestream velocity vector in the reference motion so that Vj = -Vcj. The 

purpose of the stability axis system is to introduce a simplification in the equations of 
motion, which is discussed in ref. 4 but is not used in this study. 


The reference angle of attack is used but is denoted as otj, i.e., the angle of attack in 
the reference flight condition. Lift and drag of the airplane are measured in a wind axis 
system such that in reference flight 


and 



cos sin a ^ 



= - Dj sin 0 ;^ 


- Lj cos a j 


In the perturbed motion 

F a = - D cos (o^ + ar) + L sin (a + a) 
z 

F. = - D sin (a + a) - L cos (a +a) 

7 * A 


Using these relations and the small angle approximation such that 

sin (a^ + a) sin + a cos 


and 


cos (a^ + a) cos a ^ 


a sin a 


1 
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the perturbation aerodynamic forces 


f . = - cos a* D + sin a. D.a + cos a, L a + sin a. L 

A x X p x x xi x p 

+ sin a, P a + cos ‘--La 
1 p 1 p 


( 21 ) 


f . = - cos a. B, ft-sina. D - cos a. L + sin a. L, a 

A 11 lp lp 11 

z 

- cos a.D a + sin a. L a 
1 p 1 P 


for small perturbations, products of perturbation quantities are neglected so that 


f . = - cos a, D + sin a, P a + cos a 1 L, or + sin a. L 

A x Ip lp 11 lp 

f A = - cos a D, a - sin a 1 D - cos a L + sin a, L, a 
A 11 lp lp 11 

z 


( 22 ) 


The perturbation components of aerodynamic force and moment are written in terms 
of stability derivatives as (app. B) 


D = q, S 
p 1 w 


K +2 ° D i)t; + V* l+ % 


0b 


2 V 


+ on + C 


D 


ct 


_2L£_ + c 
D- 2 V D 
a c- p 


pb 


2 V 


+ C 


qc 


D 2 V 
q c 1 


(23a) 


+ C. 


D 


rb 


2 V 


+ ? C D a 6 i 

i 6. 
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L = q, S 
p H 1 w 


K + 2 %)^ + %h + %, 


0b 


2 V 


+ C T a+ C 


cue 


+ C, 


La L d 2 V Lp 


pb 


2 V 


+ C 


qc 


Lq 2 V 


+ c. 


rb 


2 V 


i 6 ; 


( 23 b) 


f A = 


+ c yS p-,Z C 4l l 

J r • c^ x •'d. I 


- P b 
'y 2 V 
P 1 


( 23 c) 


m . =q, S b 

A.. 1 w 


2C i T" + % 0 + C Iq 2 V + C Ji 2 V 
1 p 0 p 


rb 


+ C 1 2 V + ^ “I. 

r C 1 1 4 1 J 


? C (. 4 i 


(^ 3 d) 


m . =q, S c 

A n w 


/c + 2 C \-*r~ + 

V m u ra l/ V Cl 


C5 l^l +C n 1 - 

m 0 m /3 


/ 3 b 


2 V 


+ C ct + c + c P b 

V mi 2 V l m 2 V 
“ U C 1 P 


+ C 


qc 


m 2 V 

q c. 


+ c: 


m 


rb 


2 V 


£ C 6.1 

‘ “ a i 1 I 


( 23 e) 


m . = q, S b 

A_ n l w 


2 \t;\ 

n 2 ^ + C „ + ? C „. «il 

p Ci " «! 1 «i J 


( 23 f) 
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The perturbation components of thrust forces and moments are 





b Kt + ?v) 

( 2CT m 1 ^ + ?\’’ i ) 

( 2 \V?V j ) 


= < li s 

^1 w 


= q, S c 
^1 w 


= q, S b 
1 w 


(23g) 


The above results may be combined with the perturbation equations of motion of par. 
5.1.3. The perturbation force components in the X and Z directions must be computed 
differently for large and small perturbations. Equations (21) must be used to relate 
perturbation lift and drag to the X and Z directions for large perturbations; for small 
perturbations equations (22) must be used. 


5.2 Equations of Motion for a Flexible Airplane 

5.2.1 Lumped parameterization. — In the preceding section, the motion of a rigid 
airplane was described in terms of six degrees of freedom of the airplane’s center of 
gravity-three translational and three rotational. The effects of inertia end gravity can be 
completely accounted for in terms of forces and moments applied at the airplane’s center of 
gravity. And, although the distribution of aerodynamic surface stresses depends on the 
shape of a rigid airplane, its dependence upon the rigid airplane’s motion may be expressed 
entirely in terms of the motion of its center of gravity. Finally, the effects of the 
aerodynamic surface stresses on the rigid airplane’s motion are accounted for by 
aerodynamic forces and moments acting at the airplane’s center of gravity. 
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When an airplane is flexible, however, motion which is relative to the center of gravity 
must be considered. Equations of motion written entirely in terms of center of gravity 
motion variables and in terms of forces and moments considered to act at the center of 
gravity are no longer sufficient. The additional n; u.on is called elastic motion, and it is 
common terminology to refer to “elastic” degrees of freedom. 

The flexible airplane is essentially a continuous, elastic body. As such it has a 
continuous infinity of elastic degrees of freedom. Thus, unless some simplifying approxima- 
tions are introduced the equations of motion for the elastic motion will be integrodifferen- 
tial equations (app. A, equation (6.35)). To avoid this complication the airplane is divided 
into a large number of elements. Each element is regarded to have at most six degrees of 
freedom. If their number is n, the equations of motion of a flexible airplane will govern 
motion in 6n degrees of freedom and be 6n in number. Certain approximations will be 
introduced to reduce the number of equations, but the number of degrees of freedom will 
remain effectively unchanged. The motion in some degrees of freedom will be dependent 
motion. This is the subject of several later paragraphs in this volume as well as of much of 
app. A. 

This paragraph, as well as several of those following, makes use of matrices in the 
formulation. When matrices of a particular type or matrix operations are introduced, or 
properties of a matrix are discussed, an article will be cited from Frazer, Duncan, ai d Collar 
(ref. 3 ,. The reader may refer to the cited article for the appropriate explanation. 

Letting the total volume of the airplane be represented by V and the volume of the i*h 
element by Vj, then 


n 

V = £ V. 
i~l 


(24) 


The total mass of the airplane is M and <he density is p^. The mass of the i^ element is, 
therefore, 



P A dv 


(25) 


Also, recalling that r is the position relative to the airplane’s center of mass, the position of 
the i*h lumped mass is taken to be at the point 




r P A d V 


(26) 
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This is the centroid of mj and, although refinements may be introduced, the density is taken 
to be uniform over the element of volume sc that 

(27) 

which is the location of the geometric centroid. 

The geometry of the flexible airplane is defined as shown by fig. 5. 



2 ! 

FJGU. AXIS SYSTEM FOR THE DEFORMED SHAPE AND THE UNDEFORMED SHAPE 

OF AN ELASTIC AIRPLANE 


The elastic displacement of a mass particle is a vector, d, which is a function of 
position and time. It describes the change in position relative to the airplane’s center of 
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gravity. This elastic deformation carries each mass particle from its position, 7, in the 
undeformed airplane to its position, T, in the deformed airplane. Thus, the elastic 
displacement vector is a function of position and time as * 

d = d (r , t) 


where it may be noted that the dependence 0.7 position has been written as a function of its 
undeformed position,?. This is referred to as a Lagrangian description (ref. 8, p. 29). 

For the lumped parameter formulation each volume element is taken to displace as a 
whole with the mean displacement given by 

*rzrJv. <k* dv ~ < 28 > 

1 1 


The totality of all elastic displacements is written as a column matrix (ref. 3, Art. 1 .2) as 





(29) 


Elastic rotations relative to the airplane as a whole occur for a mass particle and are 
given as 

jr v x d 


When the airplane structure is regarded as a continuous body, the rotation of the mass 
partitles does not contribute to the inertia forces. When the airplane is represented by a 
collection of elements with lumped masses, it is tempting to include rotational degrees of 
freedom and the rotational inertias of the lumped masses. But on consideration of the 
practical numerical problem of predicting the motion of the lumped masses, it is better to 
neglect rotation and include only the inertia forces generated by their mean translational 
motions. 

The number of dynamically participating degrees of freedom that may be included in a 
numerical analysis is limited. Let the upper limit of that number be m. If rotational degrees 
of freedom are included dynamically, the total number of lumped masses is m/6. If they are 
ignored, the number of lumped masses may be twice as great, i.e., m/3. The larger number 
of lumped m isses more closely represents the continuous airplane. The smaller number must 
have mean rotations that are not readily defined as well as moments and products of inertia 



that are even more difficult to define rationally. Therefore, the formulation that ignores the 
dynamics of mean rotations of lumped masses is the better approximation. 

t 

It must be clearly understood that rotational motion of the lumped masses has not 
been set to zero in the above. In fact, the aerodynamic surface stresses are dependent upon 
the surface slope of the airplane, and the surface slope changes as the airplane deforms 
elastically. The lumped surface slope (or mean slope) will be related to a mean rotation of 
the element, so that these quantities enter the problem through the aerodynamics even 
though their dynamic influence is neglected. 

It is desirable to have a convenient method of computing the mass properties of the 
entire airplane from the values of the lumped masses and their positions relative to the 
center of gravity of the airplane, rf. To do this consider the manner in which the moments 
and products of inertia are introduced in mechanics, i.e., as a consequence of the integral 
expression for angular momentum about the center of gravity: 

Ip 'U> = / v rx M /^dV 


This expression may be replaced by 


w * £ Jy r x (wx*r) p A dV 


n ^ 

= £ r. x (wx r.) m 
• i=l 1 11 


so that the entire angular momentum is replaced by the sum of the contributions to the 
angular momentum by each lumped mass. Similarly, the total linear momentum is 


»VE m V 
C i=l 1 ' 


(30) 


Introduce the diagonal mass matrix (ref. 3, Art. 1.2), [m^ , as defined in app. A, pa>\ 
6.3.1, i.e., 
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► 



and the rigid-body mode shape matrix [0] defined as 


[ 0 X ] 

10 2 ) 



(32a) (G 10 


where 



1 

0 

0 

0 

z. -y. 

1 

0 

1 

0 

- z. 

0 x. 




1 

1 

0 

0 

1 


-x. 0 
X 


(32b) 


Also, introduce the velocity matrix, written in transpose form (ref. 3, Art. 1 .2) defined as 


{V } T = |_U, V, W, P, Q, RJ 


(S3) 



so that components of momentum may be written as 


* 

MV *T 

c 

MV,'! 

MV • k 
c 

(0*<u)* i 

(?*£>)'! 

_ (0 • w) • k 

where the definition of the inertia matrix [M] is given by 

[M]=[0] T M[0] (35) 

Thus, the components of linear and angular momentum for the airplane are obtained in 
simple matrix expressions, equations (34) and (35). 

The aerodynamic and thrust forces act at the external surface of the airplane. Let the 
total external surface, S, be subdivided into the external (exposed) surfaces of the lumped 
masses, Sj, so that 


= 101 T M(01 {V}=tM]{V} (34) 


n 

s = E s. 

i=l 1 


(36) 


Then the total aerodynamic and thrust forces on the airplane may, in accordance with 
equation (3) and this subdivision of the surface, be represented as 
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~ where Fj is the <oi >v on the lumped mass. Similarly, from equation (4) 


m a + m t 



r x F dS 



Now, introduce the matrix definitions of the components of the aerodynamic forces and 
moments: 






(37) 


the components of the thrust forces and moments: 



Fq, } F„ , Irp > Mrp ) Mq, , IVlpp 

L T x y T z T x y T zJ 


(38) 


and the components of element surface forces: 



(39) 


Using the ngid-body mode shape matrix [0] , it follows that 

KHM = [*f { p } <40) 


so that a simple matrix expression relates th surface stresses to overall airplane forces a id 
moments. 



The results presented in this section form the basis^for representing a continuous, 
flexible airplane as a collection of lumped masses. However, neither of two important details 
have as yet been considered. These are the dependence of elastic constraining forces, and the 
dependence of the aerodynamic forces, on the mean elastic displacements {d}. The elastic 
constraining forces are discussed in par. 5.2.2. The aerodynamics are discussed in Sec. 7. 

The thrust forces may seem to have been slighted in the discussion. However, their 
direct effect in terms of forces applied at the engine attachment points may be readily 
introduced in an obvious manner and the induced effects on the aerodynamic forces must 
be accounted for in the stability derivatives. 

5.2.2 Structural flexibility. — The flexibility of the airplane is fundamental to the 
considerations of this report. This paragraph introduces concepts to describe what is meant 
by flexibility in physical terms and presents the equations that describe it. An exposition of 
the theory underlying these equations is the subject of par. 4.2.3 of app. B. It suffices to 
note here that the underlying theory is that for infinitesimal deformation of a perfectly 
elastic, isotropic, and homogeneous solid. It is often called the classical theory of elasticity 
and is a linear approximation of an exact, nonlinear theory of elasticity. 

The change in shape of an airplane when it is subjected to external loads is a 
manifestation of structural flexibility. The magnitude of the change, viz; deflection or 
deformation, is in approximately constant proportion to the load producing it. Thus, the 
deflection measured at a point A on the structure is linearly related to the load P at another 
point B. This is expressed as 

(S) 6 = CP (41) 

where the constant of proportionality, C, is termed the flexibility of the structure associated 
with the two points, A and B. 

A general description of deformation follows by considering two positions of each 
point of the structure. They are the position ~ before deformation and the positionTafter 
deformation, both of which are measured relative to the airplane’s center of gravity. The 
deformation or displacement vector is given by 

d=T-? (42) 

The displacement vector was originally introduced in par. 5.2.1, where the mean 
displacement of an element of the airplane was defined by equation (28) as 
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d * “ ”t /v. 


P A d dv 
A 


In the following, the mean displacement vector at the i^ 1 element is relat. . 
applied at the i^ 1 element as well as all other elements. 

The mean surface force on the jth element is: 


F. = 

J 



Fd S 


An element of the airplane also experiences body forces such as inertia and 
These are proportional to the airplane density, and the mean values for the 
defined as 


R. 

] 



7— 

where d z r /dt~ is the acceleration relative to the earth-fixed reference system. 

The mean displacement of the i^ 1 element due to forces applied to the 
defined as 


d i 




and the displacement of the i** 1 element due to forces at all elements is given h\ 


d. 

i 


n 

£ 

j = i 


O.. 

1J 



1 V 




Further, this expression is valid for all elements, so that there are n equation-' 
(45b). 



The subscript zero on the diadic C Q . is used i j indicate that the structure is clamped at 
the center of gravity, i.e., the point on the airplane that was at the c.g. before deformation is 
constrained against translation or rotation. That point may be allowed to translate as d Q 
relative to a reference axis system with origin at the airplane’s c.g. It may also be allowed to 
rotate through an angle 0 O relative to the reference axis system. Thus, the mean 
displacement dj relative to the reference axis system is given by 

d. -d -r. xfl ! C„ -Ir.- F.\ (46) 

1 1 1 " i = l “ij | 1 I 


By defining the column matrix 



the matrix definitions of par. 5.2.1 may be used to write equation (46) as 

{ d H*]W = -[ c o]({ R H p }) <48) 

where the flexibility matrix (ref 3, Art. 8.1 1), [C Q ] , is made up of n submatrices 3 by 3 in 
size and given by 


c 

C 

O 

oxx.. 

& 

O 
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ij 
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so that the flexibility matrix [C Q ] is 3n by 3n in size. The typical element c °zyy is thc 
component of deflection in the z direction at the i^ 1 element due to a unit component of 
force in the y direction at the j^ 1 element. 



Equation (48) is the result required in the equations of motion of a flexible airplane 
represented by lumped parameters. However, the flexibility matrix must be computed from 
consideration of the structural details of the airplane. The flexibility matrix so derived will 
not be in the form defined for equation (48). The structural analysis treats the airplane 
structure as a continuous body. Thus, the flexibility matrix [C 0 ] must be obtained from a 
continuous analogue using the type of averaging process that led to definitions of mean 
displacements and forces in par. 5.2.1. 

The problem of computing the flexibility matrix [C Q ] is complicated by the fact that 
the structural analysis leads to a flexibility matrix that relates displacements at points to 
forces applied at points. The method of computation of [C G ] is not immediately obvious. 

The structural analysis deals with an idealized structure made up of a collection of 
simple structural elements (app. B). For example, a typical wing structure is made up of 
panels of skin cut out by ribs and spars, the webs of the ribs and spars, as well as the rib and 
spar caps. The structural properties of each of these elements are known. The structural 
analysis proceeds by requiring the forces on the elements to be in equilibrium and by 
requiring the elements to remain joined together under load, i.e., continuity. The problem 
arises from the fact that this leads to a flexibility matrix in terms of forces and 
displacements at the points where the structural elements are joined. 

Joining points appears inaccurate for the joining of some elements, e.g., a spar cap and 
a spar web. These elements obviously join along a line. The structural analysis will not 
satisfy continuity along the entire line of intersection. It will satisfy continuity only at a 
finite number of points, called node points, on the intersection. 

One method of structural analysis, the displacement method, is readily illustrated by a 
simple one-dimensional example. For greater detail, reference is made to app. B, Consider an 
assemblage of springs: 


K K, E 

a b c 

12 3 4 


The springs have stiffnesses represented by K a , Kjj, and K c , and the nodes are numbered as 
shown. There will be four nodal forces, Fj and four nodal displacements, dj. Even though 
forces exist in the springs at all points between nodes and all interior points undergo 
displacement, only the nodal values will be required for the structural analysis. 



The element stiffness matrix for the spring from node 1 to node 2 is a 2 by 2 matrix. 
The first column is the forces at the nodes for a unit displacement of node 1 to the right, 
Le., Fj = K a , F 2 - -K a . The second column is the forces at the nodes for a unit displacement 
of node 2 to the right, i.e., Fj = -K a , F 2 = K a , so that the element stiffness is given by 



Similarly, two other element stiffness matrices may be written as 



The stiffness matrix for the entire assemblage is obtained by forming a composite matrix as 
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and may be used to write 

\F\ = [k] jd| (49) 



The stiffness matrix [K] is singular (ref 3, Art. 1.8). Thus, the displacements are not 
uniquely related to the forces. This has to be the case, since an equal addition to each of the 
displacements does not cause a change in the forces. That addition would represent a 
translation as a rigid body. 
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The structure may be clamped at node 4 such that d 4 = 0. Now, the displacements are 
determinant. Let the matrices be partitioned (ref. 3, Art. 1 .7), as 

% 


Pel' 


" K ll 

k h 


W 

_ F 4_ 


_ K 41 

K 44_ 


0 


so that 


Pel * [ K ll] l d el 

F 4 * l K 41j i d el 


The reduced stiffness matrix |Ki jj has an inverse (ref 3, Art. 1.11), i.e., 

M 1 


and it follows that 


Pel * [Kll]' 1 i F ,l < 5l > 

Finding the composite stiffness matrix to be singular in the above example is typical. 
The composite stiffness matrix for an entire airplane is also singular. The singularity is 
removed by setting enough nodal displacements to zero that the airplane structure is 
constrained against translation and rotation. The requirement is satisfied if a single point is 
clamped against rotation and translation. When this has been done the flexibility matrix 
may be obtained by inverting the reduced stiffness matrix. The displacements predicted by 
this matrix are changes in position relative to a coordinate system with origin fixed at the 
clamped point. 
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The example illustrates the use of element stiffnesses to ^escribe the elastic properties 
of a structure in ter: :s of forces and displacements at the nodes. The forces applied to the 
airplane are distributed. Also, the displacements of interest are mean displacements. Thus, 
neither the force nor the displacements used in the aeroelastic analysis are those used in the 
structural analysis. This “interface” problem must be solved when the finite element 
structural analysis methods are used to define the elastic properties of the structure. 

When the aerodynamic surfaces of the airplane are of reasonably high aspect ratio, say 
6 or larger, the interface problem discussed above may be resolved by treating the 
aerodynamic surfaces as beams. The actual structure is idealized as a beam lying along the 
elastic axis (or locus of centers of flexure, ref. 8) of the actual structure. Torsional stiffness 
of the elastic axis is defined as the GJ distribution. The flexural stiffnesses of the elastic axis 
are defined as the EI X and Ely distributions. The volume elements, Vj, of the lumped 
parameter description are attached to the elastic axis by rigid, massless members. The 
structure appears as in the following sketch of a wing. 



Elastic axis 


The lumped masses are located at the centroids and are connected to the elastic axis by the 
rigid, massless members shown as dashed lines. 

The stiffness matrix may be developed for this idealized structure by the finite- 
element, displacement method just described. The flexibility matrix obtained will be 
appropriate for the lumped parameter formulation, [C Q ] . No interface problems occur and 
slopes associated with the lumped masses are readily computed. 


5.2.3 Internal equilibrium. — The basis for writing internal equilibrium equations for 
the airplane was introduced in par. 5.2.2. There it was noted that the st.-ictural deflections 
for the lumped parameter formulation are represented by ' 

0 HI-®] ! b I ■- [ c 0 ] (l R l- l p l) < 52 > 

The matrix {R} is the matrix of body force components acting on the lumped masses. 

The body force vector on the i™ lumped mass is due to acceleration and gravity. It is 
expressed as 



P A dV 


Carrying out the integration results in 




m. 


(53) 


where d^r|7dt^ is the mean acceleration of the centroid of the i^ 1 lumped mass relative to 
the earth-fixed-axis system, x',y\z\ Letting the angular velocity of the airplane about its 
center of gravity be represented by the vector w, the acceleration of the i™ lumped mass 
may be expanded to give 


d 2 r/ av 9*d. 9<L 

J _ c , -r„ sr . x , „rr.. i 


dr 


2-r 


at +WXV C + ^2“ + 2wx at 


+ %r ■ xr, + u> x (wxr.) 

Oil 1 


(54) 
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where 


i 

i 


is the time rate of change apparent to an 
observer rotating with the airplane 



is the velocity of the airplane center of 
gravity relative to the earth-fixed-axis 
system, x' ,y* ,z' . 


It is assumed that the rotation rate of laige, flexible airplanes is sufficiently small that 
the Coriolis forces and a portion of the centrifugal forces may be neglected by comparison 
with the other inertial forces, i.e., 


9d, 


2<nx — gj- m.»0 and c Jx (cox r.) = 0 


(55) (D4 


so that the inertial force on the i™ lumped mass is approximately 




, 2 — » 

d r. 


dt" 


K - 3u ^1 d % 

= m i [ar +wxV c + "aT xr iJ +m i^ 


(56) 


In the case of steady reference motion of the airplane, as in par. 5.1.2, the inertial 
forces acting on the mass elements reduce to 


m. 



(57) 


where, again, the subscript 1 denotes the values of the motion variables in the steady 
reference motion. Also, as in par. 5.1.2, the gravity force on the lumped masses is 

gj = m. g (-sin © 1 i + cos 0.^ sin 0^ i + cos 0^^ cos k) (58) 


The perturbation equations of motion were obtained in the case of the rigid airplane in 
par. 5.1.3. It is shown in app. A, par. 6.3, that the perturbation equations of motion for the 
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( 53 ) 


rigid airplane may be expressed in the lumped parameters as 



where {Vp} = {v} - { Vj } is the matrix of perturbation velocities corresponding to the 
velocity matrix of equation (33). The perturbation matrix { r^} represents perturbation 
of the position and orientation of the airplane’s center of gravity. 



V V *pJ 


(60) 


The matrices fMj] and [M 2 I are defined in app. A by equations (6.119) for small 
perturbations and equations (6.120) for large perturbations. The matrix (Mp contains 
coefficients arising from centrifugal accelerations and f M 2 1 contains coefficients appro- 
priate to gravity forces. 

The perturbation body forces acting on the lumped masses of the airplane are given by 

H=Hjiy + W[{v p }*[» 1 ]|v p }*[M 2 ]{r;j| ,61) 

so that the perturbation form for internal equilibrium, equation (52), may be written as 


Hp(-[©] ■ [ C J M {i'Vl *[<»]{ R (} + [C 0 ] |F| ( 62 a) 


where 


{"!= ! Vp h [ M i] | V P S * t M J l r 5pl < 62b > 


This result and equation (59) constitute the perturbation equations of motion for a flexible 
airplane. 
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Internal equilibrium lor the airplane in reference motion may be expressed as 


i d l}-$3 { B i} = -[ C 0 ]{M “ ( F i>} (63) 


where 


[0 1 { a^} are the lumped mass accelerations in the 
reference motion and 

l 0 1 { S] } are the reference gravity vector components. 

A more conve. ient formulation of internal equilibrium may be achieved by 
introducing new variables for the elastic displacement {d} . Let 

{d} = [0]{u} (64) (s?) 

where the matrix! 0) is a transformation matrix (ref. 3, Art: 9.3) and the matrix { u} is a 
column matrix of generalized coordinates. 

It is desirable to require that the generalized coordinates be linearly independent. This 
is not the case with the elastic displacements {d} , which are displacements relative to the 
airplane’s center of gravity. As such, from app. A equation (6. 1 30), it follows that 


[0] T M{d} = O (65) @ 

There are 3n components of '{d} and, since equation (65) represents six linear relations, 
there are 3n-6 linearly independent generalized coordinates {u} . The transformation 
matrix [ 0 ] has 3n rows and 3n-6 columns. 

A particularly convenient choice for the transformation [ 0 ] consists of letting the 
columns of this matrix be the eigenvectors generated from an eigenvalue problem 
constructed from the internal equilibrium equations. 

Let the airplane be at rest, or at most in uniform translational motion in free space, 
and let the airplane be in a state of vibration. Internal equilibrium, equation (52), is then 
represented by 


{<*}-[0] = - [C 0 ] M{d} (66) 
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which, on utilizing equation (65) to determine the elements of { B} , may be written as 


{d} = - [C] [m] {d} . (67a) 

{ 

where 

[c] = [[i] - [ 0 ] [M] [0] [C Q ] (67b) 

Introducing the transformation, equation (64), into equation (67a) results in 

[0]{u} = - [c] [m] {0}u (68) 

where it has been assumed that the transformation matrix is not a function of time. Each 
column of the transformation matrix may be taken to be a solution to 

{0} 4 u. = - (cl frnj {0.} ii. 


which, by a separation of variables (ref. 9, p. 430), may be written as 


[[I] - u. 2 [c] [m]] {0.} = 0 


and 


2 

u. + w. u. = 0 
i ii 


(69a) 


(69b) 


where w p the separation constant. 

Equation (69a) represents the eigenvalue problem that was being sought. It has 3n-6 
linearly independent solutions, eigenvectors, { 0 j} corresponding to 3n-5 eigenvalues up. 
The eigenvectors represent free vibration mode shapes for the airplane. From equation (69b) 
it is readily apparent that u> j is the frequency of a simple harmonic motion. This is a natural 
frequency of the freely vibrating airplane. The eigenvectors {pj} are the mode shapes, i.e., 
the shapes into which the structure deforms when vibrating freely. 
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The transformation of coordinates, equation (64), may now be given more meaning. 
The displacements {d} that describe the elastically deformed shape of the airplane are a 
linear combination • 


3n-6 

{d}= £ 

i=l 1 1 


of shapes { 0 j} . The generalized coordinates uj are the amplitudes of the shapes. By 
properly adjusting the amplitudes, the elastic deformation of the airplane under the action 
of any system of self-equilibrating loads may be represented by equation (64). Thus, the 
generalized coordinates Uj are completely satisfactory for expressing elastic motion. A 
considerable advantage accrues through their use. The internal equilibrium equations, when 
transformed, assume a form much simpler than that of equation (62). 

The perturbation internal equilibrium equations may be written in terms of stiffness as 
[K] {d p } - - {dp} + {f} - [m] [0] {R} (70a) 

where 


IKl “[[cj' 1 - [c 0 ]' l [«ll { (®f [C^" 1 I®]} 1 l©] T [Cj- 1 ] (70b) 

This result is derived in app. A as equation (6.147). On introducing the transformation and 
multiplying equation (70) by the transpose of the modal matrix, it is found that 

[ic]{u} +[m] { i}= [0) T {f} (71) 

where use has been made of the orthogonality properties 

[0] T Ik] (0) = [k] (72) 

and 

(0] T [nr.] [01 =[m] (73) 
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from app. A, equations (6.157) and (6.158). In addition, it may be noted that the term 
corresponding to rigid-body accelerations and gravity { R| vanishes. This is a consequence 
of equation (65), which may be used to write • 

[0] T M[0]=O (74) 

This last result may be interpreted as an orthogonality of the free vibration mode shapes 
{0 j} , the columns of [ 0 ] , with the columns of [ 0 ] when weighted by the mass matrix. 
The columns of [ 0 ] may be regarded as rigid-body mode shapes, and that term will be 
applied in the following in calling [ p ] the rigid-body mode shape matrix. 

Two obvious advantages have accrued from the transformation. The generalized 
stiffness and mass matrices, t'Kj and fmj , are both diagonal matrices. Thus, equation 
(71) is uncoupled in the left-hand member. If the right-hand member is replaced by a 
column matrix of generalized perturbation forces, i.e., 

|Op}=[0] T lf| <™> 

internal equilibrium in terms of the generalized coordinates expands into 3n-6 scalar 
equations as 


a i a i + *i u i = Q Pl 

ft 2 a 2 + *2 “2 " \ 

♦ 

m 3n-6 U 3n-6 + Sn-6 U 3n-6 = Q P 3n _ 6 


(76) 


When internal equilibrium is expressed by writing out the equations as they are here, the 
advantage gained from the use of the generalized coordinates uj is readily evident. Only one 
of them occurs in each of the left-hand members. This is the meaning attached to the term 
“uncoupled.” The only coupling in equations (76) is in the right-hand members. The 
generalized perturbation forces Qpj are coupled. In general, each is a function of uj, uj, uj 
and the motion variables u, v, w, p, q, r and their time rates of change. That dependence will 
be discussed later in the subsection regarding aerodynamic derivatives. 



It is important to note that the motions governed by the internal equilibrium equations 
are not motions relative to an inertial reference frame. The elastic displacements {dp} in 
equation (70) and the generalized displacements {u} in equation (76) describe motion 
relative to a body axis system. The body axis system may be rotating and accelerating so it 
is not, in general, an inertial reference frame. Equations (70) and (76) are not exact 
expressions of Newton’s law of mass times acceleration equated to applied force. They are 
approximations, unless the center of gravity of the airplane is moving with a constant 
translational velocity. The approximation is valid only if the Coriolis and centrifugal 
perturbation forces are negligible by comparison with the other forces acting on the airplane 
as noted by equations (55). 

5,2.4 General equations of motion. — In e preceding subsection the lumped 
parameter formulation was used to formulate internal equilibrium for perturbation motion, 
equations (62a) and (76). Internal equilibrium for the airplane in its reference motion was 
given by equation (63). Also, the perturbation equations of motion for motion of the 
airplane as a whole were obtained by analogy with the rigid airplane perturbation equations 
of motion. They were given by equation (59). 

The system of equations listed is complete for the analysis of the stability of an elastic 
airplane with two exceptions. The dependence of the aerodynamic and thrust forces on the 
airplane’s motion has not been delineated, and equations of motion for the motion of the 
airplane as a whole have not been given for the reference motion. Those equations are 
obtained from general equations of motion. That derivation is precisely the same as the 
derivation for the rigid airplane. Equations (8) and (9), for reference motion of a rigid 
airplane, followed directly from the general equations of motion for a rigid airplane, 
equations (5) and (6). 

The only distinguisiiing feature in the application of these results is obscured at this 
point of the development. The aerodynamic and thrust forces and moments depend on the 
airplane’s flexibility. If the moments and products of inertia are denoted by the subscript 
one to indicate that they are evaluated for the airplane shape under the steady reference 
loads, the reference motion equations of motion may be written immediately as 


M (Q x W x - R x V x ) = -Mg sin ©^ + F A + F T 

X^ K 

M (R 1 - P. W^) = Mg cos ©^ sin 0 . + 

*1 

M (P^ U^) = Mg cos ©^ cos 0^ + F^ 

Z 1 


i 

(77a) 

+ F T 

*1 

(77b) 

+ f t 

Z 1 

(77c) 
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to 


+m t 

■* XX. X«^ 

(77d) 

W < p iW> + «xx,-V P A =M a +m t 

1 y l 

(77e) 

: xz, Q 1 B 1 + V V P 1 Q 1 “ M A * M T 

X XX Zj^ 

(77f) 

©1 = cos 0^- R^ sin0^ 

(78a) 

+ (Q^ sin 0^ + R 1 cos 0^) tan 



The methods for finding the reference shape of the airplane cannot be described 
without a detailed description of the aerodynamic and thrust terms. That description is the 
subject of par. 5.2.5. The problem of computing the reference airplane shape is formulated 
in par. 5.2.6. 

5,2.5 Aerodynamic derivatives. — Up to this point the aerodynamic forces have been 
introduced simply as a consequence of surface stress, and for the lumped parameter 
formulation this led to a column matrix of components of mean perturbation aerodynamic 
forces {f} . The perturbation aerodynamic forces depend on the perturbation motion of 
the airplane. It is stated in app. A, equation (6.172), and shown in app. B that the 
dependence may be approximated by the expression: 


<f> - [Aj {V,} ♦ [Aj {V p } 4 [A S J (dp} 4 [AJ {dp} + [A 5 ](dp} (79) 


As noted in app. B, the coefficient matrices may be determined by wind tunnel 
measurements. However, theoretical determinations of the matrices in equation (79) are 
based on inviscid, small perturbation flow theory. Without entering into the detailed 
derivations of the aerodynamic matrices (Ail ... [A 5] , it is easily shown that equation 
(79) is an appropriate form. 



Solutions to the inviscid flow problem are in the form of a perturbation velocity 
potential that satisfies an equation of flow and appropriate boundary conditions. The form 

b 

of the flow equation has no direct impact on the general form of equation (79). The 
boundary conditions consist of the requirements that the flow be tangent to the surface of 
the airplane and that certain vortex laws be satisfied in the airplane’s wake. As shown in 
app. B, the boundary conditions at the surface contain only the perturbation potential (not 
its time rate of change), the perturbations to the airplane’s translational and rotational 
velocity components, and the elastic displacements and their rates of change. The wu/e 
boundary condition depends on the frequency of the motion ce , but for small frequencies 
this dependence may be eliminated by an approximation. Thus, the perturbation velocity 
potential 0 must have parametric dependence as 

0 = 0(x, y, z 5 t; {V p }, {d p }, {d p }) (80) 

The perturbation pressure at the surface of the airplane is computed from the 
perturbation velocity potential using Bernoulli’s equation. That computation involves taking 
the first derivative of 0 with respect to time. Thus, the pressure force at the i^ 1 element of 
the airplane will have the parametric dependence 

p i ■ p i «y ■ <V' t a p>’ < a pi’ < a p» < n > 

Finally, the aerodynamic theory is a linear one and the pressure is normal to the surface. 
The linear dependence in equation (79) is therefore justified and the components of the 
surface force are proportional to the surface pressure force and the components of a unit 
vector normal to the surface in the reference airplane shape. 

The aerodynamic derivatives for the airplane are computed from equation (79) by 
using the mode shape matrices. The elastic displacements may be replaced by generalized 
displacements by using the transformation equation, (64). Thus, the forces and moments 
acting at the airplane’s center of gravity are given by 

{Vp} + [A 2 ] {v p } + [A 3 ] [0]{u} +[A 4 ] [0] {u} (82) 

♦DyKHa}] 
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The generalized perturbation aerodynamic forces, Qp., appearing in equation (76) are given 
by 


{Q p } =[0] T [[A 1 ] {Vp} + [a 2 ] {VpMAj] [0] (u} + [a 4 ] [0] (u) (83) 

* [A 5 ] [0] {S}] 


An aerodynamic derivative is the rate of change of a component of force or moment 
on the airplane due to a change in a motion variable while all other motion variables are held 
unchanged. Thus, for example, 



is the aerodynamic derivative, relating pitching moment to rate of change of the second 
generalized elastic coordinate. It is found from equation (82) by setting U2 equal to unity 
and all other perturbation motion variables to zero. For an airplane represented as n lumped 
masses, equations (82) and (83) contain (27n2-18n) aerodynamic derivatives. If there are 
100 lumped masses, then there may be 268 200 aerodynamic derivatives. 

Those aerodynamic derivatives which relate forces and moments at the airplane’s 
center of gravity to motion of the airplane as a whole are contained in the first two terms on 
the right of equation (82). These aerodynamic derivatives correspond to those introduced in 
the discussion of the rigid airplane and contained in equations (23). This correspondence is 
delineated in Sec. 7 of this volume and is developed in detail in app. B. 



5.2.6 Airplane reference shape pig shape determination). - The selection of the design 
shape for an airplane is the result of an optimization process involving a large number of 
design parameters. When an airplane has a significant anlount of structural flexibility, 
allowance must be introduced for the difference in the shape of the airplane as it is when 
manufactured in fabrication jigs and as it is when subjected to in-flight loads. In the case of 
a transport airplane the optimization process usually minimizes the drag of the airplane in 
steady, level, midcruise flight subject to constraints such as minimum body diameters and 
wing thickness. The result is a design shape that is the optimum shape of the airplane in a 
design point flight condition. The design shape may differ significantly from the jig shape or 
the shape of the airplane in off-design point flight conditions. Computation methods for 
determining the jig and off-design point shape ire the subject of this section. 

The basic relations required for computation of the off-design point and jig shapes of 
the airplane have already been developed. Recalling equation (63) and making use of 
equation (65), it follows that the elastic displacements in the reference flight condition are 
given by 


{d x } = KG ({a x } + {g x }) - {FjJ (84) 


The subscript 1 has been added to the flexibility matrix and the mass matrix to denote that 
these matrices are dependent upon the reference flight condition mass distribution. 

The values of the elements appearing in the acceleration matrix |aj} and the 
gravitational matrix |gjj. of equation (84) are specified by the reference flight condition. 
The aerodynamic and thrust forces {Fj} must be such that the equations of motion for 
the reference flight condition, equations (75), are satisfied. They may be written in matrix 
form in an analogy with equation (84) as 


[0] T En^] [0] ({ ai } + {g x }) =[0] T {F 1 } 

where 

{ F A l l + { F T 1 ) = 03' I >il 


(85) 


(86) 


as noted by equation (40). 



The reference motion of the airplane is a steady motion. In a steady motion the 
aerodynamic equation, equation (79), reduces to 

k 

{f } = [Aj { V p } + [A 3 ] {dp} (87) 

If it is assumed that the aerodynamic forces are linearly related to the motion variables in 
the reference motion, then equation (85) may be applied to the reference motion. This is 
not always appropriate. Thus, let 

<>y - [Aj {V t } * [A 3 ] {dj + [AJ {6j} + (F^ <88) 

where a bar has been used on the matrices [ Aj] and [ A 2 I to denote that these matrices 
may require empirical determination or, if determined from potential flow theory, may have 
forms that differ from + hose of their perturbation counterparts. Aerodynamic control 
surface deflections have been added as the elements of {6 j} using an influence coefficient 
matrix [A c ] to relate them to aerodynamic forces on the lumped masses. Finally, the thrust 
forces on the lumped masses have been included as {Fj} j with the subscript 1 exterior to 
the brackets to distinguish this matrix from the similar matrix appearing in equation ( 86 ). 

Equations (84), (85), and ( 88 ) may be combined to give 

[M] ({a 1 } + {g 1 })=[0] T |[A 1 ]{V 1 }i-[A c ]{6 1 } +{F T } -[A 3 ][[l] - 
tcj [A 3 ]j ({ ai }4{ gl }) -[Aj] {Vj-tAjlSj ■' T }J} 

(89) 

where the definition of the inertia matrix [M] has been used from equation (35). The 
control deflections { 6 \ } are the only unspecified quantities appearing in equation (89), 
with the possible exception of thrust. These quantities must be adjusted so that the six 
equations represented by equation (89) are satisfied. This balances the airplane into the 
reference flight condition. 
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Having balanced the airplane, the deformed shape may be determined. Combining 
equations (88) and (84) to find ' 


~ "[[I] " [C3] [^3]] [C-|_]|[ m j][0] ({ a 2} + {»i}) " [Aj { v ^} 

- { F t }J (90) 

all quantities on the right are known from the preceding, so { dj} is determinable. 

The airplane in the jig is not subjected to loads, including those due to gravity. Thus, 
vanishes for the jig shape. Letting the quantities on the right be evaluated for the 
design flight condition, the displacements { dj } describe the difference in the design and jig 
shapes. Off-design point reference motions lead to airplane shapes that are much more 
readily obtained from the jig shape. The displacements computed from equation (90), when 
it is evaluated for the off-design conditions, are simply added to the jig shape. 

5.2.7 Perturbation equations of motion. — The perturbation equations of motion were 
presented in a matrix formulation in par. 5.2.3. They are given by 


M ({V p } + [M x ] { V p } + [M^r ' }) =[0] T {f} (59) 

P 


and 

*K]{u} + [m]{u} = [0] T {f} (71) 


v'here the aerodynamic contribution to the right-hand members is given by equations (64) 
and (79) as 


{f} = [Aj] {V p } + [A 2 ] {V p } +[A 3 ][0]{u} + [A 4 ](0]{u} + [A 5 ][0] {11} 


(91) 



Equations (59) are six in number and equations (71) are 3n-6 in number when n lumped 
masses are used to represen the airplane. The stability of the airplane is assessed by 
determining the manner in which the variables { Vp} depend on time. Thus, equations (59) 
and (71) must be integrated to determine { Vp} . That integration requires the kinematic 
relations given by equations (20), i.e.. 


©p = g cos ©j - r sin </>^ (20a) 

0p = p + (r sin 0^ + q cos 0j) tan ©j (20b) 


The above system of equations constitutes the equations of motion for a completely 
elastic airplane. Because the number of equations and unknowns is very iarge, it is of the 
utmost importance, if at all possible, to reduce the number of equations and unknowns. 

One method of simplifying the problem is to use the residual flexibility method. The 
frequency of the motion in many of the elastic degrees of freedom is very large. When these 
frequencies are an order of magnitude or more larger than the frequencies of the rigid-body 
motion, the tv/o motions are only weakly coupled dynamically. However, the elastic 
deflections due to those elastic degrees of freedom may strongly affect the rigid-body 
motion. The effect is brought about by quasi-static deflection of the airplane structure. The 
high frequency of the motion implies that the deflection is very nearly in phase with and in 
constant proportion to the loads producing it. Thus, generalized inertial forces such as rhjiij 
are very small. 

Residual flexibility is formulated by partitioning the generalized elastic deflections 
{u} into two parts, as 



The generalized inertial forces associated with { U2 } are set to zero. Equation (71) is 
partitioned as 


[ 1 3{ u 1 } + Cm 1 J{u 1 } = [C> 1 ] T {f} (92) 
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and 




( 93 ) 


This fonnulation is discussed in considerable detail in app. A. There it is shown that 
equations (92) and (93) may be combined with equations (59) and (91) so as to eliminate 
the generalized coordinates { U 2 J completely. The result is the set of equations governing 
the motion of the airplane’s center of gravity as 


[MlAy+tMjl {V p }+(M 2 ! |r o ' 

p 

f 

i$] T [m - ia 3 i rey]- 1 ^rA x i {v p } 

+ IA 2 ] {V p }+ CA 3 ] r0 1 I {u 2 } - (A 3 1 [ Cr ] fmj [ 0 t l {u x } (6.208) 
-[Ag] tC RmJ f0)({V p }+ fM x ] {V p }+ fM 2 ]|r o ' 

+ fA 4 i[0 1 l{u 1 } + (AgHOj] {ujj 
and dynamic-elastic equations 


ffij4 {aj + nCjl {Uji = !<P 1 l T [t 1 1 - IA 3 1 tc R l] 1 [lAJ { V p } 
+ lA 2 l{V p }+lA 3 l [0J { Uj } -(A 3 ](C B lNiJt0 1 l{u 1 } 
-[Agile K'mj(0]^{v p } +[Mj){v p }+ (M 2 1 { r p }) 
+ia 4 1(0 1 i{u 1 }+ tA 5 )(0 1 ){« 1 }] 


where the matrix ( Cr ] is the residual flexibility matrix. The residual flexibility is given by 


(C R 1 =(cl " 10J [Kj]" 1 [QJ 7 (6.203) 

It may be noted that TK]J * is the flexibility associated with the dynamically included 
generalized coordinates {uj} . Thus, {Cr} is the excess of flexibility of {CJ over the 
flexibility included in the dynamically included portion of the structural behavior. The term 
“residual flexibility” is appropriate. 



The formulation presented by equations (6.208) and (6.209) has not been widely used. 
An alternate formulation used in some of the evaluations presented in this report has had 
prior use. The residual flexibility matrix is also given by 

(m r g 2 rV2] T - [ 5 r]) (-Kn<y < 94 > 

-CfjCKjWfNidpi+iti) 

The alternate formulation follows by substituting equation (94) into equations (6.208) and 
(6.209). This particular formulation has the disadvantage that all free vibration mode shapes 
must be computed with high accuracy. That is not always practical. The result is a less 
accurate evaluation of the residual flexibility than is available with the form given by 
equation (6.203). 

If we introduce generalized coordinates qj with the first six of these set into 
correspondence with the elements of the matrix { r bp } , which appears in the matrix 
formulation of the equations of motion, then to the degree of approximation appropriate to 
the small perturbation equations of motion 


<V T = L 4 r V VV V 4 eJ 

« 2 ' v v v \i 


M t = LV "s q jJ 


and 


Also, letting 


the equations of motion given by equations (6.208) and (6.209) of app. A may be merged 
into a single matrix equation. This equation represents both rigid-body and internal motion 
for a residual-flexibility formulation and is given by 


[al{q} + [b] {q}+ (c] {q } = {0 } 


(95) 
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where 


(a] a 


{b] = 


([M] - [0] r [Bj -1 (A 2 lj [0] T {B^ -1 (A 5 1 (0^ 

[(- (0/f^r 1 FA 2 l) ^m x 4- i0 1 ] T (B 1 ]- 1 tA 5 ] [0^ 

([M J - (01 T [ B x )~ l [A 1 lj [0] T [B^ -1 [A 4 1 [0^ 

(0 I i 1 [Bj] _i t Ajl^ lA J 


[c] = 


^[M] [M 2 ]J 


and 


[0] T [B^ -1 [A 4 1 [0^ 

(fol) - (0 1 ] T [B^ -1 [A 3 ] [0^ 

[B x l sfru - [A 3 1 ^[c] - [0 X 1 tKjT 1 


By letting vanish this set of equations reduces to the equations of motion for an 

equivalent elastic airplane. If all of the elastic degrees of freedom are included as 
dynamically participating, i.e., completely elastic airplane, then [Bj ] = [I] and these 
equations continue to hold. It should be noted the thrust perturbation terms have been 
deleted from equations (6.208) and (6.209) of app. A. 


When the structural dynamics are neglected entirely, only equations of motion 
governing airplane center-of-gravity motion remain. These are called the equations of 
motion of an equivalent elastic airplane. The resulting expressions are found by replacing 
{Crj by [C] in equation (6.208) and setting terms involving uj to zero. The result is 


“«y + [“i] <y + [m 2 ]{ v }> =[0 T [[I] -OsKo]]' 1 [[Ajliy 

+ t A 2 ] <Vp} -O3] [c] [m] + [Mj](v p } + [mJ {r^})] 


( 96 ) 



Some aspects of these various formulations require emphasis. The perturbation 
structural deflections for the equivalent elastic airplane gre precisely in phase with the 
externally applied loads. This is not the case in the completely elastic formulation or in the 
residual flexibility formulation. The maximum deflection, in general, does not occur at the 
instant of time when the external loads are a maximum. Also, the deflections may exceed 
the equivalent static-elastic deflections. This is of considerable importance if the structural 
motion is strongly coupled with the center-of-gravity motion of the airplane. When it is not, 
the problem becomes the concern of the flutter engineer. 

It is of interest to note that the terms on the right of equation (6.208) consisting of 
|01 T [[II - [A,] Ic^]' 1 [lA t l {V p } + IA 2 l{V p }] 


contain the stability derivatives that relate change in the forces and moments at the 
airplane’s center of gravity to changes in the airplane’s motion as a whole. They are 
equivalent to the stability derivatives appearing in equations (23). If the residual-flexibility 
matrix [Cr] is set to zero, these terms contain the rigid airplane stability derivatives. If 
[Cr] is set equal to [C] , the total airplane flexibility, so that equations (96) result, they are 
the stability derivatives for an equivalent elastic airplane. In that case the terms in equation 
(96) consisting of 

-($] T [m-tA 3 l fey]' 1 (A 3 ) fc R l[*m4[0]^{V p }+[M 1 ]{V p }j 

may be included into the equivalent elastic stability derivatives. This is sometimes referred 
to as having included inertial relief The terminology stems from the practice of constructing 
equivalent elastic stability derivatives using static-elastic representations of the airplane. 



6. STABILITY CRITERIA 


Static and dynamic stability criteria are presented aitd discussed in this section. The 
problem of establishing handling-qualities criteria is beyond the scope of this investigation. 
However, physical interpretations of stability criteria and connections with previously 
established handling qualities are pointed out. 


6.1 Static Stability Criteria 

Static stability is defined as the tendency of an airplane to develop forces or moments 
that directly oppose an instantaneous disturbance of a motion variable from a steady-state 
(i.e., equilibrium or trim state) flight condition. For example, when the nose of an airplane 
is raised relative to the flight path and, as a consequence, the airplane develops a nose-down 
moment, the airplane is said to be statically stable for such a disturbance. 

A static stability criterion is defined as a rule by which steady-state flight conditions 
may be categorized as stable, neutrally stable, or unstable. 

In other contexts, the words “static stability criterion” have been used as a 
requirement for an arbitrary minimum static margin. For example, the military specification 
for flying qualities (ref. 10, par. 3.3. 1.1) requires a negative value of C m at all times, 
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which implies a positive static margin. 

In still another interpretation, the Civil Airworthiness Requirements (ref. 11, art. 
4b.l51-155) associate stability criteria with stick-force versus speed behavior. 

The reasons for presenting and defining static stability criteria as given are as follows. 

1. The definitions of stable, unstable, and neutrally stable given are clear so that 
judgment and opinions are eliminated as factors. 

2. These definitions lead directly to important aerodynamic derivatives and show how 
they are related to static stability behavior. 

The reader will notice that these definitions are largely independent of notions of 
stability and stability criteria associated with control force or control surface displacement. 
Specifically, this report does not deal with: 

1. Stick- free stability; 

2. Stability as affected by the feel system including bob-weights; 

3. Stability augmentation systems in general. 

It is recognized that when control surfaces are allowed to float, or when springs or 
other devices are added, the longitudina. stability derivatives and associated control 



characteristics can be significantly affected. Such effects are not discussed directly in this 
-report. However, discussions of the effects of the derivatives do apply to those cases. 

kr 

The criteria for static stability are summarized in table 1. These criteria are equally 
valid for rigid and elastic airplanes. It should be noted that they are presented as expressions 
of local slope behavior. For that reason they apply (as a iocal linearization) to situations 
where aerodynamic forces behave in a nonlinear manner. This is important because airplanes 
do behave in a nonlinear fashion in many instances, for example in stall and pitch-up. 

It should also be noted that, although criteria of table 1 evolve from the definition of 
static stability criteria used here, they vary considerably in importance. For example, 
3 My / 3a(~C mtt ) is of much greater practical importance than 3Fy/ 9v (~Cy^ ). This will 
be discussed in more detail in Sec. 7. Notice also from table 1 that under the adopted 
definition of static stability, the partials, 3M x /3u (~C mu > and 3My/9v (~Cjj^ ), do not 
belong. This implies that for static stability under the current definition, the signs of C niu 
and Cj ip are not important. However, in the practical case these derivatives are important. 
An unusual feature of table 1 is that it includes moment derivatives with respect to 
rotational velocities. Such derivatives are normally associated with dynamic stability and not 
with static stability. The reason for their appearance must be found in the definition of 
statu, stability. The physical justification for including these moment derivatives in static 
stability considerations is that steady-state flight can actually involve constant rotational 
velocities.. 

The physical meanings of the criteria are stated below. 


Criterion 


An airplane is statically stable for a forward speed disturbance, u, if: 3F x /9u<0 

(97) 


The physical meaning of this criterion is that as a consequence of an increase in 
forward speed, u (along the x-axisl, a force must be generated that tends to oppose the 
increase in speed. 


Criterion 

An airplane is statically stable for a side speed disturbance, v, if: 9Fy/3v<0 

( 98 ) 




TABLE L - SUMMAR Y OF STA TIC STABILITY CRITERIA 


General form of 
static stability 
criterion 

Approximate or 
alternate form 

Importance to handling \ 
qualities 

Reference 10 
MIL-F-87 85 

3Fx „ 
3u < 

C x < 0 or 
s u 

Cq >0 
U u 

(No thrust effect) 

Needed for stable phugoid. 
Not important if throttle 
response is good. 

Par. 3.3.6 limits the 
phugoid divergence. No 
direct requirement. 

3F V 

3v ' < 0 

s <0 

(No thrust e 

ffect) 

Helps pilot in perceiving 
the sideslip. Allows 
skidding turns at low 
altitude (wings level). 

Pars. 3.4.3 and 3.4.8 
interpreted to mean 

v° 

9F Z n 

C L« >0 


Primary means for flight 
path control. Significant 
to short period. Always 
satisifed before stall. 

Par. 3.3.3* specifies 
short period requirements. 
No direct requirement. 

9M Z 

w~ >0 

% >0 


Needed to maintain 
straight flight path. 

Pars. 3.4.3, 3.4.4, and 
3.4.5 interpreted to mean 

C n >0 
"0 

3M y 

3a <o 

C m <0 
m Q 


Affects time history of 
pitch response. 

Par. 3.3.1 interpreted 
to mean C m < 0 

3M„ 

<0 

9p 

C c <0 
p 


Affects time history of 
roll response. Affects 
Dutch roll damping. 

Par. 3.4.1* specifies Dutch 
roll requirement. Par. 3.4.16 
specifies roll performance. 

3M V 

— - <0 
aq 

C m <0 

m q 

- 

Affects damping of 
short period (increases 
pitch stiffness). 

Par. 3.3.5* specifies short 
period requirements. No 
direct requirements. 

3M 7 

- <0 

dr 

Cn f <0 


Affects Dutch roll 
damping (increases 
yaw stiffness). 

Par. 3.4.1* specifies Dutch 
roll requirements. No 
direct requirements. 

3M y 

>0 

3u 

c m >o 

u 

i 

Improves speed control. 
Provides warning of 
inadvertent over or 
under speed. Affects 
stick-force behavior. 

No direct requirement, 
but par. 3.3.3 implies that 
violation is allowed 
transonically. 

3M,. 

- - <0 
9v 

V° 


Warns pilot of existence 
of sideslip. Allows 
emergency roll control. 
Affects Dutch roll. 

Pars. 3.4.3, 3.4.6, and 
3.4.7 interpreted to mean 
C, <0 


*MlL-F-8785 recognizes augmentation -on and -off cases. This document deals only with 
unaugmented cases. 




TABLE 1 . - SUM-MAR Y OF STA TIC STABILITY CRITERIA (CONCL UDED) 


General form of 
static stability 
criterion 

Approximate or 
alternate form 

b 

Reference 44 
FAR - Part 25 

Reference 1 1 
British CAR, section D 

A 

O 

C„ < 0 or 
% 

Cd >0 

(No thrust effect) 

No direct requirement 

No direct requirement 

A 

O 

v° 

(No thrust e 

effect) 

Par. 25.177(c) interpreted 
to mean Cy^<0 

Par. 7.3 interpreted to 
mean Cy^ < 0 

*|J» 

A 

O 

Q _ >o 


No direct requirement 

No direct requirement 

o 

A 

N 

S *4 

CD CD 

V 

o 


Par. 25.17(a) interpreted 
to mean C n >0 

Par. 7.2 interpreted to 

mean C n >0 
"0 

V 

A 

O 

o 

V 

3 

u E 


No direct requirement but 
Pars. 25.173 & 25.175 
interpreted to mean C m ^<0 

Par. 2. 1 requires 
< + .05 

dC L 

o 

V 

si* 

Ci <0 
p 


No direct requirement 

No direct requirement 

9M y 

9q <0 

o 

V 

/ 


No direct requirement, but 
Par. 25.181 requires all j 

short periods to be heavily 
damped 

No direct requirement, but 
Par. 8.1 requires all short 
periods to be heavily 
damped 

o 

V 

si* 

v° 

- 

No direct requirement, but 
Par. 25.181 requires all 
short periods to be heavily 
damped 

No direct requirement, but 
Par. 8.1 requires all short 
short periods to be heavily 
damped 

V 

o 

w> 


Par. 25.175 (c) implies that 
violation is not allowed 

Par. 31.2 implies that 
violation is not allowed 

9M„ 

__JL < 0 

3v 

© 

V 


Par. 25.177 (b) interpreted 

to mean C« < 0 
*(3 

Par. 7.1 interpreted to 
mean C jj < 0 




The physical meaning of this criterion is that as a consequence of a side speed 
disturbance, v (along the y-axis), a force is generated ‘that tends to oppose v. The 
approximation, v** p V C j will be used. 



The physical meaning of this criterion is that as a consequence of a positive velocity 
disturbance, w (along the z-axis), a force is generated that tends to oppose w. The 
approximation w =» a y will "be used. 



The physical meaning of this criterion is that as a result of an angle of sideslip 
disturbance, p , the airplane weathercocks into the new relative wind. 



The physical significance of this criterion is that as a result of an angle-of-attack 
disturbance, a , the airplane weathercocks into the new relative wind. 


Sf> 




Criterion 


t 

An airplane is statically stable for a disturbance in roll velocity, p, if: 9M x /9p< 0 

( 102 ) 


The physical meaning of this criterion is that as a result of an increase in rolling 
velocity, p, a moment is generated that tends to oppose the increase in rolling velocity. 


Criterion 

An airplane is statically stable for a disturbance in pitching velocity q, if: 
9My/9q< 0 


(103) 


The physical meaning of this criterion is that as a result of an increase in pitching 
velocity, q, a moment is generated that tends to oppose the increase in pitching velocity. 


Criterion 

An airplane is statically stable for a disturbance in yawing velocity r, if: 

9M z /3r<0 (104) 


Tne physical meaning of this criterion is that as a result of an increase in yawing 
velocity, r, a moment is generated that tends to oppose the increase in yawing velocity. 

Under the definition of static stability used in this report, the partial differential 
9M y /3u (~C mu ) does not qualify as a static stability parameter. However, as will be 
shown, C mu has important consequences to longitudinal stability from the viewpoint of the 
pilot. In addition, in much of the literature, this parameter is identified with longitudinal 
stability. 
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It is noted that a positive sign of 3My/9u means physically that as a result of an 
increase in forward speed, the airplane noses up, which would tend to slow the airplane 
down because of the resulting drag increase plus the increase in giavitational pull a’.ong the 
body x-axis. Therefore, an airplane would have stable pitch moment versus speed behavior 
if: 9M v /8u>0. 

Under the definition of static stability used in tills report, the partial differential 
3M x /9v ( ~ Cj does not qualify as a static stability parameter. Nevertheless, this 
derivative has an important effect on stability and on handling qualities. 

The physical significance of this is that for a positive sideslip disturbance (nose left) the 
airplane tends to roll away from the disturbance, i.e., roll to the left. If the airplane rolls 
about its stability x-axis as a result of this, it is easily seen that this tends to diminish the 
effective sideslip angle. For this reason, some investigators identify Cg^ as a lateral stability 
parameter even though strictly speaking the derivative should not be considered as such. 


6.2 Dynamic £ lability Criteria 

Dynamic stability is defined as the tendency of the amplitudes of the perturbed 
motion of an airplane to decrease to zero or to values corresponding to a new steady state at 
some time after the disturbance has stopped. For example, consider an airplane disturbed in 
pitch from a steady-state flight condition. If the resulting perturbed motion is damped out 
after some time, the airplane’s motion becomes steady. If the new state is not significantly 
different from the original one, the airplane is called dynamically stable. The subject of 
dynamic stability, then, deals with the behavior of the perturbed motion of an airplane 
about some steady-state flight path. 

A dynamic stability criterion is defined as a rule oy whi'h perturbed motions are 
categorized as stable, neutrally stable, or unstable. 

In other contexts a dynamic stability criterion has been interpreted as a requirement 
for specific response characteristics or for meeting specific frequency-damping relations. 
This type of interpretation is embodied in the military specification for flying qualities (ref. 
10) and its proposed revision as documented in ref. 12. There are important connections 
between dynamic stability criteria (viewed as mathematical statements of stability) and the 
handling-qualities criteria of refs. 10 and 12. Therefore, where the need for physical 



interpretation of the stability criteria is established in this report, the connections with 
handling qualities are pointed out. 

Dynamic stability criteria are established covering the linear and nonlinear equations of 
motion of an airplane. These criteria apply to rigid, equivalent elastic, and completely elastic 
descriptions of airplanes provided the corresponding equations of motion are written in the 
form pertinent to the criteria. Table 2 presents a summary of dynamic stability criteria and 
their relationships to the various forms of the equations of motion. The arrangement of the 
equations of motion into the required forms was discussed in Sec. 5. Those combinations of 
criteria and equations which are most commonly used in airplane stability analysis are 
identified in table 2 "ith heavy lines. 

6.2.1 Characteristic equation methods. — When airplane dynamic behavior can be 
approximated by assuming that motion perturbations relative to a steady state are small, it 
is possible to reduce the equations of motion to a set of linear, second-order, ordinary 
differential equations vith constant coefficients. It is shown in app. A that in such a case 
these equations can be reduced to the following general form: 

{x} = [A]{x} (105) 


Taking the Laplace transformation and rearranging equation (105), it follows that 


[CsJ - [Al] {* s (S)}={x t (t o + j} 


(106) 


where: 

1. S = ct ± jw = complex frequency variable 

2. The subscripts s and t are used to distinguish between the functional relationships 
{X s } and {X t } 


The quantities oj and cej are, respectively, the real and the imaginary parts of the roots 
Sj of the characteristic equation: 


IICS] -(A ]|| = 0 


(107) 
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TABLE 2. - SUMMAR Y OF DYNAMIC STABILITY CRITERIA 
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Equation (107) follows from the condition that nontrivial solutions ({X s (S)} * 0) 
to equation (106) for {Xt(t 0 + )} = 0 are allowed only if equation (107) is satisfied. 

The following dynamic stability criterion can then be formulated. 


Criterion 

If the airplane equations of motion are linear and autonomous, then the airplane 
behavior is said to be: 


• stable, if the real parts of the roots of the characteristic equation are all negative, 

• neutrally stable, if there are one or more roots of the characteristic equation with 
zero real parts and the remaining roots all have negative real parts, and 

• unstable, if there is at least one root of the characteristic equation with a positive 
real part. 

(108) 


In most of the standard literature (refs. 4 and 13 through 17), dynamic stability of 
airplanes is treated from this viewpoint. 

The possibility exists for A in equation (105) to have elements that aie known 
functions of time. This occurs in steady climbs and dives when dynamic pressure is allowed 
to vary. In such a case, equation (105) is still linear but is called nonautonomous, and the 
equations assume the form : 

{x}= (A(t)l {x} (109) 


For this type of equation, no simple stability theory is known. (The more complicated 
stability theory of Lyapunov, which is discussed in par. 6.2.4 and in app. A, could be used.) 
However, the writers of this report feel that equation (109) is valid, although no proof has 
been found justifying the approach. The approach consists of applying the characteristic 
equation method to equation (109) with the following modifications. The characteristic 
equation considered is: 


HlA(t)] -XWI-O 


( 110 ) 
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The following dynamic stability criterion is postulated. 


Criterion 

When the real parts of the roots of the characteristic equation, equation (110), are 
negative for t = 0, as well as for t = t j, where tj is the practical limit of the time 
interval considered, then the airplane is stable in that time interval. 

(Ill) 


As stated, this criterion needs proof. The proof may lie in defining the quantity E as 
equal to the total kinetic energy of the airplane in the perturbed state and also defining E = 

dE/dt. In that case, E takes the place of the Lyapunov function. If criterion (111) is 

* 

satisfied it can be interpreted to mean, according to Lyapunov’s theorem 1, that E is 
negative at the beginning and end of the time interval. 

Whether or not in a practical case the dynamic stability criteria (108) are satisfied can 
be determined by solving directly for the roots of the characteristic equation. A technique 
for determining stability behavior from the characteristic equation without solving for the 
roots is known as Routh’s criterion. The reader is cautioned not to confuse this with what is 
sometimes called the Routh-Hurwitz criterion. Routh and Hurwitz both developed similar, 
but not identical, criteria. However, from the standpoint of calculations Routh’s is the more 
direct approach and is discussed in app. A. 

There are many techniques used in systems, analyses, and synthesis techniques (control 
theory) that may be applied to the perturbed airplane equations of motion of the form of 
equation (105). Some of the more widely used, for example Bode diagrams, Nichol’s charts, 
the Nyquist criterion, root locus plots, phase trajectories, etc., can be found in the literature 
(for example, refs. 18 through 22). Most of these techniques were generated for special 
types of problems and their use is restricted because of limitations imposed by the number 
of assumptions and/or the effort required in their application. 

Applications of the characteristic equations method are limited to linear differential 
equations of motion; they cannot be used where significant nonlinearities are encountered. 
However, when the equations of motion can be linearized, the characteristic equations 
method represented by criterion (108) is a most efficient technique for determining airplane 
stability behavior. In addition to determining the stability behavior, the roots of the 
characteristic equations can be used for other analyses. For example, the frequency and 
damping characteristics, imaginary and real parts of the roots, are used intensively in 





handling-qualities analyses and stability augmentation systems design. Examples of such uses 
may be found in refs. 4, 10, and 12 through 17. 

6.2.2 Time histories. - When the equations of motion of an airplane cannot be 
linearized, it has been common practice to base judgment of stability behavior on 
observation of traces of time history solutions of the equations of motion. 

The resulting traces are judged to determine stability. 


Criterion 


When the motions of an airplane following a disturbance from steady-state flight are 
determined by a time history, the behavior is: 

• stable, if if. 0 motions remain in proximity to the steady state, 

• neutrally stable, if the motions are undamped and oscillatory about some steady 
state. 


• unstable, if the motions diverge from the steady state either linearly, exponentially, 
oscillatorily, or in any combination. 


( 112 ) 


It should be recognized, however, that for nonlinear equations of motion there may be 
combinations of locally stable, neutrally stable, or unstable regimes. An example of what 
can be encountered is airplane pitchup above a certain angle of attack. If motions were 
investigated for small angles (caused by a small initial disturbance), stable motion would be 
observed. However, if a large disturbance were introduced which would force the angle of 
attack above the pitchup angle, divergence would occur. It is easily realized that one stable 
case does not guarantee airplane stability throughout its operating limits. Practically, this 
problem is resolved by running several time histories with various disturbance magnitudes. 

6.2,3 Energy decay methods. — A relatively new and unknown type of stability 
analysis is the energy decay method, which is discussed in refs. 23 and 24. The fundamental 
idea behind energy decay methods is that in dynamically stable systems, energy is being 
dissipated. In the case of linear differential equations with constant coefficients, it is 
possible to show the inverse, i.e., that if energy is being dissipated the corresponding system 
is dynamically stable. Extension to nonlinear equations of motion can be justified by 
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applying the Lyapunov stability theory discussed in app. A. The approach can be stated in 
the following steps. 

1. Derive expressions for the total perturbed energy, E, of the airplane. 

2. From (1) derive the AE required to make the airplane appear to be a conservative 
system in the first half cycle of oscillation, i.e., neutrally stable. 

3. Examine the sign of AE. 


The following dynamic stability criterion can then be formulated. 



A theoretical discussion on how to apply this criterion is given in ref. 24. However, 
because of algebraic complexities it is not considered practical to use the above criterion in 
cases involving nonlinear equations of motion. For these cases, Hahn (ref. 25) suggests an 
energy decay method based on an idea by Lebedev. This idea was further developed by 
Roskam (ref. 23). There, stability was connected with energy decay through the parameter: 



where T is the perturbed kinetic energy, ti is the beginning of a time interval during which 
the motion of the airplane is being studied, t 3 is the end of that time interval, and t 2 is the 
midpoint of that time interval. The criterion for stability in this case follows. 
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It is shown in ref. 23 that the condition F< 1, indicating stability, is satisfied in the 

case of stable, linear, small perturbation equations of motion. 

» 

» * 

The advantage of the latter criterion is that it applies to nonlinear equations of motion. 
A disadvantage is that either considerable numerical work or a computer program is 
required. 

Because energy decay methods have not been widely applied, their limitations, 
advantages, and disadvantages have not been assessed, although it is felt that there should be 
few limitations because of the general nature of the approach. However, for linear, 
autonomous, small perturbation equations of motion, this approach will probably prove to 
be less efficient than the characteristic equations method. 


6.2.4 Lyapunov’s method. — The time history method was suggested as a way to 
determine the stability behavior of an airplane when the equations of motion are nonlinear. 
However, with the time history method it is necessary to solve the equations of motion. 
Lyapunov has devised a stability theory for both linear and nonlinear differential equations 
of perturbed motion that obviates the necessity to solve these equations. 


Lyapunov’s stability theory has received little attention from airplane stability and 
control engineers. For this reason, an introduction to this theory as well as some pertinent 
definitions and theorems are given in app. A. The potential applications of the analysis 
techniques devised by Lyapunov, and by others who have followed his approach, are quite 
numerous becuase of the generality of the approach. Rather than solving any particular 
problem, Lyapunov realized that the stability of dynamic systems (including moving bodies) 
was a problem of studying the behavior of differential equations in general. He devised two 
classes of approach, one for equations whose solutions are known functions of time and 
another for equations of motion written in the perturbation form. The first approach, using 
known solutions, is similar to the use of the stability criterion for time histories. 


The second approach, called the “direct” or “second” method of Lyapunov, requires 
choosing a “Lyapunov function” and relating the behavior of this function to the behavior 
of the differential equations of motion. Owing to a theorem attributed to Zubov, 
Lyapunov’s direct method becomes a particularly attractive approach to the problem of 
nonlinear airplane stability behavior. Because of its similarity to the more familiar 
characteristic equation approach it is felt that Zubov’s theorem should appear as a logical 
extension of that approach. In fact, as shown in app. A, it is possible to prove criterion 
(108) (stable roots for characteristic equations) using Zubov’s theorem for linear, 
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autonomous equations. However, the application of Zubov’s theorem would be more useful 
for nonlinear equations. . 

The large perturbation equations of motion of an airplane can be written in the form: 

{x} = (F ({x}, t)]{x} (116) 

Also, nonlinear, small perturbation equations with nonlir ’ar, aerodynamic, cross-coupling 
terms, can be written in this form. Before stating the stability criterion for these nonlinear 
equations, however, some definitions are required. First, the equation 

|||[(F({x r }, t R )] T + [F ({ X R }> ty]]- A[l]|| = 0 (117) 


will be called the “quasi-characteristic equation” where {xp} and tp are defined as values 
belonging to a representative set of x and t. Next, by a “representative set,” the following is 
meant. For given initial disturbances, the solutions to the equations of motion (116) yield a 
time sequence of values of the motion variables {x} . In most practical cases the engineer 
will have an idea of the practical limits of the perturbed motions the airplane can 
experience. In other words, the engineer can make a reasonable estimate of the “cylindrical 
neighborhood” surrounding the time axis, within which the motion takes place. Choosing 
discrete values of {x} and t, called {xp} and tp, within practical limits related to the 
steady-state flight condition in accordance with these ideas, generates a representative set of 
{ x } and t. Analogous to this choice of a representative set is the selection of combinations 
of Mach numbers, dynamic pressures, angles of attack and angles of sideslip for which wind 
tunnel data is to be obtained or stability is to be assessed in the usual analytical approach. 

Finally, the eigenvalues, \ , that will satisfy equation (1 17) are called the eigenvalues of 
the quasi-characteristic equation. 

Using the above definitions, the application of Zubov’s theorem as a dynamic stability 
criterion is postulated as follows. 


Criterion 

If the eigenvalues of the quasi-characteristic equation are nonpositive (s 0) for each 
{Xp} and tp in a representative set of {x} and t, then the airplane is stable. 


( 118 ) 




Unlike the other stability criteria presented in this section, this criterion has no neutral 
or unstable counterparts. It is shown in app. A that Lyapunov’s direct method has certain 
limitations and that the existence of positive eigenvalues does not necessarily imply 
instability. In other words, criterion (118) is necessary but not sufficient. 

It is emphasized that the use of Zubov’s theorem as a basis for determining stability has 
limitations and disadvantages. Particularly important is the consideration that proving 
Zubov’s theorem requires the use of a particular Lyapunov function, which may lead to very 
rough stability analyses. Another disadvantage is loss of physical feel for the problem until 
familiarity with and understanding of the direct method are achieved. Where the engineer 
can visualize the predicted motions, the time history approach has a distinct advantage. 

The particular method presented here, using Zubov’s theorem, has no restrictions on 
the types of perturbation equations to which it can be applied. However, the. criterion only 
pertains to dynamic stability and, as stated above, there are no neutral or unstable 
counterparts. Also, it will not always predict stability for stable airplanes. This approach has 
not been sufficiently explored to adequately assess its value. 
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7. METHODS FOR DETERMINING STABILITY DERIVATIVES 


7.1 General Considerations 

Aerodynamic derivatives are used to relate changes in the aerodynamic forces and 
moments to changes in the airplane’s attitude, motion, and shape due to elastic deformation 
and control deflections. Those derivatives associated with airplane translation and rotation 
are generally called the airplane’s stability derivatives. Here, those associated with elastic 
deformation are also included in that classification. 

The stability derivatives must be estimated in order to evaluate the stability 
characteristics of an airplane configuration. Several methods for obtaining estimates are 
usually used in combination. These methods divide roughly into three categories: (1) 
estimates based .on numerical solutions to the equations of fluid dynamics, (2) estimates 
based on semi-empirical handbook data, and (3) estimates based on experimental wind 
tunnel data. The objective here is to describe the best available analytical techniques. These 
techniques are then evaluated by comparing computed values for stability derivatives with 
their values measured in wind tunnel and flight tests and by carefully investigating the 
approximations involved in the theoretical and handbook methods. 

There are 24 stability derivatives associated with the translational and rotational 
degrees of freedom that are usually regarded to have some importance to the stability 
characteristics of large, flexible airplanes. They are listed in table 3 along with an evaluation 
of their relative importance. The degree of confidence that can be placed on their predicted 
values v/hen wind tunnel measurements are used in conjunction with handbook and 
theoretical techniques is also shown. This table reflects a consensus existing among 
experienced stability and control engineers at the time of writing. It illustrates the need for 
improved techniques for estimating lateral-directional stability derivatives. 

The advent of large digital computers has greatly enhanced the theoretical approach to 
estimating stability derivatives. One conclusion of this study is that digital computer 
programs based on solutions to the equations of fluid dynamics constitute the best method 
for estimating stability derivatives. These programs can incorporate empirical corrections in 
somewhat the same manner as the handbook methods. They can also use data from wind 
tunnel pressure models and the results of flow visualization studies as empirical corrections 
in a way that is impossible in the current handbook methods. 

The computer programs generally use influence coefficient theory. Aerodynamic 
influence coefficients give the change in pressure at a small region of the airplane’s surface 
due to a change in surface incidence at that and any other small region of the airplane’s 
surface. Structural influence coefficients give the change in incidence and position at a small 
region due to a change in force applied to that and any other region. 



TABLE 3- RELATIVE IMPORTANCE AND PREDICTION ACCURACY OF 
STABILITY DERIVATIVES 



Longitudinal 


Lateral-directional 

Derivative 

Relative 

importance 

Estimated 

accuracy* 

Derivative 

Relative 

importance 

Estimated 

accuracy* 

C L a 

Primary 

Good 

c„ 

n 0 

Primary 

Good 

C m a 

Primary 

Good 

C i/3 

Primary 

Good 

c m 

m q 

Primarv 

Good 

% 

Primary 

Good 

r 

m u 

. Primary 

Good 

S 

Primary 

Fair 

c m . 

m a 

Primary 

Fair 

S 

Primary 

Poor 

C D 

Primary 

Good 

c* 

r 

Primary 

Poor 

Cr\ 

u a 

Secondary 

Good 

C Y/3 

Primary 

Good 

C L 

u 

Secondary 

Good 

\ 

Cy p 

Secondary 

Poor 

C L& 

Secondary 

Poor 

S 

Secondary 

Poor 

C L 

L q 

Secondary 

Good 

C V r 

Secondary 

! 

Fair 

Cd* 

u a 

Minor 

Poor 

% 

Minor 

Poor 

Cr\ 

u q 

i 

Minor 

Fair 

| 

Minor 

Poor 


""Estimated prediction accuracy assumes use of theoretical, handbook, and 
wind tunnel data. 
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The influence coefficients can be used in a formu'ition of the stability derivatives for a 
rigid airplane, an equivalent elastic airplane, or for the elastic degrees of freedom of an 
airplane represented as completely elastic. The aerodynamic^ influence coefficients are based 
on either lifting surface theory or lifting line theory, although lifting surface theory is the 
more general and is not restricted in application to high aspect ratio wings as Is lifting line 
theory. 

As noted, thu aerodynamic influence coefficients can either be combined with the 
structural influence coefficients to formulate equivalent elastic stability derivatives or be 
used by themselves to formulate rigid stability derivatives. The ratio of these two stability 
derivatives or the difference between them can then be used to correct handbook estimates 
of rigid stability derivatives for the effects of elasticity. These procedures can also be used to 
correct derivatives measured in the wind tunnel using rigid or nearly rigid models. It 
represents the best available method for correcting rigid stabiii y derivatives. 

The most widely known source of information for estimating stability derivatives is the 
USAF Stability and Control Handbook (ref. 6). It was the primary handbook source used in 
the evaluation. This handbook, as well as all others, utilizes empirical and theoretical data 
determined over a range of certain parameters associated with airplane configuration, e.g., 
aspect ratio and tail volume, and with the airplane reference flight condition. The stability 
derivatives of a particular airplane are found by interpolation of the handbook data. 

A serious deficiency in the handbook is that the effects of structural flexibility are 
either not accounted for or are accounted for only through coarse approximations. The 
handbook methods are not sufficient in themselves to treat flexible, low aspect ratio wing 
and tail surfaces. 

Tables 4 and 5 summarize the results of the evaluation of methods for calculating 
stability derivatives. The estimates based on ..'ting line theory (of. 26) contain empirical 
corrections and were made using a computer program that is typical of those available 
throughout the industry. Those based on lifting surface theory (ref. 27), however, were 
made without empirical corrections. It is important to note this because the lifting surface 
technique used in this evaluation was developed for wing optimization studies. Leading edge 
suction is poorly represented in this program; hence, induced drag and yaw due to roll are 
poorly estimated. The program will only represent an airplane in symmetric flight or a flat 
airplane in roll. In addition, a first approximation to derivatives that involve unsteady 
aerodynamics is not an immediately accessible feature. Ail of these aspects reflect against 
the lifting surface estimates, although they may be overcome in future programs. 

Tphles 6, 7, and 8 summarize the limits of applicability ol the techniques for obtaining 
stability derivative estimates. Tables 6 and 7 show the applicability of techniques using 
aerodynamic influence coefficients-table 6 concerns techniques based on lifting surface 
theory and table 7 those based on lifting line theory. The applicability of handbook 
methods is shown in table 8. 
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TABLE 4.- LONGITUDINAL STABILITY DERIVATIVES 




1 

■ 

Lifting 


Handbook 

Handbook + 




■ 

line 





lifting 




(computer,- 

(computer) 




surface 


Derivative 

Rigid and 

Rigid and 









equiv. elas. 

equ«v. elas. 

Rigid only 

Equiv. elas. 



■■ 

Sup 

Sub 

Sup 

Sub 

Sup 

Sub 

Sup 



707 

SST 

SST 

707 

SST 

SST 

707 

SST 

SST 

707 

SSI 

SST 


C D 

a 

B 

B 

P 

B 

P 

■ 

P 

P 

■ 

fl 

■ 

■ 

Ot 

Cl « 

F 

D 

G 

G 

G 

■ 

G 

P 

G 

G 

fl 

P 


C m a 

G 

p 

F 

F 

F 

■ 

P 

P 

G 

P 

II 



C D* 

a 

fl 

fl 

■ 

B 

■ 

■ 

P 

P 

P 

B 



• 

a 

c L . 

B 

■ 

■ 

B 

fl 

■ 

P 

P 

II 

fl 




c 

m a 

fl 

■ 

■ 

B 

■ 

fl 

P 

P 

P 

B 




C D 

B 

II 

P 

fl 

■ 

■ 

B 

■ 

■ 

B 



u 

Cl 

L u 

fl 

G 

G 

B 

■ 

■ 

p 

P 

P 

B 




c m 

m u 

G 

G 

G 












P 

P 

P 







B 

■ 


q 


F 

F 

F 

B 

■ 

■ 

p 

II 

P 

B 

p 

P 


C m q 

F 

F 

F 

fl 

■ 

■ 

p 

P 

P 

p 

p 

P 


G (good)— method compares favorably with experiment, with some exceptions. 
F (fair)-method compares favorably with experiment, with exception of some 
M, A , etc. 

P (poor)— method does not compare favorably with experiment. 

[BLANK] not calculated. 
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TABLES.- LATERAL-DIRECTIONAL STABILITY DERIVATIVES 




Lifting surface 
(computer) 

Lifting, line 
(computer) 

Handbook 

Handbook + 
lifting surface 

Derivative 

Rigid and 
equiv. elas. 

Rigid and 
equiv. elas. 

Rigid only 

equiv. elas. 



Sub 

Super 

Sub 

Super 

Sub 

Super 

Sub 

Super 



707 

SST 

SST 

707 

SST 

SST 

707 

SST 

SST 

707 

SST 

mm 


% 







F 

F 

F 

fl 



0 

C *0 







n 

fl 

P 

■ 




% 







p 

' 

F 

P 

■ 




m 







i 

P 

P 

■ 



0 

C*0 







p 

P 

P 

■ 




c . 
n 0 




■ 



p 

B 

P 

■ 




a 




1 



B 

fl 

■ 

P 

fl 

D 

P 

% 

F 

F 

F 

F 

F 


B 

P 

■ 

P 

P 

P 


% 

■ 

■ 

■ 

■ 

1 


D 

P 


D 

P 

P 


m 

■ 

■ 


1 

■ 


B 

P 

P 

■ 

■ 


r 

% 

■ 

■ 


■ 

■ 


p 

P 

P 

B 

■ 

■ 

■ 


s 

■ 

■ 

■ 

1 



p 

P 

P 

■ 

■ 

iHH 


G (good)-method compares favorably with experiment, with some exceptions. 
F (fair)— method compares favorably with experiment, with exception of some 
etc. 

P (poor)-method does not compare favorably with experiment. 

[BLANK] —not calculated. 
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TABLE 6. - AFPLIC ABILITY AND LIMIT A TIONS OF LIFTING SURFACE THEOR Y 




Thiit-body 
formulation, 
thickness ratio 
6« c, b/2 

Non-coplanar, thick-body formulation 

Deriva- 

tive 

Wing 

Wing-body 

Non- 

coplanar 

surfaces 

Thick 

surfaces 

■ 

Wing-body 
Horiz-tail 
Vert, tail 

Wing+tail 

dihedral 




T 



T 



T 



T 



T 



T 







S 

R 

S 

s 

R 

S 

S 

R 

s 

S 

R 

S 

s 

R 

S 

S 

R 

S 

S 


S 



U 

A 

U 

U 

A 

U 

U 

A 

u 

U 

A 

U 

u 

A 

u 

U 

A 

u 

U 


u 



B 

E 

P 

B 

El 

P 

B 

□ 

p 

B 

□ 

P 

B 

□ 

p 

B 

N 

D 

B 


p 

a 

C D a 

C L 

r ^a 

m<Y 



B 

fl 

| 

H 

fl 

1 

B 

fl 



H 


H 

H 



N 

N 

N 

m 

a 

C Da 
Ct v 
r L c* 

ntQ, 

1 


B 

fl 


B 

fl 

1 

fl 

fl 



fl 


fl 

1 



N 

N 

N 

u 

C D U 

CL U 

C m u 

D 


B 

fl 

1 

fl 

fl 

I 

B 

D 



H 


fl 

H 



N 

N 

N 

q 

CD q 

Cl2 

Cm’ 

fl 


B 

B 


B 

fl 

1 

fl 

fl 



fl 


fl 

fl 



N 

N 

N 

0 

np 

B 


B 

B 


fl 

fl 

l 

fl 

fl 



1 


fl 

H 



H 

1 

fl 

0 

pi 

Pi 

«0 

B 


B 

fl 


B 

1 

l 

B 

1 



1 


1 

1 



1 

1 

1 

1 

C Yp 

p* 

n P 

B 


B 

fl 


B 

D 

l 

fl 

fl 



1 


1 

1 



1 

1 

1 

1 

C Y r 

C *r 

Cn F 

B 


B 

H 


fl 

1 

l 

1 

1 



1 


1 

1 



1 

1 

1 


Ability to calculate various 

flow conditions 

Flow 

Capability 

Compressible 

Yes 

Steady 

Yes 

Quasi-steady 

Yes 

Unsteady 

Yes 

Viscous 

No 

Separated 

No 


• Capability exists or can be 

developed from existing theory. 

N Not applicable to first-order 

approximations of this derivative. 

Blank No capability exists. 

? Capability may exist, but further 

development is required. 









































TABLE 7 .- APPLICABILITY AND LIMITATIONS OF LIFTING LINE THEORY 


Deriva- 


Wing-body 
horiz. tail 

Non- 
coplanar 
sv -faces 

Wing-body 
horiz-tail 
vert, tail 

Wing-tail 

dihedral 

Thick 

surfaces 

Thick 

body 

tive 

S 

T 

R 

S 

S 

T 

R 

S 

S 

T 

R 

S 

S 

T 

R 

S 

S 

T 

R 

S 

S 

T 

R 

S 

S 

T 

R 

S 



U 

A 

U 

U 

A 

u 

U 

A 

U 

u 

A 

U 

U 

A 

U 

U 

A 

U 

U 

A 

U 



B 

N 

P 

B 

N 

p 

B 

N 

P 

B 

N 

P 

B 

N 

P 

B 

N 

P 

B 

N 

P 

a 

Crv 

c? 01 

c„ a 

m Q 

• 



• 



• 






N 

N 

N 







• 

a 

Cr\. 

r ° 
C L* 

c m a . 




• 



• 






N 

N 

N 







u 

r Du 

c Lu 

m u 

• 



• 



• 






N 

N 

N 







q 

C n 

c£j 

C m q q 

• 

* 


• 



? 






N 

N 

N 







0 

p 0 

np 

• 






? 






? 









* 

0 

pi 

«0 

? 



? 



? 















P 

Cv 

c> 
c p 

^n p 

? 



? 



? 















r 

c 

n r 

? 



? 


* 







I 


i 








Ability to calculate various flow conditions 

Flow 

Capability 

Compressible 

Yes 

Steady 

Yes 

Quasi-steady 

Yes 

Unsteady 

Yes 

Viscous 

No 

Separated 

No 


• Capability exists or can be 

developed from existing theory. 

N Not applicable to first order 

approximations of this derivative. 

Blank No capability exists. 

? Capability may exist, but further 

development is required. 
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TABLE 8. -APPLICABILITY AND LIMIT A TIONS OF THE HANDBOOK METHODS 


Deriva- 

Wing 

Wing-body 

Non- 
cop Lanar 
surfaces 

Thick 

surfaces 

Thick 

tjpdy 

t 

W-B-Ht 
vert, tail 

Wing-tail 

dihedral 

five 

S 

T 

R 

S 

S 

T 

R 

S 

S 

T 

R 

S 

S 

T 

R 

S 

S 

T 

R 

S 

S 

T 

R 

S 

S 

T 

R 

S 



u 

A 

u 

U 

A 

U 

U 

A 

U 

U 

A 

U 

U 

A 

u 

U 

A 

u 

U 

A 

U 



B 

N, 

p 

B 

N 

P 

B 

N 

P 

B 

N 

P 

B 

N 

P 

B 

N 

p 

B 

N 

P 

a 

Cd 

c? a 
c a 

Lm a 

• 

• 

• 
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• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

N 

N 

N 


Cr>* 
r u a 
w - 










• 


• 










+ 

a 

• 

• 

• 

• 

• 

• 

• 

• 

• 


• 


• 

• 

• 

• 

• 

• 

N 

N 

N 



• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 




u 

Cl u 

C u 
m u 



















N 

N 

N 




* 
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C A 

m q 

• 
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• 

• 

• 

• 

• 

• 

• 


• 


• 

• 

• 

• 

• 

• 

N 

N 

N 


• 

• 

• 

• 

• 

• 

• 

• 

• 


• 


• 

• 

• 

• 

• 

• 




ft 

r Y S 

£*<> 

*P 

• 

• 

• 

• 

• 

• 

i 

— . .1 

• 

• 

• 

• 

• 

• 

i 

• 

• 

• 

• 

• 

• 

i 

i 

• 

• 

• 

• 

P 

C Y a 

cS! 

C*J 

n P 






1 

| 




• 

• 

• 

• 

• 

! 

! 

I 

• 

• 






Cy p 






? 









9 




? 

? 

? 

p 

c i 

c *p 

• 

• 

• 



• 




• 

• 

• 



• 

• 

• 

• 

• 

• 

• 




I 



• 









• 






• 


S Y r 

? | 


? 1 



- 









1 




9 


? 

r 
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• 


• 







• 

• 

• 



i 

. i 

• 

• 

• 

• 


• 


S 

• 

• 

• 











1 

i 




? 


? 


Ability to calculate various flow conditions 

Flow 

Capability 

Compressible 

Yes 

Steady 

Yes 

Quasi-steady 

Yes 

Unsteady 

Yes 

Viscous 

No 

Separated 

No 


• Capability exists or can be 

developed from existing theory. 

N Not applicable to first order 

approximations of this 
derivative. 

Blank No capability exists. 

? Capability may exist, but 

futther development is required. 




The. preference for lifting surface theory expressed in this report is heavily supported 
by the results presented in these tables. Three classes of capability are shown. That which 
exists or can be acquired readily is shown by a dot. That which is within the fundamental 
limitations of the technique is shown by a question mark. A blank space indicates capability 
which exceeds fundamental limitations of the technique. 

Estimation techniques using lifting surface theory clearly have the greatest basic 
potential. When this theoretical capability is coupled with empirical corrections to the 
lifting, surface technique, the preference is further justified. Empirical correction methods 
may be utilized in the aerodynamic influence coefficient methods in a very efficient and 
rational manner. This feature has been incorporated into existing computer programs based 
on. lifting line theory but is currently in only a rudimentary stage of development in existing 
programs based on lifting surface theory. 

The formulation of stability derivatives in terms of aerodynamic influence coefficients 
involves the theory of small perturbation inviscid fluid flow. A general formulation of the 
aerodynamics of stability and control based on a mathematical theory of fluid dynamics 
appears nowhere in the literature. Many specialized investigations may be found, and 
examples appear in the works by Miles (ref. 28), Chester (ref. 29), Ward (ref. 30), and Van 
Dyke (ref. 31). The only general development found in the literature is that by Bryson (ref. 
32) t which applies to slender rigid bodies. However, Bryson’s development is valid for 
slender bodies only, to the exclusion of wing-body combinations with aspect ratios larger 
than 2 or 3, since his formulation is based on the constant-density, cross-flow assumption. 

It is also important to note that the assumptions applied to the conservation laws of 
ref. 33 to develop the inviscid fluid dynamic equations do not apply to hypersonic flow. 
Precise limits of applicability cannot be placed upon these equations, because the onset of 
hypersonic flow depends upon the shape of the body (ref. 34). Hypersonic flow 
characteristics are observed at Moo = 3.0 on some blunt bodies, but not until Moo = 10.0 on 
some thin bodies. For supersonic jet-type airplanes, which are relatively slender, important 
hypersonic flow effects do not occur below Moo = 5.0. Therefore, the flow equations used in 
this report are assumed valid up to that speed for such slender configurations. 

A formulation of the stability derivatives in terms of lifting surface aerodynamic 
influence coefficients is presented in par. 7.2. The basis for that formulation is described in 
physical terms in par. 7.2 and is developed in detail in app. B. Summarizing, the pertinent 
equations for the stability derivatives are written in terms of five aerodynamic matrices, 
[Aj], [A2I , [A 3 j , [Aq], and [A 5 ], which are derived from the equations of fluid 
dynamics as described briefly in par. 7.2 and in detail in app. B. Their precise form depends, 
'i pan, on the particular aerodynamic influence coefficient theory chosen. The equations in 
•err. i ' the aerodynamic matrices are as follows. 




The derivatives contained in the boxes above are lateral-directional; the rest are 
longitudinal. The aerodynamic influence coefficients cannot be used to compute drag due to 
viscous forces. However, those drag derivatives which are for the most part due to 
nonviscous forces can be computed. The values of Cj) u , and Cq can be computed 

for use in obtaining elastic corrections to rigid values obtained separately from a computer 
program system. 


Equations (119) and (120) contain certain approximations of the aerodynamic forces 
in order to be consistent with the small perturbation equations of motion. The 
lateral-directional unsteady derivatives Cy^ , Cj • , C n ^ contain the most severe approxi- 
mations and assumptions. The longitudinal steady derivatives appearing in equation (119) 
contain the least severe approximations. Besides those assumptions previously listed in 
deriving the small perturbation equations of motion, assumptions (DU) arR l (A8^,are 
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applied in calculating the longitudinal stability derivatives. The evaluation of the 
lateral-directional stability derivatives depends upon . Unsteady stability derivatives 


are based on approximation (A9 j in lieu of (A8 


Additional stability derivatives are required which correspond to aerodynamic forces 
generated by the dynamic-elastic motion, i.e., motion in the elastic degrees of freedom. 
Stability derivatives corresponding to the i 1 * 1 elastic degree of freedom are 
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'D^ 

1 


yj 

C Lv 

1 

i 

C m, 

1 


. t?) T [A 5 ]{ ei } 


“r 


where { 0 j } is the i in elastic free vibration modal shape. 


( 123 ) 


The derivatives appearing in boxes are lateral-directional and are nonzero if i 
corresponds to an antisymmetric mode shape and zero if i corresponds to a symmetrical 
mode shape. The remaining derivatives are longitudinal and are zero or nonzero conversely 
with tiie lateral-directional derivatives. 


Stability derivatives for changes in generalized aerodynamic forces in the elastic degrees 
of freedom due to changes in the motion velocity components are given by 



T 


-{.jfV,] 


( 124 ) 


Stability derivatives for changes in generalized aerodynamic forces due to rates of 
change in incidence and sideslip are: 



( 125 ) 


Stability derivatives relating changes in generalized aerodynamic forces to changes in the 
generalized coordinates and velocities in the elastic degrees of freedom are: 
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{ C J 1 }'i*j} T |A 8»*i> 

Kl’W* ,A 4H»i} 


(126) 


(127) 


Stability derivatives relating changes in generalized aerodynamic forces to changes in 
the generalized accelerations in the elastic degrees of freedom are: 

{ C ^{0j} T t A 5 1 {«,.} (128) 


The derivatives in equation (128) may be computed on the basis of aerodynamic 
influence coefficients for reduced frequencies that are small and for symmetrical mode 
shapes. 

In summary, equations (119) through (128) contain all the stability derivatives for a 
symmetric, elastic airplane. If the appropriate aerodynamic influence coefficients are 
available, then all of the stability derivatives may be computed. 

The stability derivatives computed on the basis of aerodynamic influence coefficients 
represent a valid first-order approximation for the pressure coefficient of an inviscid, 
isentropic fluid flow (ref. 30). In addition, it is assumed that the unsteady effects are due to 
motion of sufficiently slow variation that the reduced frequency is of zero order of 
magnitude (ref. 28). 


7.2 Theoretical Methods 

Stability derivatives relate changes in the aerodynamic forces and moments on an 
airplane to changes in airplane attitude and motion. This section describes a basis for 
computing the stability derivatives from the mathematical theory of fluid flow. The 
aerodynamic influence method is used in the development. Results sought are the 
aerodynamic matrices [ A j ] , [A 2 ], (A 3 ], [A 4 ], and [A 5 ], which appear in the matrix 
formulations of the equations of motion appearing in Sec. 5. The detailed developments are 
contained in app. B. This section merely summarizes and discusses the physical aspects of 
the problem, but proceeds from rather fundamental concepts. 
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.rplane in flight may be looked upon as a solid body in motion in a fluid. Its 
motion disturbs a portion of the atmosphere. However, its shape md motion are usually 
such that the perturbation of the disturbed fluid is very flight. In fact, the magnitude of 
perturbation velocity of the fluid, when divided by the magnitude of the airplane’s velocity, is 
of order of magnitude zero except in the vicinity of stagnation points. Also, the viscosity of 
the atmosphere is so small that viscous forces are negligible in comparison with dynamic 
forces in the fluid, except in a thin layer of flow at the airplane’s surface and in its wake. 
These conditions allow small perturbation im'iscid flow theory to provide an accurate 
theoretical representation of the flow about an airplane. 

As shown in app. B, as well as in most textbooks on aerodynamics, in the theory of 
small perturbation inviscid flow the principles of conservation of mass and momentum and 
the equation of state are reduced to a single linear partial differential equation. This flow 
equation contains a single dependent variable, perturbation velocity potential <p . The flow 
equation is written in terms of a fluid axis system whose coordinates, along with time, are 
the independent variables in the flow equation. 

The velocity of the fluid is given by V = U(i + V <t> ). This is the velocity of the fluid 
particles measured relative to the fluid axis system and U is the component, along the x-axis 
of the fluid axis system, of the velocity of the airplane relative to fixed space. 

The perturbation velocity potential must be such that it satisfies the flow equation and 
such that the velocity field V(x,y,z,t) does not produce flow through the surface of tne 
airplane. This last requirement is called a boundary condition. An additional requirement of 
the perturbation velocity potential follows from the vortex laws of Helmholtz. The wake 
behind the airplane must be such that they are satisfied. However, the pressure in the fluid 
calculated from Bernoulli’s equation must be continuous across the wake. Discontinuities in 
the pressure field in the fluid can only exist across solid r odies. This is all that will be said 
about the flow problem here; for a more detailed description, the reader may refer to app. 
B. The special form appropriate to the question of aerodynamic forces in stability and 
control considerations will now be introduced. 

The boundary condition at the surface (see app. B) may be written as 

n.(V s - H) = 0 (129) 

where n is a vector normal to the surface of the airplane. The velocities V $ and 'd are the 
velocities of fluid particles and elastic particles of the airplane at the surface. They are 
measured relative to the stability axis system that rotates and translates with the airplane 

<*s> Y s> Z s>- 



The stabi’ity axis system moves relative to the fluid axis system with the translational 
velocity = Vj^+ Wk^ and rotational velocity te = pi s + Qj s + RlT s . Thus, the two fluid 

velocities V and V. at the surface of the airplane are related as 

* • 

k 

V — V + Vrj + ux ? (130) 

S X\ S 

where 7^ is the position vector of a point on the surface relative to the center of gravity, the 
origin of the stability axis system. 

Using equation (130), the boundary condition at the airplane surface becomes 


n • [uT + UV0-(Qz s 


R V i s 


- (V + Rx 


pz s )j s -(w+py s 


-Qx ,) k, -a]= 0 (131) 


where it is important to note that some quant'ties are described in the fluid axis system 
Wiiile others are described in the stability axis system. Also, each of the velocity 
components appearing in equation (131) is treated as an independent variable in the 
equations of motion, Sec. 5. Further, the perturbation velocity V0 must be of order of 
magnitude zero. Thus, each individual term involving a velocity component when divided by 
U must be of order of magnitude zero for a valid approximation based on small 
perturbation, inviscid flow theory. The approximation in the small perturb tion equations 
of motion becomes linked to those in the flow theory through this boundary condition. 


In the reference motion of the airplane, which is steady, the boundary co^’ T. 
reduces to 


n, 


ni • V0 t • [U 1 i s + (Q.z s - V.) i, + (V, + V,) i, 


' Q i x Sl i S +3 d 


(132) 


In the disturbed motion, the airplane may be so orient d relative to imed space that 
the Euler angles are perturbed as 0,0,0 from those of the reference motion. These 
variables are treated as independent, so that their perturbations are independent of the 
velocity perturbations. The flow problem must be reorien.cd with the -.Irplane. The only 
change in the boundary condition resulting from the orientation perturbation arises from 
the elastic deformation perturbation due to the perturbation of the gravity forces. Thus, 
there is a change in the direct; l of the normal vector iT relative to the stability axis system, 
but there is no other change. 
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Neglecting products of perturbation quantities anu using the reft 'ence motion 
boundary condition, the perturbed boundary condition is found to be 
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The velocity potential o p is the perturbation to the perturbation velocity potential in the 
reference motion such that $> = 0j + 0 p . The vector fiTTj is the change in the vector 
normal to the surface caused by the change in elastic shape of the airplane from reference to 
disturbed motion. To a first-order approximation, it is given by the elastic rotation at the 
surface ©g as 

(134) 

and Gg is due to the perturbed gravitational, inertial, and aerodynamic forces. It must be 
found from the i r* *al equilibrium equations of Sec. 5, as must the perturbed elastic 

displacement rate, u_, 

The velocity component Un*V0 is the perturbation to the flow from its freestream 
direction that is required to satisfy the slope boundary condition at the airplane’s surface. In 
general, the perturbed flow may be unsteady, but the angle between th^ direction of the 
streamlines of the undisturbed flow and the streak lines of the disturbed how is defined by 
tan = TT*V0 at the airplane’s surface. This is the local incidence angle at any instant of 
time and because it is small 


& « n* V 0 (I 35 ) 

By definition of the reference motion, the incidence angle does not change with time 
in the reference motion, so that 


fi~ ”i • 


V9> 1 


(136) 


is ihe angle between unperturbed and perturbed streamlines. The perturbation, local 
incidence angle 



( 137 ) 
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P 


is in general a function of time. 


The local flow incidence angles must be related to aerodynamic pressure at the 
airplane’s surface before the above results may be used to generate stability derivatives. 
Surface pressures are readily computed from the velocity potential using an appropriate 
small perturbation form of Bernoulli’s equation as shown in app. B. Thus, the problem 
centers on finding perturbation velocity potentials which satisfy the flow equation and the 
boundary conditions. 


The most widely used means for finding perturbation velocity uses the linearity of the 
flow equation in its small perturbation form. Elementary solutions to the flow equation are 
obtained. These are perturbations to uniform flows produced by sources and doublets. Their 
distribution throughout the flow and their strengths are so adjusted that the uniform flow is 
perturbed into conformance with the airplane’s shape and motion. I he sources and doublets 
are located in the volume that is interior to the airplane and doublets are distributed in the 
wake. In lieu of the solid body, an interior flow is formed; the boundary between this 
interior flow and the exterior flow is made to conform to the airplane surface. Further, the 
circulation required by the Kutta condition is imposed and the wake is made to satisfy the 
vortex laws. 


The linearized flow equation is essential to this type of flow representation. The linear 
combination of the eflects of the sources and doublets is only possible if the flow equation 
i a linear partial differential equation. 

In practice, the flow boundary condition cannot, in general, be satisfied at every point 
of the airplane’s surface. A finite number of surface points are chosen as control points, and 
the boundary condition is satisfied at those points to within some cnosen numerical 
accuracy. At points in between, the boundary condition may or may not be satisfied. 
However, when the number of control points is large and the surface is a smooth 
aerodynamic surface, the approximation may be of high quality. 

Additional approximations may be introduced. These aie associated with a lineariza- 
tion of the boundary conditions in a manner consistent with the linearization of the flow 
equation. As shown in app. B, for thin bodies such as wings and tail surfaces the thickness 
may be considered separately from the incidence. For slender bodies such as fuselages the 
axial component of flow may be considered separately from the cross-flow component. 
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The thickness of thin bodies may be represented by a distribution of sources at a 
middle surface, and incidence is represented by distributed doublets. For the axial flow over 
slender bodies, line distributions of sources are placed on the fuselage axis. The effect of 
cross flow is met by line doublets along the fuselage axis. 

There are several schemes used in distributing the sources and doublets. They may have 
specified distributions over small regions called panels. In this case their strengths may be 
constant or have some specified variation over the surface of the panel. The distributed 
doublets on the panels are vortex sheets. In the method by Woodward, designated as 
computer program TA-67A in this report (ref. 35), the Kutta condition is satisfied along 
with the vortex laws on each panel (see app. B). This method does not require separate 
treatment of the wake. 

An alternate scheme for treating the vorticity distribution on thin bodies utilizes a 
vorticity distribution represented by a series of loading functions that vary over the entire 
middle surface. The coefficients of the series are adjusted so as to satisfy the flow boundary 
condition at a finite number of control points. This is a collocation method; numerous other 
methods have been used. Reference to the books by Ashley and Landahl (ref. 27) and Ward 
(ref. 30) will provide an introduction to many of the alternate approaches. 

The approach used by Woodward (ref. 35) works well with the method that includes 
the empirical corrections required by stability and control applications. It relates a localized 
cause to a localized effect so that localized empirical corrections may be introduced. 
Aerodynamic influence coefficients are generated and admit of the following aerodynamic 
equation 

{F A } = q([A]{j&}+fA T ]{^ T }) (138) 


where the elements F A ^ are the aerodynamic panel pressure forces. The elements are the 
incidence angles of the panels and i£y. are the panel incidences due to thickness. The 
matrices [A] and [Ay] are matrices of the corresponding aerodynamic influence 
coefficients. 

Up to this point the discussion of the aerodynamic problem has ignored the effects of 
unsteady flow; the result represented by equation (138) is valid only for steady or 
quasi-steady flow. As shown in app. B, however, an approach suggested in Miles’ monograph 
(ref. 28) may be used to include the effects of unsteady flow in a manner appropriate to the 
aerodynamics of stability and control. This approach leads to a direct extension of equation 
(138) when the time dependence of the flow incidence angles is slowly varying. As shown in 
app. B, the time variation must be such that its reduced frequency is in the region of 
reduced frequencies of order of magnitude zero. This is a good approximation for 
consideration of the stability characteristics of large, flexible airplanes even though it is not 
satisfactory f or flutter problems. 



The consequence of mechanizing this approach following the method by Woodward 
will lead to an aerodynamic equation of the form 

{F a } =q ([A] {$ } + [6A]{i}+ f A t 1 {^ T }) (139) 

fc 

The elements of [ 6A] will relate the panel pressure forces to the rates of panel incidence 
change. 

By evaluating the boundary condition given by equations (132) through (136) at the 
aerodynamic control points and substituting it into equation (139), a complete formulation 
of the aerodynamic problem of stability and control is obtained. In the reference motion 
case, time dependence drops from consideration, and the pressure force due to thickness 
enters the problem only through the elastic deformations. However, an additional 
consideration is necessary in the case of disturbed motion. This involves the coordinate axis 
systems used to describe the problem. 


The aerodynamic equation, equation (138) or (139), is derived in terms of the fluid 
axis system. It may be recalled that the freestream velocity, U, and the Mach number, M, are 
perturbed. The fluid axis system in the reference motion differs from that in the disturbed 
motion. The aerodynamic equation is perturbed, and this may be expressed by writing 
equation (139) as 


{ I 'A 1 } + < f A }'<Vq p > 


I6A){i} + ^[ a tJ + ["air] 


(140) 


By discarding products of perturbation quantities and using the form of equation (139) for 
the ’.eference motion, it follows that 



+ q x [A 1 ] {j^pl + q! f6A]{^ p } 


(141) 
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This expresses the perturbation to the aerodynamic forces in terms of the perturbation 
motion variables. The first-order approximation, qp =» 2qj u/Uj, has been used. Also, it 
may be noted that the first term contains products of perturbation quantities since the 
reference motion incidence angles must be of perturbation order. The justification of this is 
somewhat loose. It falls on the requirements of stability and control. The u perturbation 
variable must be treated independently of the other variables. The first term, which contains 
u solely, must be considered separately in the order of magnitude consideration. 

This last consideration points up a problem that requires a great deal of further 
investigation. Equation (141) has been written nearly a priori, but the basis of the equation 
is far from self evident. A very careful and fundamental order-of-magnitude analysis is 
required in order to obtain a consistent, valid formulation of the influence coefficients in 
the equation. This is illustrated by the development of the results, which appear in app. B. A 
derivation of equation (141) is carried out for a thin body undergoing longitudinal 
perturbations. The basis for a sour'* extension of that analysis to thin wing, slender body 
combinations undergoing longitudinal disturbances also appears to exist (refs. 29 and 35). 
However, the lateral-directional problem for wing-body combinations has not been 
developed. An intuitive approach to developing a computer program for solving this 
lateral-directional problem might be instituted. However, without an order-of-magnitude 
analysis there would be no way to assess the limitations of such a program. A failure of the 
program to accurately predict lateral-directional derivatives might defy further intuitive 
judgment and a large engineering effort might be needlessly lost in numerical 
experimentation. 


In the preceding equations, (133) through (141), elastic deformations must be 
expressed in terms of the free vibration mode shapes. The perturbation to the elastic 
displacement vector is given by 

' ( d p} - (M2) 

where the generalized perturbation displacements Up. are functions of time alone. The 
displacement rate dp in equation (133), therefore, can be written in matrix form as 


{%}= I ♦H»p} 


The triple scalar product (0Ep ;; n l ) *i s may be written as 




( 143 ) 


(144) 
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•where 4 p and 0p are the components of elastic rotation about the z„ and y c axes, 

L 'p ”p ^ 3 

respectively. Two additional modal matrices are introduced such that 

|^E p j = ( V ^ j^J^ 0 © 1 ^p* (145) 

These represent the elastic slope deflections required for equation (144), but it must be 
noted that the free vibration mode shapes do not describe the elastic deformation 
completely when the mode shapes are generated in the manner described in app. A. The 
elastic deflections arising from perturbation changes to the orientation of the gravity force 
and rigid-body inertial forces must be introduced separately. This is described in the 
discussion leading to equation (6.86b) in app. A as well as in the disucssion given in par. 
6.3.4 of app. A. 

Denote the elastic rotations arising from inertial relief and gravity perturbation forces 
as 

{%} “ d Kl 


The elastic rotations may be found from specialized forms of the flexibility matrix denoted 
by [C 0 ] and [Cq] . Then, in accordance with equation (6. 1 65) of app. A 

j*E p | = - [s]M I* ) (iMy + lM l' < V p > +[M 2'<r,;>) 

and (146) 

| 6 E | = -[' 6 ]WI'|(hIV * :M ! 1 (V p > + IM 2 l{r' } 

i P ' 


Further, letting the matrices of the components of the normal vectors nj. be defined as 

cy ■ 7 S - fy . cy • i - cy . Pd • \ ■ c-y 

the expression for the perturbed flow incidences, equation (143), may be written in matrix 
form as 

<V by ^ f u p>) < 147! 
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( 147 ) 

By introducing appropriate matrix definitions, equation (147) may be written as 

(Jt p }=[G 1 ] {y+[G 2 ){y + tG 3 l } + W{ e E 1(148) 

p p 

where 

{V } T = |u, v, w, p, q, rj 

r 

as in Sec. 5. 

The perturbation aerodynamic equation, equation (141), may be written in an 
abbreviated form by letting 

(Ajl - [ 2 (A 1 ) + Mj [||]] {»j} 4 [2 [A T J 4 Mj {if T } (!«) 

so that on introducing equation (148) 

{f A } = {G i ] {y + V 6A1{ V 

+ q x [Aj] l G 2 ] {u p } + q x [fA^ [Gg] + [6A] (G 2 l] { u p } 

+ q 1 [<5A] (G 3 1 {ii p } + q x (A^ £n y ] }• + {nj { 0 ^ (150) 
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where the effects of time rates of change of and {0'Ep}h ave been con-idered 

negligible. The aerodynamic matrices introduced in the equations may now be recognized as 

* 

fc 


[Aj] =q 1 [{Aj}-^, {0}, {OHOHO} {0}]+^ [Aj [ 

(lbla) 

(A 2 1 =q 4 [l6Al + [A 1 l ^-[n y j [C*] -r ' [C 0 ]^mJ (0] T J 

(151b) 

(Agl A^] f G 2 ] + [A^] + L U Z J ( C Q C f^] T 

(MjjJ (151c) 

[[AJ [G 3 ] + (6A![G 2 ]J 

(151d) 

(A g ] =q 1 [6Al [Gg] 

(151e) 


Hence, with these definitions equation (150) becomes identical to equation (126) of Sec. 6, 
i.e., 


{f A } =[A iHV p }+(A 2 ] {V p } + [A 3 ](0] {u p } + [A 4 ][0] {u p } 


+ [A g ] 10] { U p } 


(152) 


7.3 Semi-Empirical Methods 

The semi-empirical or “handbook” methods have evolved from numerous wind tunnel 
tests and theoretical analyses of various aircraft configurations to evaluate the effects of 
geometry on the stability derivatives. These techniques, as do both the purely experimental 
and purely theoretical techniques, seek to solve fie exact flow equations of ref. 33. In 
addition, the handbook techniques seek to improve the theoretical solutions based on 
inviscid fluid dynamic results by including experimentally measured forces and moments 
due to leading edge suction, boundary layer, and other nonlinear effects not included in the 
inviscid approximation of the aerodynamics. The inviscid flow theory is then modified with 
the experimental results and plotted as functions of aerodynamic and geometric parameters 
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(ref. 6). An engineer interested in calculating stability derivatives for a particular 
configuration must seek the appropriate tables and select the separate effects due to wing, 
body, tail, canard, etc. Appropriate interference factors between wing-body, wing-tail, etc., 
are also presented in the tables to account for the nonlinearity of the flow field. Appendix 15 
indicates which of the stability derivatives can be calculated using the USAF Stability and 
Control Handbook. 


The procedure used to calculate the stability derivatives by the handbook technique is 
best illustrated by an example. Consider the calculation of Ci : 

n _ Of 


C L = 


a ( C L a ) e S' [ IC N + ^VfB) + IC B(W)] + (153) 

WING-BODY TERM / \"L„ „ V' / 9€ \j_ S" K 

\ % ) p/(B) K B(W)J 


HORIZONTAL TAIL TERM 


where (Cj )' 


is the lift curve slope for the wing of appropriate aspect ratio, 
taper ratio, section profile (camber, thickness), twist distribution, 
etc. 


Kj^' is the ratio of aircraft nose lift to aircraft wing lift. 


K B(W) 


is the effect of lift carryover on the wing due to the body. 

is the effect of lift carryover on the body due to the wing. 

is the lift curve slope of a tail with the appropriate geometric 
parameters. 


Kw(b) is the effect of lift carryover . n the tail due to the body. 
Kg(W) . is the effect of lift carryover on the body due to the tail. 


d€ 

da 


is the change in downwash on the tail due to the change in 
angle of attack of the airplane. 


q M is the ratio of dynamic pressure at the tail to freestream dynamic 

q« pressure. 

s; is the ratio of wing area to the “effective area” of the tail. 

Se" 
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Se M 

S" 


is the ratio of effective tail area to the actual tail area. 



Equation (153) represents summed knowledge based on the personal experience of the 
authors of the handbook. It should be emphasized that the approximation of C^q, in the 
equation is not unique because of the “experience factor”'dnvolved. Therefore, there is no 
guarantee that the approximation is correct for all configurations. The handbook technique 
is the only semi-empirical technique available at this time that can be used to estimate 
nonlinear effects in the viscous fluid dynamic equations, Accuracy by this technique is only 
a function of the ingenuity, insight, and experience of the engineer performing the analysis. 

One other disadvantage of the handbook technique is that it can only be used to 
evaluate th • rigid and equivalent elastic mathematical models of an elastic aircraft. It cannot 
be used on the completely elastic and residual elastic models because there is no known way 
to develop consistent structural derivatives dependent only upon external geometry 
parameters. 

Handbook data for a study SST configuration and the 707-320B are presented in app. 
B. 


7.4 Experimental Methods 

The two principal experimental methods used to extract the stability derivatives are: 
(1) wind tunnel tests of a scaled model of the proposed aircraft and (2) flight tests of a 
similar aircraft or the actual aircraft. Both of the experimental methods are potentially more 
accurate than either the theoretical or semi-empirical methods, but neither offers a useful 
and convenient way of optimizing a configuration for good stability and control 
characteristics. 

The main disadvantage of both experimental methods is the high cost of each test 
point and the difficulty of correcting the collected data for errors due to experimental 
procedure. In addition, the only derivatives generated by a flight test technique are for 
completely elastic airplanes, and the extraction of rigid airplane values for comparison with 
the results of other methods is difficult, if not impossible. 

7.4.1 Wind tunnel tests. — Wind tunnel tests are a reliable way of measuring rigid 
aircraft static stability derivatives. In addition, several experimental techniques to measure 
the a , q, and p dynamic stability derivatives are currently under development. 

Recently, a new technique of generating the static stability derivatives for an 
equivalent elastic aircraft was developed and is outlined in app. B. A semi-elastic model is 
built with structurally and aerodynamically scaled wings and tail. Force and moment data, if 
properly adjusted for scale effects and tunnel wall effects, can then be reduced to give the 



equivalent elastic stability derivatives. These derivatives are then used in the dynamic 
equations describing the airplane’s motion through space. Inertial effects must be included 
in the analysis to account for load factor changes. •, 

i 

In some cases the equivalent elastic stability derivatives can be calculated from rigid 
wind tunnel data. In this method the lifting surface computer program is used to analyze the 
rigid and elastic aircraft First, rigid lift curve slopes are calculated; then, using the 
appropriate structural influence coefficients, the elastic lift curve slopes are calculated. The 
ratio of Cg a g/Cg a R = Lg/L^ is formed and applied to all subsequent rigid wind tunnel 
d.«.a to calculate the effect of elasticity. For some configurations, such as the SST of this 
study, the effect of elasticity on rigid wind tunnel values can be estimated by an elastic 
increment (see par. 7.5). 

Several difficulties arise in the use of wind tunnel data. Because the model parameters 
used to evaluate stability derivatives have not yet been outlined, the engineer must 
determine the correct trip-strip size and location for each stability derivative to duplicate 
the proper nonlinear viscous effects. The effects of Reynolds number and “rigid model 
elasticity” are usually unknown. Also, some wind tunnel tests may be plagued with high 
tunnel turbulence or other forms of wind-tun..el-induced effects in the desired force and 
moment data. 

Data for a study SST configuration and the 707-320B are presented in app. B. 


7.4.2 Flight test techniques. — The evaluation of stability derivatives from flight test 
data should, by definition, give the most accurate elastic airplane values. The engineer can 
use such data in much the same way that he uses the handbook technique, provided 
individual component contributions to each stability derivative are known or can be 
extracted from the data measured in flight testing, 

The primary problem associated with stability derivative evaluation from flight test 
dau is the difficulty of performing maneuvers holding all motion variables, except one, 
equal to zero. Obviously, even if it were possible to fly only at an angle-of-attack variation, 
the aircraft would elastically deform at some frequency determined by the structural 
properties, and the force and moment data measured at the center of gravity would contain 
both a and structural motion (dp) contributions. 

In addition, there are a limited number of maneuvers available from which to extract 
the stability derivatives. The usual methods for the longitudinal stability derivatives are 
' control column steps and pulses, windup turns, and thrust steps and pulses. Lateral 
derivatives may be extracted from maneuvers arising from steady sideslip and from wheel 
and rudder-pedal steps and pulses. The resulting forces and moments measured at the 
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aircraft center of gravity by accelerometers re.iect not only aerodynamic effects bi ,f also the 
effects of the pilot's transfer function (la^s, response time, insensitivity to sm.-ll but 
measurable errors o r trim, etc.) and the control system elasticity (cable stretch, internal 
friction, and nonlinear effects due to the wear of hvdraufic systems and control surface 
linkages). Thus, the stability derivatives extracted from such flight test data also reflect 
these errors. 

Data for the 707-320B are presented in app. B. 


7.5 Comparison of Methods 

Figures 6 through 17 present comparisons of the theoretic ’ ^emi-empirical, and 
experimental evaluations of Cj_^, ( mQi , C m ^, Cj^, and Cj^ stability derivatives. The 
experimentally determined values (wind tunnel and flight test) are to be taken as the basis 
for comparing accuracy. A comparison of the estimates c f all stability derivatives listed in 
table 3 is contained in app. B. Discussion of the approximation?, used in the estimates also 
appears in app. B. 

The stability derivative evaluations presented here were chosen because a complete 
comparison of theoretical and semi-empirical methods cou !i be made against both wind 
tunnel and flight test data. As shown by tables 6, 7, and 8, estimates are not possible fc 
both techniques for all stability derivatives of interest. Further, testing techniques do not 
exist or have not been used to evaluate all of these stability derivatives for the 707-3 20B and 
the SST configuration. Thus, a basis for comparison does not exist for all of the desired 
calculations. 

Wind tunnel evaluation of and C mQ! from both rigid and elastic wind tunnel 
models is shown in figs. 16 and 17. These data are shown in comparison with estimates 
obtained from lifting surface theory. The results are functions, of both dynamic pressure and 
Mach number, as noted on the figures. Consequently, the interconnecting lines between 
values for the elastic model are for visual purposes only; .nterpolation for intermediate 
values is not possible. 

Primary impc tance has been placed on the calculat'on of the static and quasi-steady 
longitudinal derivatives, since for these derivatives all three techniques are applicable. Four 
of the six pri nary longitudinal derivatives are presented in figs. 6 through 13 to compare fhe 
techniques. For both the 707 an 1 SST, the lifting surface values of Cl c r (figs. 6 and 7) are 
closer than the handbook values to the wir... *unnel values. The value of C mQ; ^ for the 7C7 
(fig. 8) calculated by lifting surface theory is closer to the wind funnel value than arc vain*, s 
calculated by either lifting line or handbi. 'K techniques. In the case rf the SST, C ni( r is 
also calculated more accurately by the lifting surface technique than by the hanubook 



technique. No experimental data were available to evaluate the accuracy of C m r (figs. 10 
and 11) for the 707 and SST. The value of Cjn u ^ (figs. 12 and 13) calculated by lifting 
surface theory compares very well with wind tunnel data .for both the 707 and SST; the 
handbook technique will not handle th : s derivative. 

The only significant lateral-directional derivative listed in table 3 that can be calculated 
by both lifting surface and handbook techniques is Q^. Figures 14 and 15 show the 
comparison for Qp^; unfortunately no rigid experimental values were available. 

The comparison of the techniques used to calculate elastic stability derivatives is 
hampered by an inconsistent set of data. Of the five primary derivatives shown in figs. 6 
through 15, Qpg for the 707 (fig. 14) provides the best test of accuracy. Both the 
handbook and lifting surface techniques give the same value of Qpg at 3050 meters 
altitude, tut at 10 675 meters the lifting surface technique is more accurate. For this reason, 
and because lifting surface theory is more accurate for rigid derivatives, it is concluded that 
the remaining equivalent elastic derivatives will also be calculated more accurately by the 
lifting surface technique than by the handbook technique. 

As additional justification, figs. 16 and 17 show that the incremental error between 
wind tunnel and lifting surface values of and C m ^ does not increase appreciably 
between the rigid and the equivalent elastic cases. In fact, figs. 16 and 17 indicate that an 
empirical correction to the rigid wind tunnel data may result in a more accurate elastic value 
of the stability derivative. For the study SST at 72° sweep, these corrections are of the form 
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Change in lift coefficient with angle of attack, Ci , deg 




Change in Lift Coefficient with Angle of Attack, Cj_ , Deg 








Change in pitching moment coefficient with angle of attack, C m , deg’ 




Change in pitching moment with angle of 
attack, C^.deg' 1 
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FIGURE 10. -COMPARISON OF ESTIMA TION METHODS FOR C m - 707-320B 
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Change in pitching moment with pitch rate,C m q, rad^ 
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FIGURE 12. -COMPARISON OF ESTIMA TION METHODS FOR C - 707-320B 

u 


131 


Change in pitching moment with i rward velocity 
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FIGURE 14. -COMPARISON OF ESTIMATION METHODS FOR C* --707-320B 
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£.0 METHODS FOR CALCULATING STABILITY 
AND RESPONSE CHARACTERISTICS 

% 

8.1 General Considerations 

The determination of static stability and control characteristics may be considered as 
an interne diate step in the solution of the total problem of dynamic stability and control. 
Taking this intermediate step has a number of advantages. 

1. Steady-state trim, balancing, and control information are provided. 

2. Comparison with static stability criteria can be made. 

3. Through examining the static picture it is possible to predict some dynamic 
characteristics. 

However, the analyst is limited in trying to predict the dynamic motion from static 
stability considerations for the following reasons. 

1. Static longitudinal stability is, in general, a prerequisite for dynamic longitudinal 
stability. There ate, however, instances where static stability is not required for 
dynamic stability. 

2. Static stability is usually, but not necessarily, required for good handling qualities. 

3. Structural dynamic motions cannot be investigated. Static stability analyses are only 
concerned with the rigid and equivalent elastic airplanes. 

Static stability characteristics are found by an analysis of the stability derivatives (using 
calculation methods discussed in Sec. 7), and by solving for the neutral point, maneuver 
point, etc. through fairly elementary expressions containing the stability derivatives. 

Dynamic stability characteristics are generally determined by analyzing the roots of the 
characteristic equation when the perturbations of the airplane motion from the steady state 
are small. Approximate formulas for the frequency and damping of the motion are at times 
applicable for small perturbation motion. When nonlinear effects are significant, time 
history solutions of the arbitrary or large pertuii ration equations are used to study the 
stability characteristics. 

The three airplane models treated in this study, viz., rigid, equivalent elastic, and 
completely elastic airplanes, were derived from consideration of the degree of airplane 
flexibility. The completely elastic airplane has been analyzed with various elastic degrees of 
freedom, and the effects of residual flexibility have been studied. The various structural 
mathematical models hav* been reviewed in Sec. 5. Static and dynamic stability character- 
istic.. were determined by using stability derivatives from several sources: handbook 
formulas, aerodynamic lifting surface theory, amt wind tunnel data. 



General results obtained from the static and dynamic analyses are presented in tables 9 
and 10. Table 9 summarizes the accuracies obtained from lifting surface theory and hand- 
book techniques in predicting the rigid airplane static stability and control characteristics 
when compared with wind tunnel predictions. A limited amount of substantiation with 
flight test data was obtained. Also shown in table 9 is the relative effect of elasticity on the 
various static characteristics. Although most of the substantiation of methods was for the 
rigid airplane, one could expect to obtain similar accuracies for the equivalent elastic air- 
plane. A method based upon ratioing the elastic and rigid stability derivatives that can be 
expected to improve the prediction of static stability characteristics is discussed in par. 8.2. 
The poor accuracy obtained for some configurations and characteristics is almost entirely 
due to poor prediction of the stability derivative . Appendix B discusses the calculation 
of this derivative and the expected improvements to the lifting surface theory mechanisa- 
tion program to improve the prediction of C^. 

Table 10 presents the same general results for the dynamic characteristics of the 
airplane studied.* In addition, general results on the usefulness of approximate formulas for 
predicting frequency and damping are summarized, and the number of elastic modes needed 
for an accurate dynamic elastic analysis is given. Appendix C should be consulted for details 
of the derivatives that were used in an individual method and for a complete discussion of 
all the dynamic stability results. A limited discussion appears in par. 8.3.5 for some of the 
more important results. 


8.2 Static Stability Characteristics 

The static stability characteristics of an airplane are strongly dependent on the indi- 
vidual stability derivatives and on how the stability derivatives combine. The effects of the 
individual stability derivatives can be judged by comparing their signs with the static stabil- 
ity criteria, as discussed in Sec. 6, and by noting the magnitude of the derivative. The sign of 
ttie derivative simply indicates whether the airplane is stable, unstable, or neutrally stable 
with respect to a certam motion variable (requirements for stable motion are summarized in 
table 11). The magnitude of the derivative gives an indication of the degree of stability or 
instability. 

The derivatives in combination can also be used to evaluate airplane stability character- 
istics. The static stability and control characteristics usually investigated are: 

1. Elevator and stabilizer trim angles; 

2. Stick-speed stability; 

3. Elevator and stabilizer angles per g; 

4. Neutral point; 

5. Maneuver point. 

Expressions used for calculating the above quantities are given in equations (154) 
through (161). 



TABLE 9.- STATIC STABILITY CAL CULA TIONS-GENERAL RESULTS 



Relative accuracy of calculation method 2 

Relative 
effect of 
elasticity' 1 

Computer lifting surface, 
rigid and equivalent elastic 
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a. Reflects almost entirely ability to calculate derivative C mQi and resulting effect on 
characteristic. 

G (good)— method compares favorably with wind tunnel predictions (exception allowed.) 
F (fair)— less favorable correlation with predictions. 

P (poor)-method does not compare favorably with predictions. 

b. L (large)-elasticity considered a significant effect. 

M (moderate)— elasticity considered moderately important; not quite as significant as 
differences due to stability derivative calculation methods. 

S (small)-elasticity considered a minor change to stability characteristic; changes due 
to stability derivative calculation methods usually much more important. 

c. No data available. 

d. Correlation with flight test, but based on a very limited amount of data. 
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TABLE 10 - DYNAMIC STABILITY CALCULATIONS-GENERAL RESULTS 
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8.2.1 Elevator and stabilizer trim angles. -For rigid and equivalent elastic airplanes, the 
horizontal stabilizer angle to trim is given by 







Similarly, the elvator angle to trim is given by 
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1 






(154) 


(155) 


8.2.2 Stick-speed lability. -The elastic airplane stick-speed stability for stabilizer trim 
is given by the following expression: 
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For the rigid airplane all variations with dynamic pressure vanish, i.e., 3C mo /3q = 
8Ci^/3q = 0, etc. Therefore, for the rigid airplane equation (156) becomes 
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As shown in table 1 1 , a stable gradient of longitudinal control displacement versus 
speed is defined as one for which dijj/dVj n= .j >0- Stick-speed stability is usually referred to 
as a handling-qualities parameter. The equations for elevator speed stability are, of course, 
identical with equations (156) and 157), except that 6g replaces ijj in all terms. 

8.2.3 Elevator angle per g, — An expression for elevator angle per g in a normal pullup 
can be written as 
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Equation (158) holds for both rigid and equivalent elastic airplanes. A stable gradient 
of elevator displacement versus load factor is defined as one that satisfies 
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(159) 


The term containing the stability derivative Cg is usually insignificant. 

q 

8.2.4 Neutral point. — The neutral point is calculated from the expression 
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where h is the center-of-gravity position. The requirement on the neutral point is that it 
should be aft of the aft center-of-gravity limit. Equation (160) is good only for small angles 
of attack in the linear range. „ 


8.2.5 Maneuver point.-An expression to calculate the maneuver point is 




( 161 ) 


The effect of the derivatives Cj^ on the maneuver point is negligible, as it was for the 
elevator-angle-per-g characteristic. The requirement on the maneuver point is that it should 
be aft of the aft center-of-gravity limit. 


Complete derivations of the above equations for the static stability characteristics can 
be found in app. C. 


Five methods were used in tb ! : study to determine the longitudinal derivatives. 



Applicability 

Rigid 

Airplane 

Equivalent 

Elastic 

Airplane 

1. Computer, using lifting surface theory 

yes 

yes 

(aerodynamic influence coefficient 



method) 



2. Handbook 

yes 

no 

3. Handbook combined with computer 

no 

yes 

lifting surface theory 



4. Wind tunnel 

yes 

yes 



(limited) 

5. Flight test 

no 

yes 

— — — 


(limited) 


The characteristics of the rigid airplane were compared using the previously listed 
expressions for static stability and control. 


Elastic stability characteristics were determined using equivalent elastic stability deriva- 
tives of the Formulation II type (app. B). These derivatives were generated by methods 1. 
and 3. above. 
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Some flight test data were available to correlate theoretical predictions of longitudinal 
static stability. The data were obtained as part of the certification requirements for longi- 
tudinal stability and control for the model 707-320B airplane; Results for both maneuvering 
stability, dd^/dn and stick-speed stability, dfig/dV, were obtained from these tests. The 
elevator-angle-per-g information was reduced from windup turn maneuvers. 


Figure 18 illustrates some typical static margin results for the 707-320B and the 
42® -wing-sweep SST configuration. Variations in the results due to both stability derivative 
calculation technique and elasticity can be noted. Appendix C contains detailed results for 
all characteristics in graphic form. 


Lifting surface theory (aerodynamic influence coefficient method) was found from the 
static stability analysis to give better predictions than handbook methods for some cases 
(fig. 18 and table 9). It gives direct, acceptable results for some stability characteristics for 
equivalent elastic and rigid airplanes. However, if wind tunnel data are available a more 
accurate way of- predicting elastic effects would be to compute an elastio-to-rigid ratio or 
increment referenced to the wind tunnel value. This was substantiated in Sec. 7. For 
example, 



EqEl 



m, 


a. 


Elastic 


m. 


“Rigid 


Computed 


8.3 An Analysis of Dynamic Stability by 
Characteristic Equation Rooting 

8.3.1 Applicability of characteristic equation methods. — Airplane equations of motion 
can be reduced to a set of linear, second-order differential equations v/ith constant coeffi- 
cients when dynamic behavior can be approximated by assuming that motion perturbations 
relative to the steady state are small (see Sec. 5). These equations are called small perturba- 
tion equations of motion and are amenable to generating characteristic equations whose 
roots can be examined to determine motion characteristics. 

The dynamic stability criteria for all characteristic equation roots are given in par. 
6.2.1. These criteria can be satisfied by inspection, i.e., by checking the sign of the real parts 
of the roots or their absence. However, these yes/ no-type answers relate very little informa- 
tion about airplane motion characteristics. Some of the parameters that can be deduced 
from the roots and are more physically oriented are discussed in par. 8.3.2. 



Static margin, h n • h, X/Cr Static margin, h n • h, X/c 



FIGURE 18.~ 707'320B and 42° SS 1 SI «/G/D KS. EQUIVALENT ELASTIC 
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8.3.2 Rigid and equivalent elastic airplanes.— The small perturbation longitudinal and 
lateral-directional equations of motion are given as equations (17) and (18). It is possible to 
take the Laplace transformation of the equations of motion and solve for the roots of the 
resulting characteristic equations. However, the program that generated the roots and asso- 
ciated data in this study used a different technique, which will now be described briefly. 

The lateral-directional characteristic equation, with the indicated rigid/equivalent 
elastic definitions, evolves from equations (17) and (18) as indicated in the equation in table 
, 12. This equation is a result of the requirement that (for 6^ = <5^= 0) the equation 

GM D >] { x } = {°} (162a) 

has nontrivial solutions, i.e., 



The longitudinal characteristic equation in determinant form evolves as indicated in the 
equation in table 13. Inertia relief is handled implicitly as indicated in table 13 (called 
Formulation II in app. B). The effects due to 0 cannot be treated explicitly and are handled 
as summarized in table 13. 


The expanded form of the characteristic equations of tables 12 and 13 for the small 
perturbation equations of motion written in the parameter X is 


AX 4 + BX 3 + CX 2 + DX + E = 0 


(163) 


where the coefficients A through E are determined by the case (longitudinal or lateral- 
directional). Equation (163) is obtained by assuming the solutions 
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(0=V xtA * 

L v xt/t * 

( f 0 e Xt/t * 


(164) 


(165) 


and substituting these into the equations in the form of equation (162a). After carrying out 
the differentiation, e^*/ 1 * ''an be eliminated, leaving 



TABLE 12.- LATERAL-DIRECTIONAL CHARACTERISTIC EQUATION IN DETERMINANT FORM 


(2 M - Cy£ ) D - C y ^ -(C yp D + C Lj ) (2 m- C yr ) 


-< C | j g +C ^ D ) 

^ C n„ + C n;, D) 


-C| p D 


-(iED + C.) 


= 0 


-(ipD^ + C- D) i p D - C n 
“fl “p t n p c n r 

where: M = M/pS w l; M = W/g; i=b/2 

D = t* d/dt; t*= l/V. 


Note: For the equivalent elastic 

airplane, Cy£ p Cy^, Cf jr , and C n ” 

- ‘ are not accounted for. 


Rigid airplane 

© 

Equivalent plastic airplane 
© & © 

All derivatives are conventional 
rigid derivatives. 

Derivatives include elastic effects due 
to aerodynamic loading (Formulation I) 
and would have E subscripts. 

*A * 

«c * WpV 3 

*E = WPS w <3 


*A = ^xx' ^ip{ S w b)/pS w i3 
iC = (^zz * Cnfj Q S\yb)/pS w 4^ 

iV _ (*xz + 1/zC lfi ^ s wb + '/iCnpj q S w b) 

E PSw 13 

W = Actual Gross Weight - W j 

W = g(M - Cyyj q S w ); M=W,/g 


Assumptions required: 

D1 and 2, D5, D7, D1 1,G1-6,G10, SI and2,S5 and Rj = 0j * = 0 




TABLE 13.- LONGITUDINAL CHARACTERISTIC EQUATION IN DETERMINANT FORM 


2pD + Cj) u + 2Cj)j 

c D a - c L, 

CLi 1 

C L U + 2C r.i 

(2p + C L JD + (Cl^+ Cj) j) 

-(2m-C l )D 


"fl-m igP + 

ijjD 2 - C D 
m q J 


where: P - M/PS W 4; M = W/g; I = c/2 


D=t* d/dt; t*=l/V cl 


Note: For the equivalent elastic 
airplane, CL^and all uj 
are not accounted for. 


Rigid airplane 

© 

Equivalent elastic airplane 

(S) & (S7) 

All derivatives are conventional 
rigid derivatives. 

Derivatives include elastic effects 
due to both aerodynamic and 
inertial loading (Formulation II) 
and would have E subscripts. 

ifi “ ^yy/^w^ 


= (lyy = ^mQj 5 S w c)/pS w i^ 

W = Actual airplane weight 

W = Actual airplane weight 


Assumptions required: 


Dl-2,5,7,8, 11.G1-6, 10, Sl-2,5 and Ri =0! = = C Tmi +C mi =0 




IA («]{*„} ={0} 


(166) 


where [A(A)J is identical with [A(D)J for A = D. The solutions Aj/t*, i = 1, 2, 3, 4 of the 
fourth-order polynominal equation (163) are the “roots of the characteristic equation.” 
Since a rather complete discussion of the coefficients A through E and the occurrence and 
significance of various combinations of real and complex roots may be found in refs. 4 and 
36, it is not repeated here. 

The systems analysis analogies used for airplane characteristic equations lead to the 
conventional “mode” definitions which follow (in equation (163) let A j/t* = S and divide by 
A). 


First, for longitudinal equations (for two complex pairs of roots): 

(S 2 + 2f p u 2 p ) (S 2 + 2 f sp u sp S + u 2 sp )-0 


where: p ~ phugoid mode 

sp ~ short period mode 


the roots are 


-£p(sp) w np( sp) 

Wn p(sp)‘\/ 1 * ^p(sp) = w p(sp) 


is the real part of the 
phugoid (short period) 
root pair 

is the imaginary part of 
the phugoid (short period) 
root pair 


s l,2 = ^p =*= J^p 
S 3,4 =0 sp*i w sp 

= a P(sp) 


^p(sp) * s the phugoid (short period) damping ratio 

Wn p(sp) * s ph u S°*d (short period) undamped 
natural frequency 
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Second, for lateral^diiectional equations (for one complex pair and two real roots): 

(S 2 + 2f D w nD S + W 2 D )( S + T | : )(S+^) =0 

where: D ~ Dutch roll mode 

r ~ro’ling convergence root 
s ~ spiral root 

-U.L 

T r T s 

the roots are Sj 2 = <t d ± J w D 

S 3,4 = _ 4^»*i^ 55 °3> °4 

- ?£) a) n jj = 0j) is the real part of the Dutch 
roll root pair 

«n D Vr^= W D is the imaginary part of the 
Dutch roll root pair 

is the Dutch roll damping ratio 

td n p is the Dutch roll undamped natural 
frequency 

The relationship between the damping factor (£), undamped natural frequency (w n )> 
damping (a) frequency, and damped frequency (w) is illustrated in the diagram below. 



♦Sometimes the term (S+4*) is written (S-=r). 

l s h 
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o - -£w. 


1/2 


w = Wn (1 - t ) 


£ = - cos 6 


Assume that there are two complex pairs of roots that result from the solution of the 
rigid or equivalent elastic longitudinal characteristic equation. The real-time solution fora 
will ther ! > ve the form 


a = e p sin (w t 
P 


- P l> 


CTgt 

e 8p sin fe sp t - P 2 ) 


( 167 ) 


where Pjand P 2 are constants determined by initial conditions. 
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One pair of roots, the highest frequency pair, determines the short period mode. 
The damped frequency is given by 

w S p~ radians /second ’ 


The period is given by 


P„ =TT ~ seconds /cycle 
S P w sp 


(168) 


Time to damp to half amplitude for cr S p <0 is given by 

T i/2 = ^ p ln{ i ) ~ seconds 

For cr S p'- 0, time to double amplitude is given by 

T 0 = ~ In (2) ~ seconds 
* °sp 

Cycles to damp to half amplitude is given by 


(169) 


(170) 


r T l/2 M sp T l/2 
1/2 P sp ‘ 2" 


cycles 


(171) 


Cycles to double amplitude is given by 


C = — 
2 P 


sp 


Cl 

2tt 


> cycles 


(172) 
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These parameters apply also to the longitudinal phugoid mode, the lateral-directional 
Dutch roll mode and, in some cases, a lateral phugoid mode when the lateral-directional 
characteristic equation yields two pairs of complex roots. . 

8.3.3 Completely elastic airplane. - For the completely elastic airplane a less restricted 
approach is required than that used for the rigid and equivalent elastic cases. This involves a 
mathematical model that will account for the structural dynamic motions of the airplane as 
well as the flight path and rotational motions. It is accomplished using a model with the 
option of including or excluding residual flexibility, thus giving a system with an arbitrary 
number of variables. The airplane has, then, the usual six degrees of freedom plus an 
arbitrary number of degrees of freedom that involves the structural dynamics. 

Subject to the type of problem to be solved and the degree of accuracy required, the 
engineer has a choice as to the number of variables (degrees of freedom) to include in any 
one analysis. This is a considerable departure from the philosophy of the well-defined, six 
rigid-body degrees of freedom associated with the rigid and equivalent elastic mathematical 
models previously discussed. 

The equations of motion that represent both rigid-body and internal motion have been 
developed as equation (95) in Sec. 5 and are repeated here 

[a]{q}+ [Bl {q} + tC]{q} = {0} (173) 


Expressions for the coefficients [a] and [B] have been presented in Sec. 5 for a residual- 
flexibility formulation and for the case where all elastic degrees of freedom participate 
dynamically. For the latter case, i.e., a completely elastic airplane, the coefficient [a] 
represents the generalized mass, [B] includes the aerodynamic damping, and [C] the gener- 
alized stiffness and generalized displacement dependent aerodynamic coefficients. 

Taking the Laplace transformation of equation (173) yields 

[fa) S 2 + [B] S + ICl]{q(S)}={0} (174) 

NOTE: From the definition of perturbation variables {q}, 

{q (t = 0) } = {q (t = 0)} = {o} 
which will have nontrivial solutions {q} only if 

|[fa]S 2 + (B]S + (C ]]| = 0 (175) 
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The characteristic equation (175) then yields a determinant with elements ajjS^ + BjjS 
+ Cjj. When equation (175) is expanded, it yields a polynominal of degree < 2n, where n is 
the order of the determinant. The roots of this polynominal are the roots of the character- 
istic equation (175). These roots are obtained using an eigenvalue approach. The theory is 
described in refs. 37 and 38. 

8.3.4 Approximate solutions.— Approximate solutions for rigid airplane frequency 
and damping characteristics have long been in existence. It is assumed that the appropriate 
solutions also apply to the equivalent elastic mathematical model because of its similarity to 
the rigid model. 

An extensive discussion of approximate characteristics, transfer functions, etc., can be 
found in ref. 36. For a two-degree-of-freedom ( aand 9 ), longitudinal short period mode 
approximation, the expressions for frequency and damping are: 
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Data in both refs. 4 and 36 show these to be accurate expressions for certain rigid 
airplanes when compared with the exact quartic solution of equation (166) for longitudinal 
equations. These expressions hrve been considered in the light of the study airplanes and the 
results summarized in table 10. Detailed results may be found in app. C. 

The longitudinal phugoid mode can also be approximated by two degrees of freedom, 
u and#. Approximate frequency and damping expressions for the phugoid mode are: 



and 
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A set of approximate expressions for the lateral-directional moues is given in ref. 36. 
The Dutch roll, rolling covergence, and spiral modes are given as 



(180) 



^or a less complicated expression 
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which, together with equation (180), gives good results isi some cases) 
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The merits of using equations (180) and (182) are discussed in detail in app. C and 
summarized in table 10. The complexity of equations (183) and (184) makes their use 
almost ineffectual when compared with the use of the small perturbation program. These 
expressions for 1/T r and 1/T S were not used in this study. 

Similarity between the rigid and equivalent elastic mathematical models implies that 
the use of the approximate expressions could also be used for equivalent elastic cases with 
some minor redefinitions of terms involved. 

Some studies have been made for the completely elastic airplane from the viewpoint of 
approximate transfer functions. In one study (ref. 39), some approximate frequency and 
damping expressions were obtained for rigid-body modes with one and two elastic modes 
for three dissimilar configurations. It is apparent from the study that the form and accuracy 
of approximate expressions are sensitive to configuration, number of elastic modes being 



considered, and dynamic pressure. The treatment of special problems for known significant 
isolated elastic effects appears to be in the approximation category. The availability, speed, 
and versatility of digital computer techniques tends to preclude the use of approximate 
expressions for solving general problems in elastic airplane dynamics. 

8.3.5 Discussion of results of the characteristic equation methods,— The most impor- 
tant conclusions arrived at from the longitudinal dynamic analyses were that the stability 
characteristics are more sensitive to aerodynamic derivative accuracy than to elastic effects 
for the study airplane cases. In addition, the effects of elasticity were relatively small; this is 
illustrated in fig. 19 for the 707-320B short period frequency and damping characteristics. 

Figure 20 shows that the addition of dynamically participating elastic modes to the 
SST configurations has much less effect on the short period frequency (2x2 versus 22 x 22 
modes) than does the static-elastic type of correction (3 x 3, rigid versus equivalent elastic). 
The effect of dynamic pressure is also illustrated in fig. 20 at M = 2.7. It appears that an 
increase in q at "that condition has an overall stiffening effect, as observed by comparing the 
elastic increments between the comparable models. 

A particularly disturbing quality of the data in fig. 20 is the lack of consistency in the 
effects of elasticity. For example, the truncated, completely elastic data show increases, 
decreases, and no changes in the frequency. In addition, the change between rigid and 
equivalent elastic frequency shows increases in frequency for 42“ sweep and decreases at 
72® sweep. This precludes guessing or making any general statements as to the overall 
effects of elasticity. This is even more evident as illustrated in fig. 21 where the undamped 
natural frequency is presented for the SST. In each case shown, the effect of adding elastic 
degrees of freedom (generalized coordinates) is illustrated. 

In general, the frequency increased when the two lowest frequency elastic modes were 
added, then decreased when the next sets of two were added out to eight total modes. From 
there on adding modes had rather unpredictable effects, except that in all cases the fre- 
quency tended to approach a constant value as more modes were added beyond 12. The 
apparent inconsistency is that the constant value is not always less than or more than the 
rigid 2x2. 

The study also showed that adding many elastic degrees of freedom consistently de- 
creased the damping of the short period mode for the study airplanes. 

The adequacy of a particular mathematical model (structural or aerodynamic) for the 
longitudinal dynamics (table 10) would be an important consideration of a handling- 
qualities study. 


159 



Time to damp to half amplitude, Ty 2 , sec Persod sp , P sp , 



Mach number 

FIGURE 19.- 707-320B LONGITUDINAL SHORT PERIOD RIGID VS. EQUIVALENT ELASTIC 


Short period frequency, o> sp , rad /sec 



FIGURE 20. - RIGID, EQUIVALENT ELASTIC AND COMPLETELY ELASTIC SHORT 

FREQUENCY - 42° AND 72° SST 


161 




Short period undamped natural frequency, ft> 
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FIGURE 21- SST, A LE =42° & 72° VARIATION OF SHORT PERIOD FREQUENCY 
WITH NUMBER OF MODES 


For the study airplanes at the flight conditions analyzed, the elastic effects on the 
Dutch roll frequency are quite small. This is illustrated in fig. 22 where the undamped 
natural frequency and damping are shown for two 707-320B flight conditions for various 
numbers of elastic modes. The effects of residual flexibility are also shown. A static analysis 
predicts the period of this mode accurately enough for stability and control purposes. The 
damping of the Dutch roll mode decreases with the addition of the first few elastic modes, 
but then increases very slightly as more elastic degrees of freedom are added. A static-elastic 
analysis would appear to predict the damping with sufficient accuracy for this 
configuration. 

The various truncated, completely elastic airplane models gave good correlation with 
flight test data for the 707-320B for the Dutch roll mode. Figure 23 shows the correlation 
of the damping. Also shown are the poor results obtained from the equivalent elastic 
handbook method. The method has certain deficiencies, which are discussed in apps. B and 
C. 


For many cases the variations in dynamic characteristics for the rigid airplane due to 
the use of different methods for calculating the stability derivatives are as large as any elastic 
effects (for all modes). This points to the fact that a sophisticated, completely elastic 
airplane mathematical model is only as good as the basic rigid stability derivatives. (Table 10 
summarizes the effects of elasticity on the dynamics for all configurations.) Therefore, a 
need exists for an accurate analytical approach to generating all stability derivatives in 
conjunction with lifting surface and lifting line aerodynamic theories. The longitudinal rate 
derivatives (C m . , C m , Cl , etc.) and all lateral-directional stability derivatives need to be 
mechanized. 


8.4 Dynamic Stability Characteristics 
by Time History Solutions 

8.4.1 Applicability of time history solutions.— There are today several practical cases 
where nonlinearities in the equations of motion (dynamic or aerodynamic) are large enough 
that they cannot be neglected. It has been common practice in such cases to base judgment 
of stability behavior on time history solutions of the equations of motion. A time history is 
a set of data that describes airplane motions as a function of real time, i.e., {X} = {X(t)}. 

Time histories have the advantage of providing a clear physical picture of the motion of 
the airplane. In addition, they have the merit of allowing a direct comparison of analytical 
with experimental data. 




FIGURE 22. - 707-320B DUTCH ROLL FREQUENCY AND DAMPING FOR 
COMPLETEL Y ELASTIC AIRPLANE 


10 


o> 

£ 

a 

6 

a 


8 -5 


o 

±j> 

o' 


a. 

E 

co 

O 


-10 


-15 


*20 


-25 










■ 














■ 






[ 

Stable 

A 



i 




■ 


— 





j 


31 



■ 


. 10 6 

75 m 


f 

1 

- 





r 



■ 

L 


i 

i 

j 



3050 n 

i 







< 

r< 

m 

■ 

■ 


o 









i 

s 

/ 




V 

L 







■ 

■ 

HHBHH 



rr 

h = 3050 m 

" 










■ 

■ 

■ 

■ 

■ 


in 

Jnstabl 

- 

D 






■ 

■ 

■ 

■ 



L- 







■ 

■ 

■ 

■ 




u 







■ 

■ 

■ 

■ 







t 

r 



■ 

■ 

■ 

■ 







* 




■ 

■ 

3 

■ 



! 





i 



| 

■ 

B 

■ 



| 

Symb 

ol IVU 

ithod 






! 

■ 

■ 

g 



A Rigid wind tunnel 

[j Equiv elas handbook 

^ Flight test 

Q 14 static elastic modes 

^ 3 modes dynamic, 1 1 modes static ela 

£ 5 elastic modes 

- <3 17 elastic modes 

l 1 ____ I 1 I J 


■ 

■ 

■ 

■ 


■ 


■ 

■ 

■ 

■ 


1 

stic 

■ 

■ 

■ 

i 


1 

S 

i 


■ 




■ ■ 


0.4 0.5 0.6 0.7 

Mach number 


0.8 


0.9 


1.0 


FIGURE 23.~ 707-320B DUTCH ROLL DAMPING, ANALYSES RESULTS VS. FLIGHT TEST DATA 
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Time histories can be generated by integrating with respect to time the complete 
airplane equations of motion or, for that matter, any of the equations of motion shown in 
Sec. 5. The integration technique may vary, but the approach is generally the same for any 
type of computer. The equations must be trimmed (equilibrated) either separately from or 
in conjunction with the problem to be solved, i.e., the solutions {Xj} of the algebraic 
steady-state equations must be obtained and used as initial conditions. The program is 
started with t = 0. At some time when t 0 ^ 0, a disturbance, {AX}, is introduced and the 
response, {X(t)}, calculated for t 0 < t< tj, where t j - 1 0 is usually a time interval sufficiently 
long to establish stability behavior but not so long as to involve mass or other changes that 
would significantly affect assumptions made in deriving the airplane equations of motion. 

In this fashion, it is possible to determine stability behavior by observation, i.e., by 
judging the behavior of the variables of the resulting time history. The stability criteria 
associated with time histories have been stated in Sec. 6. 

An important observation must be made. For nonlinear equations of motion (see app. 
A), several different cases involving disturbances different in both kind and magnitude must 
be run to obtain sufficient information to establish the stability behavior. The reason for 
this must be found in the property of nonlinear differential equations, i.e., that their 
response behavior can be a function of the initial disturbance. 

For the linearized, uncoupled, small perturbation equations of motion (app. A), only 
one arbitrary disturbance is required for each mode (longitudinal or lateral-directional). 
Linearity implies that the response behavior is independent of the size or type of disturb- 
ance in that mode. However, time history generation for the linear, small perturbation 
equations is not necessarily the most efficient approach to stability analysis. 

The major advantage of the time history (integration) approach is that it is in terms of 
real time. The analyst has more physical feel for the problem, since he observes motions 
similar to those which the airplane would experience in flight under the same conditions. 
Most of the disadvantages of the time history method are not really pertinent to the 
problem of stability behavior. Instead, they are of an economic nature, involving such things 
as acquisition, unkeep, and availability of hardware and facilities; and man-hour expendi- 
tures in programming, data preparation, and reduction. 

8.4.2 Rigid and equivalent elastic airplanes. — The time history technique for rigid and 
equivalent elastic models is essentially that written for the rigid airplane. The mechanized 
solution will subsequently be referred to as the “rigid-body, six-degree-of-freedom pro- 
gram,” even though it is also used for equivalent elastic solutions. The equations of motion 
solved by the program, as described in app. A, are the “equations of arbitrary motion.” 
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For this program these equations can be nonautonomous. Thrust forces and moments 
may be input as explicit functions of time. In addition, aerodynamic forces and moments 

j 

due to controls may be explicit functions of time. The aerodynamic data may be nonlinear 
and aerodynamic cross-coupling may be included, e.g., C^, C m ^, and Cj^. Certain negli- 
gible derivatives such as Cj^and Cy^ have been neglected. 

A description of the particular integration scheme used to solve the equations of 
motion for this program may be found in app. C. 

The rigid-body, six-degree-of-freedom program is capable of analyzing handling- 
qualities problems. Part of the basic program output is the veil .'ty (Up, Vp, Wp) and 
acceleration (Up, Vp, Wp) at the pilot’s station. Also, because engine thrust may be input 
separately for each engine as an explicit function of time, the program has the capability to 
analyze engine-out-type time history solutions. 

8.4.3 Completely elastic airplane. - The time history solutions of the completely 
elastic airplane equations of motion for this study were obtained Uoing a special program- 
ming language called MIMIC. This special technique is documented in ref. 40. The time 
histories are merely the time-dependent analogs of the frequency-dependent equation (174) 
with initial equilibrium conditions to which the perturbations are added along with a dis- 
turbance to excite the system. The scheme is simple. From the equilibrium conditions and 
the disturbance, the accelerations are calculated, e.g., 
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q^q}. {q}> {q}) 

i = 1, 2, . . . , n 


(185) 


these are integrated by making the statements 


and further by 



(186) 


q. = INT^q., q.(0)^ 

= q i (0)+ / q i dt i = l»2, . 


(187) 
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Included as a subroutine in the MIMIC program is a Runge-Kutta numerical integration 
technique that accomplished the integrations, equations (186) and (187). This approach is 
easy for the engineer since the programs are then physically oriented. However, as one might 
expect, the easier the program is to write, for a given problem, the more time is required to 
execute it. 

As was mentioned already, this solution is an analog of the frequency-dependent 
equation (174). Indeed, aside from errors inherent in the numerical integration technique, 
there should be no difference between the MIMIC solutions and those obtained from the 
explicit expression 


([[a]S 2 + [B]S+ [C]]" 1 {q 0 }) (188) 

where indicates the inverse Laplace transformation. 

8.4.4 Discussion of time history solutions.- For the time histories generated for this 
report, the aerodynamic coefficients were of a linear nature only. Also, the lateral- 
directional and longitudinal modes were constrained to be uncoupled. The time history 
methods were used here to obtain a graphic presentation of the motion in response to 
various disturbances. The frequency and damping characteristics are u-entially those given 
by the characteristic equation method. 

An example of the use of the MIMIC program to generate longitudinal time histories is 
shown in fig. 24. The response of the 707-320B to an elevator pulse is shown for three 
mathematical models. The rigid model corresponds to the truncated, completely elastic 
model resulting in the use of a 2 x 2 matrix. The inclusion of 14 dynamic modes in a 
static-elastic manner results in the greatest pitch amplitudes. If four modes are allowed to 
participate dynamically, the result is an effective increase in the damping as reflected in the 
decrease of amplitude. 

Appendix C presents r -e examples of the use of MIMIC and the six-degree-of- 
freedom program to generate . U ue histories. 
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Pitch perturbation, A© , deg A<5g, deg 
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FIGURE 24.- 707-32Q3 PITCH PERTURBATION IN RESPONSE TO ELEVATOR PULSE 





9.0 SOLUTION OF THE COMPLETE PROBLEM 

to 

to 


9.1 Arrange ment of a Computation System 

Previous discussions have dealt with individual parts of the overall problem of deter- 
mining the stability characteristics of an elastic airplane. The pertinent methods can now be 
arranged to form a general calculation procedure. Certain preparatory calculations and 
inputs must also be included as part of the procedure. 

To handle elastic airplane problems completely and in sufficient detail for most appli- 
cations, computer mechanization of the methods is necessary. On this basis, a flow diagram 
of recommended computer program system is given in fig. 25. This system is now under 
development. 

The airplane definition section computes the characteristics of an airplane in terms of 
aerodynamic and structural influence coefficients, free vibration normal modes, and inertias. 
This section will not accept empirical data as input. The airplane stability evaluation section 
evaluates the stability characteristics of the airplane, using the results from the airplane 
definition section either aloi.e or in conjunction with empirical data which would apply 
corrections to or replace the computed results from the airplane definition section. 

The airplane definition section consists of four computer programs identified as: geom- 
etry definition (GD), aerodynamic influence coefficients (AIC), structural influence coeffi- 
cients (SIC), and normal modes (NM). The airplane stability evaluation section consists of 
three computer programs identified as: stability derivatives and static stability (SD&SS), 
characteristic equation rooting (CER), and time histories (TH). Specifications and com- 
ments oi each of the program elements are given in following sections. Element specifica- 
tions were arrived at after careful consideration of the detailed results and developments 
contained in the appendixes as well as of information in this document. 


9.2 Use of the System 

The system as presented above could be operated in a number of ways to suit the 
particular needs of the user. It has the capability of accepting very detailed inputs as well as 
empirical and test data to obtain very accurate answers. It can also be used in a less complex 
manner to find preliminary-design-type data. The major use classifications can be stated as: 

1. Preliminary design; 

2. Configuration development; 

3. End product development 



Geometry 

definition 


Mach no. 


Input: Description of airplane 

surface; panel density, control 
surface definition 
Output: Panel areas; panel cen- 

troid and control point locations; 
thickness, camber and dihedral 
slopes; panel corner point loca- 
tions; control surface panels 


Output: Symmetric and/or 

non symmetric AlC's for 
thickness; incidence 


Input: Elastic axis data or 

flexibility matrix; mass dis- 
tribution 

Output: Mass distribution, 

inertias, c g location, com- 
puted or interpolated deflec- 
tion and slope SIC's 


Airplane 

definition 

section 


Input: Frequency cutoff or 

number of modes 
Output: Slope and displace- 

ment mode shapes 


Direct analytical 
form 

Empirically 

corrected 


Input: Dynamic pressure 

Output: Stability derivatives for 

rigid, equivalent elastic and residual 
flexibility, also, static stability data 
Input: Pressure and/or force data 

Output: Corrected stability deriva- 

tives and static stability data 


Airplane 

stability 

evaluation section 


Small Large 
pert. pert. 


A/P and structural 
damping and frequency^ 


Output: Motion variables versus time 


FIGURE 25.- COMPUTING PROGRAM SYSTEM-FLOW DIAGRAM 



9.2.1 Preliminary design use. — This category is one for which only relatively rough 
input data are usually available. Structural characteristics and mass distribution are very 
approximate, and no wind tunnel or other aerodynamic test data are available. Under these 
circumstances, coarser airplane paneling would ordinarily be used than for configuration or 
end product development. It would also be usual to concentrate on finding stability deriva- 
tives and static stability characteristics from the SD&SS program. Use would not ordinarily 
be made of the NM and TH parts of the system nor, possibly, of the CER program. 

Therefore, by using the GD, AIC, SIC, SD&SS, and CER elements, preliminary-design- 
type answers could be obtained on the elastic effects. Comparison of the computed elastic 
and rigid stability derivatives and static stability would provide a good first look at the role 
th't flexibility plays for a particular configuration. It should be remembered, however, that 
the specific rigid and elastic values may be in considerable error, although better than 
obtainable in the past, even when the elastic-to-rigid ratios or increments are adequate. 

Use of the system in this manner allows relatively easy calculation of effects due to 
configuration changes, loading, flight regime changes, etc., for both the rigid and the elastic 
airplane. 

9.2.2 Configuration development use. — For this use, input data for the system are 
more accurate and complete than for preliminary design. By this time, detailed structural 
analyses have generally been accomplished, weights determined, and wind tunnel tests run 
on the basic configuration. Usually aerodynamic force data are available, but not pressure 
distribution data. 

The GD element is now used to give a more accurate description of the airplane. 
Denser paneling is specified. Details such as dihedral and detailed geometry are included. 
The AIC program, using the better description, can now calculate more accurate aero- 
dynamic data. Items such as nacelle and wing-body-tail interference, which may have been 
neglected in the preliminary design studies, are now included. 

In the SIC element the early matrix obtained with an approximate beam analysis is 
replaced with a more exact beam analysis or with an externally developed matrix. If the 
externally developed matrix is based on different paneling than the AIC’s, an interpolation 
routine in the SIC program changes the matrix to be consistent with the aerodynamic 
paneling. 

It may be elected at this point to include the structural dynamic effect on the stability 
characteristics. This is initiated by exercising the NM program to provide data for inputing 
programs in the airplane evaluation section of the system. 
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In the stability evaluation section, the user has the opportunity to supplement or 
replace previously calculated data with vest, handbook, or empirical data. A modest amount 
of adjustment is normally desirable for cases where viscous effects and nonlinearities are 
small. The adjustment would usually take the form of applying calculated elastic-to-rigid 
ratios to wind-tunnel-obtained stability derivatives. For example, 


m a 

, Elastic 

m a . ~ C m 

e< l- el - “Rigid 


• c 


m a . . 

Computed Rigid 


W.T. 


Where viscous effects and/or nonlinearities are important, more extensive adjustment 
should be made. 

The SD&SS program then calculates the derivatives and static stability using the 
previously calculated data plus desired empirical inputs. 

To obtain the dynamic stability characteristics, one chooses either the CER or TH 
program. The CER program is faster but is more restricted, since it is based on the linear, 
small perturbation equations. The TH program would be required for cases where large 
perturbations describe the motion. It can also be used for the small perturbation case, but 
would not give any more information than could be obtained from the CER program. 

The structural dynamic effects can be included in either program by inputing the NM 
program results. Comparisons of rigid, equivalent elastic, -id completely elastic airplane 
representations are then possible. 

9.2.3 End product development. - The methods for this area of use are an extension 
of those of the previous section. More and more test data would be fed into the system so 
that ultimately the system might be considered as a vehicle for calculating perturbations to 
the test information. In particular, pressure distribution data, local separations, complete 
surface stall, viscous wake effects, and other effects would be empirically included. The 
extent of this type of investigation is usually limited by the time and experience of the user. 
The main value of this approach is to investigate problem areas that may come to light late 
in the development of a configuration or in flight test. 
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9.3 System Element Descriptions 

» 

fc 

9.3.1 Geometry definition (GD). — This program accepts as input the basic geometric 
description of the airplane’s surface either in its cruise condition or in its unloaded condi- 
tion (jig shape). In addition, the program accepts an aerodynamic paneling density selection. 
Constraints placed on the paneling selection reflect, primarily, the requirements of the 
aerodynamic representation. The program will have the option of letting the program select 
the paneling based on a selected density or of operating with a user-selected paneling. The 
number of aerodynamic panels will be open-ended. 

The program will compute thickness, incidence, and dihedral slope at each panel con- 
trol point and panel centroid, and body surface paneling required for wing-body-tail inter- 
ference flow on a cylinder of mean body radius. It will also compute the coordinates of the 
comers, centroids, and control points of the aerodynamic panels as well as their areas. 

9.3.2 Aerodynamic influence coefficients (AIC). — The surface of the airplane is 
divided into panels in the geometry definition (GD) program. The function of the AIC 
program is to compute the change in the pressure force coefficient at each surface panel due 
to a unit change in inclination to the flow at each panel. The AIC’s will include the effects 
of a slender body, thin wing, and tail as well as wing-body-tail interference. The approach is 
essentially that of ref. 35. The computation is Mach number dependent and will handle 
subsonic and supersonic flows up to Mach 5. Transonic flow in the range of about Mach 0.9 
to 1.2 is not handled rigorously, so calculations should not be made in the regime. Methods 
are satisfactory for altitudes up to 30 000 meters (about 100 000 ft). 

The AIC’s for changes in surface panel inclinations to the flow which are symmetric 
with respect to the airplane’s plane of symmetry are computed separately from those which 
are nonsymmetric. The method of ref. 35 consists of representing the perturbation of a 
uniform irrotational flow (due to the presence of the airplane) by line singularities (at inner 
surfaces of the wing and tail). 

Body thickness, camber, and incidence are represented by line sources and line 
doublets on the body axis. The wing and tail surface thickness slopes are represented by 
surface distributions of sources with linearly varying strengths. The effects of wing and tail 
surface incidence and camber, as well as body interference on the wing and tail surfaces, are 
represented by vorticity distributions. 

Unsteady aerodynamic effects can be accounted for in a manner consistent with the 
requirements of airplane stability evaluation. The method suggested by Miles (ref. 28) for 
reduced frequencies that are less than unity will be used to reduce the aerodynamics of 
unsteady flow to one of steady flow. 



An aerodynamic influence coefficient theory for nonsymmetrical flow past wing- 
body-tail combinations has not previously been developed. This development is now being 
carried out. The techniques involved in the development and theoretical justification are 
outlined by Van Dyke (ref. 41), Ashley and Landahl (ref. 27), and Chester (ref. 29). 

9.3.3 Structural influence coefficients (SIC).- This program computes the elements of 
four flexibility matrices that have the following properties: 

• Matrix (1) gives the displacement of each panel control point due to a unit load at 
each panel centroid. 

• Matrix (2) gives the displacement of each panel centroid due to a unit load at each 
panel centroid. 

• Matrix (3) gives the rotation about the y body axis (z body axis for vertical tail and 
body) at each panel control point due to a unit load at each panel centroid. 

• Matrix (4) gives the rotation about the x body axis due to a unit load at each panel 
centroid. 

• All computed flexibility matrices are such that a self- equilibrating system of loads 
applied at panel centroids does not give rise to a displacement of the airplane’s 
center of gravity or a rigid-body rotation of the airplane about its center of gravity. 

The program will also compute the location of the airplane’s center of gravity and the 
components of its inertia in the body axis system. 

The program accepts as input: 

• Distribution of mass based on paneling computed in GD program. 

• Distribution of bending and torsional rigidities along an elastic axis system. 

• Geometric description of the elastic axis system, 
or 

• Elements of a flexibility or a stiffness matrix computed exterior to this system using 
an arbitrary selection of structural modes. 

If elastic axis data are input to the program, the flexibility matrices are computed from 
beam theory for an arbitrary number of beams. The flexibility matrices are computed by an 
interpolation method if a flexibility or a stiffness matrix is input to the program. The 
program is such that an arbitrary number of mass distributions may be combined using 
multiplying factors. 

9.3.4 Structural normal modes (NM). - This program is essentially an eigenvalue pro- 
gram. The free vibration of the structure is given by the eigenvalue problem 

(M ♦<?[*]) <♦„>-" < l89 > 


176 



1 The eignevalues, a> n , are the natural frequencies. The eigenvectors, {<p n }, when normal- 
ized, are the free vibration normal mode shapes. They give the deflected shapes of the 
airplane at the resonant natural frequencies, 0 J n . Slope normal mode shapes can be calcu- 
lated from the deflection shapes. These results are combined into three mode shape 
matrices. One relates elastic deflections at the panel control points to generalized coordi- 
nates (elastic degrees of freedom), another relates angle-of-incidence changes at the panel 
control point to elastic degrees of freedom, and the final one relates dihedral changes at the 
panel control points to elastic degrees of freedom. 

9.3.5 Stability derivative and static stability (SD&SS). - This program accepts the 
following as input from the airplane definition section: paneling geometry, normal mode 
shapes, reference flight condition attitude (0 {, 9 1 ), Ma:h number, dynamic pressure and 
load factor, flexibility matrices, aerodynamic influence coefficients, mass distribution, and 
control surface ‘deflection angles. It will compute all significant stability derivatives based on 
aerodynamic influence coefficients; as well as control surface trim angles for the reference 
flight condition, airplane shape in reference flight condition, stick-speed stability, elevator 
angles per g, neutral point, and maneuver point. 

The matrix equations for the stability derivatives formulated for residual flexibility are 
fundamental to this program. These were presented and discussed in Sec. 7. Recall that the 
matrices [AjJ, [A 2 ], [A 3 ], [A 4 ], and [A 5 ] appearing in the equations of motion intro- 
duced the aerodynamic forces. These matrices, when premultiplied by the matrices [p]T 
[Bi]"l or [*lJ T [Bj]' 1 , result in the airplane stability derivatives for the residual- 
flexibility formulation. If the matrix ['Kid in the equation of motion is replaced by a 
matrix whose elements are all zero, the stability derivatives reduce to the equivalent elastic 
stability derivatives. These are derived in app. B from the flow boundary condition at the 
airplane’s surface. If, in addition, the flexibility matrix [C] is set equal to a zero matrix, 
then the stability derivatives are the rigid airplane stability derivatives. Thus, the stability 
derivatives in the form contained in equation (47) are the most general that may be chosen, 
since rigid and equivalent elastic derivatives may be readily obtained from them. 

Inclusion of a corrector matrix technique in the SD&SS program based on the analysis 
presented in app. B is recommended. That technique utilizes a diagonal corrector matrix 
that corrects the aerodynamic influence coefficient matrix using wind tunnel pressure model 
data. A correction is made at each test condition, e.g., angles of incidence (a ,(3) and Mach 
number. The SD&SS program will accept pressure model data as input. The pressure data 
cannot be constrained to correspond to the aerodynamic panel centroid Hection as they are 
in the GD program. Therefore, the user is required to setect those aerodynamic panel 
centroids which are to be corrected by each pressure data point. The SD&SS program 
partitions the aerodynamic influence coefficient matrix in accordance with the user’s 



selection and performs the corrections on the basis of the partitioned matrix. The program 
will compute the stability derivatives and static stability results for each test condition for 
which pressure data is input. The results are then tabulated in a form acceptable as input to 
the TH program. 

The program will also accept as input stability derivative data obtained from wind 
tunnel and flight testing in tabular form. At the selection of the user, ths data would be 
used in combination with corrected or uncorrected aerodynamic and structural influence 
coefficients to obtain stability derivatives that are corrected for elastic effects. 

The corrections that can be computed are as follows: 

LE/LR = ratio of lift, elastic to rigid 
DE/DR = ratio of drag, elastic to rigid 
YE/YR = ratio of side force, elastic to rigid 
ME/MR = ratio of pitching moment, elastic to rigid 
RE/RR = ratio of rolling moment, elastic to rigid 
NE/NR = ratio of yawing moment, elastic to rigid 
Aac = change in aerodynamic center, elastic from rigid 

The program will be such that these elastic-to-rigid corrections can be computed for 
complete configurations to correct rigid stability derivatives and control derivatives obtained 
exterior to the program. The program can also compute the elastic-to-rigid corrections for 
contributions due to the components of a complete configuration, i.e., wing, body, and tail. 

9.3.6 Characteristic Equation Rooting (CER).— This program operates only on the 
basis of the small perturbation equations of motion. The stability derivatives will, therefore, 
be constants and the equations of motion will be linear, ordinary differential equations with 
constant coefficients. 

The program will accept as input the coefficients of the motion variables appearing in 
the equations of motion. It will combine the coefficients in the appropriate form for the 
rooting method. The method should be programmed in open-ended form so as to accept as 
many degrees of freedom as desired. 

The program computes: 

• times and number of cycles to damp to half and to one-tenth amplitude, 

• frequency and period of modes, 

• undamped natural frequency of modes, 

• damping ratios of modes, 

• phase and amplitude of model coupling terms, e.g., argument, and 

• (0/j3) and magnitude 
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9.3.7 Time Histories (TH). ~ This program integrates the large perturbation equations 
of motion of the airplane by a variable, step-size Runge-Kutta method. The computed 
motion variables are as follows: • 

. • • 2 2 2 ^ • 

u, v, w, u, v, w, V c = (u + v + w ) , p, q, r, p, q, r, 0, ip, 0, n z and h 


The program will be such that the stability derivatives are input as constants or as 
variables. When the stability derivatives are variables, their values as functions of the motion 
variables, will be input from tables and interpolated or extrapolated linearly. The interpola- 
tion method can handle three independent variables; where more than three are required, 
superposition will be used to obtain the value of the stability derivative. 

The forces .and moments from engine thrust will include the effect of engine location 
and attitude. Inertial effects from rotating engine parts can also be included. 

The mass and components of inertia are constant. The airplane will be considered to be 
flying over a flat, nonrotating earth. 
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10. CONCLUSIONS AND RECOMMENDATIONS 


The conclusions and recommendations listed below have been formulated following an 
examination of the results of work reported both in this document and in the appendixes. 

1. In general, reliable pndictions of elastic airplane longitudinal stability and control 
characteristics can be made using current state-of-the-art theoretical methods. These 
methods are applicable and practical for preliminary design purposes. 

2. The general equations of motion, along with either the large perturbation or small 
perturbation equations of motion, were found to be suitable for determining refer- 
ence motion and stability characteristics of an airplane when it is disturbed from the 
reference motion. The small perturbation equations are applicable to nearly linear 
systems. The large perturbation equations are necessary in examining nonlinear 
systems and large disturbances from the reference motion. 

3. The mathematical formulation of stability criter. is the same for rigid, equivalent 
elastic, and completely elastic airplanes. 

4. Static stability is usually, but not always, a prerequisite for dynamic stability and 
good handling qualities. 

5. Energy decay methods and Lyapunov stability theory have potential stability cri- 
teria application. However, more research is needed to establish their practical 
application. 

6. Handbook methods in some cases give low-quality results for estimation of rigid air- 
plane stability derivatives. In addition, taken alone they cannot adequately predict 
elastic effects. 

7. The aerodynamic influence coefficient method using lifting surface theory gives 
generally acceptable results for estimation of rigid and equivalent elastic longitudinal 
derivatives. The method is applicable to lateral-directional problems but has not 
been mechanized as yet; this should be done. 

This method, when arranged for accepting some empirical and test data, was 
judged to be the best available approach for predicting airplane stability derivatives. 
Incorporation of leading edge suction, more accurate shed vortex field representa- 
tion, and other improvements are possible and should be developed where greater 
accuracy is desired than would be obtained with the basic method. 

8. If wind tunnel data are available, the method of (7) above can be used to provide 
more accurate predictions of equivalent elastic effects. This is done by computing an 
elastic-to-rigid ratio or increment, which is then applied to the wind tunnel value. 
For example, 


m 


i - ^Elasti c | 

'm a ~ C 

^■ El - “Rigid 



Rigid 


Computed 


W.T. 



inis approach is recommended as a practical and reasonably accurate way to obtain 
equivalent elastic effects. 

9. The remarks in (6), (7), and (8) above also apply "‘to prediction of airplane static 
stability. 

10. Both static and dynamic stability characteristics were found to be sensitive to 
inaccuracies in estimating the rigid stability derivatives. This uncertainty is usually as 
large as the elastic effects. 

11. Dynamic flexibility effects were usually found to be modest and smaller than the 
equivalent elastic effects. About 20 elastic modes are required to obtain good re- 
sults; using less than this number can lead to significant error if residual flexibility is 
not used. 

12. For many cases the equivalent elastic formulation represents the airplane accurately 
enough so that there is no need to go to the extra complication of the completely 
elastic formulation. Engineering judgment is required to decide on this for dynamic 
stability evaluation of any particular configuration, however. 

13. Approximate formulas for determining damping of longitudinal dynamics were satis- 
factory provided C m< ^ is dominant. However, these methods are unreliable for gen- 
eral use and should be avoided. 

14. Characteristic equation rooting and time history methods are adequate and are 
recommende for use in determining both longitudinal and lateral-directional 
dynamic stability characteristics. 
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